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CHAPTER 1 
INTRODUCTION 


• 


• 


• 


1-1. General. 


a. This 
series 
of manuals, 
entitled 
Designing 
Facilities to Resist Nuclear Weapon Effects, is 
organized as follows: 
TM 5-858-1 
Facilities System Engineering 
TM 5-858-2 
Weapon Effects 
TM 5-858-3 
Structures 
TM 5-858-4 
Shock Isolation Systems 
TM 5-858-5 
Air Entrainment, Fasteners, Pene- 
tration Protection, Hydraulic-Surge 
Protective Devices, EMP Protective 
Devices 
TM 5-858-6 
Hardness Verification 
TM 5-858-7 
Facility Support Systems 
TM 5-858-8 
Illustrative Examples 
A list of references pertinent to each manual is placed 
in 
an 
appendix. 
Additional 
appendixes 
and 
bibliographies are used, as required, for documenta- 
tion 
of 
supporting 
information. 
Pertinent 
bibliographic material is identified in the text with 
the author's name placed in parentheses. Such 
bibliographic material is not necessary for the use of 
this manual; the name and source of publications 
related to the subject of this manual is provided for 
information purposes. 
b. The purpose of this series of manuals is to pro- 
vide guidance to engineers engaged in designing 
facilities that are required to resist nuclear weapon 
effects. It has been written for systems, structural, 
mechanical, electrical, and test engineers possessing 
state-of-the-art 
expertise 
in 
their 
respective 
disciplines, but having little knowledge of nuclear 
weapon effects on facilities. While it is applicable as 
general design guidelines to all Corps of Engineers 
specialists who participate in designing permanent 
military facilities, it has been written and organized 
on the assumption a systems-engineering group will 
coordinate design of the facilities. 
c. Technical Manual 5-858 addresses only the 
designing of hardened facilities; other techniques to 
achieve survival capacity against nuclear weapon at- 
tacks 
are 
deception, 
duplication, 
dispersion, 
nomadization, reconstitution, and active defense. A 
facility is said to be hardened if it has been designed 
to directly resist and mitigate the weapon effects. 
Most of the hardening requirements are allocated to 
the subsidiary facilities, which house, support, and 
protect the 
prime 
mission 
materiel/personnel 
(PMMP). This manual is applicable to permanent 


facilities, such as those associated with weapon, 
systems, 
materiel 
stockpiles, 
command· centers, 
manufacturing centers, and communications centers. 
d. The nuclear weapon threats considered are listed 
below. Biological, chemical, and con~entibnalweapon 
attacks are not considered. 
- Weapons aimed at the facility 'itself or at nearby 
targets 
A range 
from 
many, 
relatively 
small:yield 
weapons to a single super:yieldweapon 
Weapon yields from tens of kilotons to hundreds 
of megatons 
- Weapon delivery by aerial bombing, air-to- 
surface missile, surface-to-surface missile, or 
satellite-launched vehicle 
- Detonation (burst) of a weapon in the air, at the 
ground surface, or beneath the ground surface 
- Direct-overhead bursts for a deep-buried facility 
- Near-miss bursts for a near-surface facility, pro- 
ducing peak 
over-pressures from 
tens 
to 
thousands of psi at the facility 
e. The designing of facilities resistant to nuclear 
weapon effects is an evolving specialty that uses a 
relatively narrow data base that incorporates both 
random and systematic uncertainties. The range of 
these uncertainties may vary from significant (order 
of 1 to 2 magnitudes) to normal (10% to 100% varia- 
tion from average values). The applicable uncertainty 
value depends on the specific weapon effect or harden- 
ing objective under consideration. Loading uncer- 
tainty is generally more significant than resistance 
uncertainty. Awareness of the appropriate uncer- 
tainty factor is essential not only for system engineer- 
ing trade-offs but in the utilization of available 
analysis or test procedures. Studies and experiments 
are being conducted to improve methodology, to better 
define random uncertainties, and to reduce systematic 
uncertainties. This manual willI be revised as signifi- 
cant improvements occur in either methodology or 
data base. 


1-2. TM 5-858-2: Weapon effects. 


This manual presents the 13 weapon effects that act 
as the facility "loads." Each effect is described by its 
expected value, its random coefficient of variation, 
and its systematic coefficient of variation. Uncertain- 
ties in the weapon, the weapon delivery, the site, and 
other factors are considered. The effects data are 
presented both graphically and analytically. 


1-1 


" 
I , 


• 


• 


• 


• 


TM 5-858-2 


CHAPTER 2 
NUCLEAR-WEAPON EFFECTS 


2-3. Burst elevations. 


to 1,000,000 times more destructive. It is this much 
larger magnitude of potential destructiveness that 
creates disparity between hardening against nuclear 
weapons and hardening against conventional HE 
weapons. 


o 
° 
< -1000 ftIMtO,25 


< 100,000 ft 
< 3500 ftIMtO. 4 
° 


100,000 ft 


3500 ft/Mt°.4 


2-2. Nuclear weapons. 


Nuclear weapons include atomic and thermo-nuclear 
weapons. The explosive energy of an atomic weapon 
results from the fission of uranium and plutonium; 
the yields of atomic weapons are tens of kilotons or 
less. A thermo-nuclear weapon employs both fission 
and fusion reactions to obtain a very large yield. An 
ordinary nuclear weapon of either basic type can be 
modified to produce an "enhanced" weapon, e.g., the 
output of the weapon can be channeled to produce 
relatively large amounts of nuclear radiation, or if the 
weapon is so designed, it may enter the earth prior 
to detonation resulting in a buried or partially buried 
detonation. 


a. The burst elevation of a nuclear weapon can be 
classifed as follows: 
High-altitude 
< 
- 
< 
burst 
< Air burst 
< Low-altitude 
burst 
Contact burst 
:::: 


-1000 ftIMtO.25 < Shallow-buried < 
burst 
Deep-buried 
burst 
Note the following characteristics. Air burst: Fireball 
does not touch the ground. Contact burst: Center of 
gravity of weapon is slightly above the surface; some 
penetration required with true surface burst in which 
C.g. is at ground surface. Deep-buried burst: Com- 
pletely 
contained 
and 
does 
not 
vent 
to 
the 
atmosphere. 
b. State-of-the-art weapons cannot penetrate to suf- 
ficient depths to result in a deep-buried burst, i.e., a 
contained burst. High-altitude and air bursts general- 
ly are not effective burst elevations for the higher 
hardness levels considered in this volume. However, 
there may be threat scenarios that would require the 
consideration of these burst elevations; their effects 
are selectively included in this volume. Low-altitude 
bursts, contact bursts, and shallow-buried bursts pose 
the greatest threat. They are collectively referred to 
in this volume as near-surface bursts. 


2-1. High 
explosive 
and 
nuclear-weapon 
effects. 


a. A 
nuclear 
explosion produces 
13 potentially 
destructive effects, or environments, that must be ac- 
counted for during facility design: 
- Initial nuclear radiation 
- Electromagnetic pulse (EMP) 
- Fireball 
- Thermal radiation 
- Airblast 
- Airblast-induced ground shock 
- Crater-induced ground shock 
- Crater 
- Ejecta/debris impact 
- Ejecta/debris depth 
- Dust 
- Firestorm 
- Residual nuclear radiation 
With the exception of the first and last, the nuclear 
radiations, these effects are also produced by a high- 
explosive (HE) explosion. Even an HE explosion can 
produce the electromagnetic pulse and the mushroom- 
shaped cloud that are commonly identified with 
nuclear explosion. The gross non-nuclear effects pro- 
duced by the two types of explosions are so similar 
that the explosive energy of a nuclear explosion is also 
expressed in terms of an equivalent TNT-charge 
weight. 
b. This equivalent TNT-charge weight is referred 
to as the explosive energy release, or yield, W. A facil- 
ity designed to withstand the effects produced by an 
HE explosion would be quite similar to a facility 
designed to withstand the effects produced by a 
nuclear explosion of equivalent yield. This statement 
is true for peak overpressure larger than a few tens 
of psi, which are the hardness levels addressed in this 
volume. 
c. Since the protective measures required for the 
non-nuclear effects also provide, in general, sufficient 
protection against nuclear radiation, the protective 
measures required for the two facilities would be of 
the same kind-but would differ drastically in degree 
when the effects are produced by weapons of war. HE 
weapons (e.g., conventional bombs) have yields of only 
several tons or less, whereas state-of-the-art nuclear 
weapons range in yield from a few kilotons to a few 
tens of megatons. The potential destructiveness of a 
weapon at a given distance varies between Wl/3 and 
W, depending on the particular weapon effect under 
consideration. Thus, a small state-of-the-art nuclear 
weapon is 10 to 1000 times more destructive than the 
largest HE weapon, and a large nuclear weapon, 1000 
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~ 
2-4. Weapon effect, a random variable. 


a. Parameters. The specification of a weapon effect 
for use during the designing of a facility requires that 
the following basic parameters be quantified: 
- Weapon design 
- Weapon yield 
- Burst elevation 
- Slant range' . 
- Weather and atmospheric conditions 
Geologic factors 
Because of the uncertainty associated with these basic 
parameters, it is necessary to treat them, and thus 
the weapon effects, as random variables. 
b. Definition o( a random variable. 
A random 
variableX is a numerical variable, the specific value 
of which cannot be predicted with certainty. Every 
specific value of the random variable X, denoted by 
x, represents a possible outcome of the random 
phenomenon, under considenition. 


c. Density (unction. The probability density func· 
tion of X,denoted by fx (x), describes the relative 
likelihood of the various outcomes of X. It is defined 
such that, 
f: fx(x) dx = 1.0 


where fx(x) dx represents the likelihood of X being 
between x and x + dx. This relative likelihood or 
relative freql}ency of occurrence is known as the pro· 
bability of x < X.< x + dx, i.e., 


P(x < X < x + dx) = fx(x) dx 
(2-2) 


Also, note the probability that Xl < X < x2 is 


P(xl < X < x2) =J 
X 2 fx(x) dx 
(2-3) 
Xl 


d. Mean and variability o(a random variable. The 
mean (the average, or expected) value of the random 
variable X is 


e. Evaluation o(probabilities. In general, obtain the 
probability of Xl < X < x2 from equation 2-3. For a 
number of special density Junctions, these pro- 
babilities are tabulated in reference books. (see, for 
example, Abramowitz, 1972). As an example, for a 
Gaussian or normal density function 


fx(x) = 
ox~ ~xp [ -0.5(~ - IAXI210~] 
(2-7) 


Letting z denote the number of standard deviations 
above mean, one obtains 


lAx + zOx 


P(z)= 
PrX <lAx + zoxJ=1 ~ exp 
" 
ax 
2rr 
' 
. 


_00 


[-05 (x-"X)2/0~}x ~J:;'i21rr ex+ ~2] du 


(2-8) 


where transformation u = (x - lAx) lax is used. Values 
of this integral for a number of values of z are given 
below. 


z 
Prz) 


0 
0.500 
0.1 
0.540 
0.25 
0.599 
0.50 
0.691 
0.75 
0.773 
1.00 
0.841 
2.50 
0.994 
5.00 
0.999 


The above table may be used to evaluate other pro- 
babilities for the Gausian distribution; for example: 


• 


• 


lAx =J 00x fx(x) dx 
-00 
(2-4) 
P(-z) = 1 - P(z) 


P(lAx + zloX <X< lAX + z2 ax) = P(z2) - P(zl) 


(2-6) 


A useful measure of the randomness, or variability, 
of the random variable X is its coefficient of variation: 


(2-5) 


where ax is the standard deviation of X and is equal 
to 


Ox ~ I[f 


oo 


oo(x - "xl' [X(x) dX]1I2 


2-2 


f 
Accounting (or uncertainties: A single variable. 
The uncertainty associated with a random variable 
is due to two factors: (1) the natural or inherent ran- 
domness of the random variable itself, and (2) errors 
which result from the inability to accurately estimate 
the mean, the standard deviation, or the density func- 
tion of the random variable. The latter may result, 
for example, from lack of sufficient data. To account 
for both uncertainties the random variable X is writ· 
ten as 
• 


However, if an unbiased set of data is used to evaluate 
x, then Nx = 1.0 and 


where X is the random variable as represented by the 
available data and has a mean xand a coefficient of 
variation Ox which are estimated from the data. Nx 
is a random correction factor which is included to 
model the estimation errors. Its mean value Nx 
represents the bias in data, whereas its coefficient of 
variation 
/),x 
is 
a measure 
of the uncertainty 
resulting from estimation errors. As a first approx- 
imation /),x is considered as the error in estimation 
of the mean value x. The mean value of the random 
variable X is evaluated from 


of 
variation 
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Total 
coefficient 
of 
variation 
associated with the functional form of 
X 


where 
Qf 


= I d2 = 
/),2 
"f 
f 
Total 
coefficient 
associated with Yn 


= 102 
+ /),2 
"Yn 
Yn 


Qnm 
= 
Correlation coefficient of Ynand Yrn' 
-1 
~ Qnrn ~ 1 
The subscript "0" in equation 2-16 denotes that 
8f18Yn is to be evaluated at the mean values ofthe 
variables. Implicit in equations 2-15 and 2-16 are the 
assumptions that Q~ «1 and the nonlinearity in 
f near the mean valugs is not large. In this manual, 
both symbols Q or Q2 are called "uncertainty" to 
avoid the repetitive use of the term "uncertainty 
squared." 
h. Distributions. If a particular distribution is de- 
fined, the uncertainty can be calculated in terms of 
the upper and lower bounds assigned to the data. For 
example, suppose that data are presented within 
bounds that are said to include 95% of the data. Ifthe 
data are assumed to have a Gaussian distribution, the 
bounds correspond to the ± two-standard deviation 
limits (see e above). Procedures for a number of 
specific distributions are described in the following 
paragraphs. 
(l) Normal: The coefficient of variation of a nor- 
mal distribution is 


(2-9) 


(2-10) 


(2-13) 


(2-11) 
I-lx = x 


i.e., for an unbiased set of data, xis the best estimate 
of the actual mean. The total coefficient of variation 
of X, designated by Qx, is obtained from 


Q 
- 
r 02 + 
.2 
(2 12) 
X -...,; 
X 
Ll X 
- 


Here, Ox models the natural randomness of X, 
whereas /),x represents the uncertainty arising from 
errors in estimation. As an example, when a set of 
n unbiased samples are available, it is usually as- 
sumed that 


• 


• 


g. Accounting for uncertainties: Multiple variables. 
If X itself is a function of several random variables, 
1(R - 1) 
Ox = n 
R + 1 
(2-17) 


(2-14) 


then using first-order linear approximations (Ang- 
Cornell, 1974) 


where R = L2!L1' and L1 and L2 are the lower and 
upper bounds corresponding to n standard deviations 
on either side of the mean. 
(2) Log-normal: The coefficient of variation of a 
log-normal distribution is 


and 


(2-15) 


[ 
[ 
UnR)2 ] 
]1/2 
a = 
exp -- 
- 1 
x 
4n2 
(2-18) 


(2-19) 


where R = k2, k = x!L 1 or k =L2/X, and L1 and L2 
are the lower and upper bounds of the distribution. 
The quantity n is that number of standard deviations 
of the bounds from that log mean. 
(3) Uniform: The coefficient of variation of a 
uniform distribution is 


1~ 
V3~~ 
a = 
x 


N Y2 
2 


L 
_n2 (8~ 
) 
Q~ 


n=1 
X 
non 


Q2 = Q2 + 
X 
f 
• 
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where L1 and Lz are the lower and upper bounds of 
the data. 
(4) Log-Uniform: The coefficient of variation of a 
log-uniform distribution is 


d· =..~ (lIn Lz )( Lz + L1) -11 
112 
(2-201 


x 
( 
2 
L1 
Lz - L1 
~ 


where L1 and Lz are the lower and upper bounds of 
the data. 
i. Uncertainty in nuclear weapon effects. Uncertain· 


. "r-' 


ties in nuclear weapon effects arise from yield, 
delivery accuracy, weapon capability (air bursts, 
penetrating 
bursts), 
as 
well 
as 
from 
imperfect 
knowledge of specific effects for known yield, range, 
etc. In the following chapters an attempt is made to 
identify the level of uncertainty of our knowledge of 
weapon effects. The uncertainties relative to weapon 
yield, range, and burst conditions are defined, but 
since they are considered a part of threat definition 
they are not quantified; however, provision is made 
for incorporating uncertainties in these parameters 
from threat definition applicable to specific problems. 
• 


• 


• 


IS 


spectrum and the attenuation characterization ofthe 
standard atmosphere. 
c. Personnel. Personnel injury resulting from radia- 
tion exposure is complex and numerous (Glasstone, 
1977). If the type of weapon treat is known and Ifthe 
slant range is less than 5000 yd (4572 m), it is possi, 
ble to quickly estimate the neutron does to person- 
nel regardless of the spectral content of the pulse'. For 
low-altitude atmosphere bursts, where personnel are 
located on or near the earth's surface, the neutron 
dose absorbed by human tissue is (DNA, 1972). 


.:: 


(3-1) 
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where 
W 
R 
FN 


Weapon yield, kt 
Slant range, yd (m) 
Normalized 
neutron 
dose, 
yd2 
radlkt[m2 (Jlkg)lkt] 
HI 
Burst height correction factor 
(1) The mean value of equation 3-1 is calculated 
by referring to figures 3-2 through 3-5 (which show 
the relationship between FNand R), and to figure 3-6 
(which shows the relationship between HI and the 
height of burst) for the standard division (DNA, 1972) 
of we'apon types, I through VIII. 
(2) The uncertainty associated with equation 3-1 


CHAPTER 3 
INITIAL NUCLEAR RADIATION 


3-1. Gross phenomenology. 


a. A nuclear explosion is accompanied by the emis- 
sion of nuclear radiations consisting of gamma rays, 
neutrons, beta particles, and alpha particles. The 
range of beta and alpha particles is so short that they 
can be neglected in the estimate of the initial radia- 
tion. The initial nuclear radiation is considered to be 
composed of neutron and gamma ray radiation from 
both initial effects (t < 1 minute) and secondary in- 
teractions. Because induced radioactivity (a threat to 
personnel for many hours) begins during t < 1 min, 
it is included in this chapter. 
b. The basic variables affecting the destructiveness 
of initial nuclear radiation are: 
- Weapon design 
- Weapon yield 
Burst elevation 
Slant range 
Atmospheric conditions 
The destructiveness of the initial nuclear radiation 
will be enhanced by any of the following: Larger 
weapon yield, higher burst, smaller slant range. 
c. The primary direct effect of nuclear radiation is 
an antipersonnel effect. Direct nuclear radiation ef- 
fects on equipment and materials are less significant, 
but certain detector materials and electronic com- 
ponents can be damaged. The designer must either 
design for these effect.,; (by shielding, for example) or 
he must preclude them (by burial, for example). The 
two significant aspects of initial nuclear radiation, 
neutron radiation and gamma ray radiation, are con- 
sidered separately. 


• 


• 


'I~ ." " 


3-3. Gamma radiation. 


a. Prompt and delayed. Gamma radiation can be 
divided into two components: Prompt (within 1 or 2 
shakes, a shake being equal to 10-8 sec), and delayed 
(up to 1 min). Electronic equipment is sensitive to the 
energy spectrum of the radiation-only prompt dosage 
is calculated; personnel are influenced primarily by 


where Q[ is the uncertainty of the functional form of 
equation 3-1, Qw is the uncertainty of the yield, etc. 
Consider that Q[ :::: 0.2, QFN = 0.30; using equation 
2-19 and the assumption that figures 3-2 through· 
3-5 
represent 
mean 
values 
of data 
uniformly 
distributed with upper and lower bounds within 
±25%'ofthe mean (from DNA, 1972) andQHI :::: 0.3. 
Thus, equation 3-2 reduces to 


• 


3-2. Neutron radiation. 


a. Initial. The weapon type will have an influence on 
the energy spectrum of the source radiation. The loca- 
tion of the weapon source (and the target) will also 
affect the spectrum. (Specific situations relating these 
parameters are too numerous and specialized for 
presentation here; however, such information can be 
found in the report prepared by SAl, 1971). Only in- 
itial neutron radiation effects (occurring within 1 
minute after detonation) are considered in this 
chapter; residual neutron effects (after 1 minute) are 
discussed in chapter 15. 
b. Equipment: The vulnerabilities of electronic 
equipment are expressed in terms of fluence spectra 
(i.e., incident neutrons/in. 2 or neutrons/cm2 at various 
energy levels expressed in MeV). Typically, fluence 
spectra are presented in charts of the kind shown in 
figure 3-1, which are based on the weapon (source) 


Q2 
= 
0 22 + Q2 
+ 4 Q2 
DN 
. 
W 
R 
(3-3) 


• 


2000 
1600 


yd 


800 
1200 


SLANT RANGE, 


Figure 3-1. 
Typical Neutron Fluence Incident on a Receiver 
Located on or near the Surface of the Earth (DNA, 
1972) 
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Figure 3-2. 
Neutron Dose F.v as a Function of Slant Range from a l-kt Surface Burst, Weapon Type's I through IV. Short Ranges (DNA, 
1972) 
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Figure 3-3. 
Neutron Dose FN as a Function ofSlant Range from a 1-kt Surface Burst, Weapon Types I through IV, Long Ranges (DNA, 
1972) 
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Figure 3-4. 
Neutron Dose FN as a Functwn ofSlant Range from a 1-kt Surface Burst. Weapon Types V through VIII, Short Ranges (DNA, 
1972) 
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Figure 3-5. 
Neutron Dose FN as a Function ofS/ant Range from a l-kt Surface Burst, Weapon Types V through VIIT. Long Ranges (DNA, 
1972) 
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Figure 3-6. 
Burst Height Adjustment Factors H j for Neutrons and Secondary Gamma Rays (DNA 1972) 
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the total dose absorbed-both prompt and delayed ef- 
fects are important because the effect is cumulative. 
b. Prompt gamma ray dose on equipment. Prompt 
gamma radiation that effects electronic equipment 
can be expressed by an equation of the form 


(3-4) 


where 
K 1, K2 
Constants 
W 
Yield in kt 
R 
Slant range in cm 
e 
Air density in gm/cm3 
Other sources of gamma rays include inelastic scat- 
tering of neutrons and the decay of fission products. 
These sources depend primarily on weapon design 
(SAl, 1971). Equation 3-4 will generally apply to 
materials with atomic number 20 or less. 


(1) The magnitude of the constants in equation 
3-4 should be determined from the upper bound 
estimates of the prompt gamma ray flux equation 
(4) 
). Moreover, the source energy density MeV/cm2 


~ 
d' 
should be converted to absorbed er.ergy expresse 
In 
rays (by the conversion factor presented in the 
reference), and the result multiplied by 2 shakes to 
obtain the prompt gamma ray dose. 
(2) The uncertainty associated with equation 3-4 
is expressed by 


where Q[ is the uncertainty in the form of equation 
3-4, Qw is the uncertainty in the yield, etc. Assume 
that Q[ :::: 0.2, Q e :::: 0.05, and that Q K 
= 0.3 (using 
equation 2-19 and the assumption eDNA, 1972) that 
the bounds are within ± 33% of the mean). 
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c. Secondary gamma ray dose on personnel. In ad- 
dition to the neutron dose described in paragraph 3-2 


C, secondary gamma rays and fission-product gamma 
rays will be significant to personnel. If the type of 
weapon is known and if the slant range is less than 
5000 yd (4572 m), it is possible to estimate the secon- 
dary gamma ray dose. For a low-altitude burst in a 
standard atmosphere where the personnel are located 
on or near the earth's surface, the secondary gamma 
ray dose is 


(3-6) 


where 


Weapon yield, kt 
Normalized 
secondary 
gamma 
ray 
dose, yd2 rad/kt ([m2(J/kg)/ktD 
Detonation height correction factor 
Slant range, yd (m) 


(1) The 
mean 
values 
of equation 
3-6 
are 
calculated by referring. to figures 3-6 through 3-10 
for the standard division of weapon types (DNA, 
1972), I through VIII. 
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Figure 3-7. 
Secondary Gamma Ray Dose Fp as a Function ofS/ant Range from a 1-kt Surface Burst, Weapon Types I through IV, Short • 


Ranges (DNA, 1972) 
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Figure 3-8. 
Secondary Gamma Ray Dose F,s as a Function ofSlant Range from a 1-kt Surface Burst, Weapon Types I through IV, 
Ranges (DNA. 1972) 
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Figure 3-9. 
Secondary Gamma Ray Dose Fys as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types V through VIII, 
Short Ranges (DNA, 1972) 
• 


3-10 


• 


1M 5-858-2 


SLANT RANGE, 
m 
2000 
2S00 
3000 
3S00 
4000 . 


107;:=:::I::::::;::::;:==;::::::=I::;===:::;:::==~I====;::::::==:::;::I==::::;::==::::;=I==:;===~~::; 


..... 
106 
104 


~ 
..... 
........ 
~ 
., 
"0 
........ 
rc 
l- 
V 
C'l 
,:,t, 


N 
........ 
"0 
..., 
>- 
- 


LiJ 
lOS 
10 3NE 


V'l 
0 
LiJ 
Cl 
Vl 
0 
>- 
Cl 
.e::t 
c::: 
;)- 


e::t 
c( 
c::: 
::E: 


102 , 
::E: 


104 
e::t 
• 


e::t 
::E: 
c.::l 
::E: 
e::t 
>- 
c.::l 
c::: 
e::t 
:>- 
Cl 
. c::: 
z 
e::t 
, 0 
Cl 
U 
Z 
LiJ 
0 


Vl 
10 1 
u 


103 
LiJ 
Cl 
Vl 
Z 
=> 
Cl 


0 
Z 
c::: 
=> 
c.::l 
0c::: 
c::: 
t.::l 


e::t 
c::: 


l/1 


102 
100 
e::t 
lJ...r- 


., 


101 . 
10- 1 


2000 
2S00 
3000 
3S00 
4000 
4S00 
SOOO 


R SLANT 
RANGE, yd 
• 


Figure 3-10. 
Secondary Gamma Ray Dose Frs as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types V through VIII. 
Long Ranges (DNA, 1972) 
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(2) The uncertainty associated with equation 3-6 
• 


except for megaton weapons where R < 1000 yd 
(914 m), then 


H2W2 
Q2 = 0.09 + Q2 + _I_ 
T 
W 
D2R4 
T 


[0.09 F~ + (0.09 + 4Q~)(FN + F ys) 2] 


0.36 W2F2 H2 H2 E2 
+ 
f 
yt~R 
W 
(3-12al 
D2 
T 


e. Total personnel dose. The total dose to personnel 
is the sum of the neutron dose, the secondary gamma 
ray dose, and the fission-product gamma ray dose. 
(1) The total mean dose to personnel is the sum 
of equations 3-1, 3-6, and 3-9. 
(2) The uncertainty Q? of the total dose to person- 
nel is the uncertainty of the sum of equations 3-1, 
3-6, and 3-9. Thusly, 


(3-8) 
Q2 
= 0.22 + Q2W + 4QR2 
Dys 


where Q[ is the uncertainty in the form of the equa- 
tion, Qw is the uncertainty of the yield, etc. Consider 
that Qf ::::: 0.2, QF 
= 0.3 and using equation 2-19 
with the assumpti~n that figures 3-7 through 3-10 
represent mean values of data uniformly distributed 
with bounds of ±25% (DNA, 1972), and QH1 ::::: 0.3. 
Thus, equation 3-7 reduces to 


d. Fission-product gamma ray dose on personnel. If 
the weapon yield, height of burst, slant range, and 
air density are known, the fission-product gamma ray 
dose absorbed by human tissl1e can be estimated. The 
fission-product gamma ray dose Dyf is 


(3-9) 


then equation 3-10 reduces to 


where Qf is the uncertainty of equation 3-9, QWf is 
the uncertainty of the yield, etc. Using equation 2-19 
and the assumption that figures 3-11 through 3-13 
represent mean values of data that are uniformly 
distributed and whose bounds are within ±25% of the 
mean values (DNA, 1972); and assuming other values 
of uncertainty, 


• 


• 


(3-12bl 


H2W2 
Q2 = 0.17 + Q2 + _I_ 
T 
W 
D2R4 


[0.09 F~ + (0.09 + 4Q~)(FN + FYSrl] 


0.36 W2F2 H2 H2 E2 
+ 
f 
yt-R 
W 
D2 


f 
Neutron-induced gamma activity in soils, residual 
effect on personnel. Neutron-induced soil activity must 
be considered for a scaled burst height below 75 
ftlktO,35 (23 mlktO.35) but above the ground. Soils that 
contain sodium, manganese, and aluminum will con- 
tribute most to the induced radioactivity. For times 
after burst of less than one-half hour, aluminum is 
most important; between one-half hour and five 
hours, manganese is most important; at later times, 
sodium becomes increasingly more important. 
(1) Figure 3-21 presents upper and lower bounds 
on the neutron-induced gamma dose rate at a 
reference time one hour after burst. At any given 
slant range, assume that the mean is the logarithmic 
average of the upper and lower bounds. For a time 
other than one hour after burst, this mean must be 
multiplied by the factors shown in figure 3-22 for the 
four types of soil presented in table 3-1. Additionally, 
all data must be multiplied by the factors shown in 
table 3-2 to obtain a fiqal estimate of the mean dose 
rate. 
(3-11) 


where 
Wf 
Fyf 
Fission yield (kt) 
Normalized 
fission-product 
gamma 
ray dose, radlkt ((Jlkglkt)] 
Slant-range-dependent 
adjustment 
factor 
Hw 
Yield-dependent adjustment factor 
E 
Hydrodynamic enhancement factor 
(1) The mean value of equation 3-9 can be 
calculated by referring to figures 3-11 through 3-20. 
In these figures, Q= Q/Qo' where Q and Qo (= 1.2250 
x 10-3 gm/cm3) are, respectively, the mean actual 
and ambient standard air densities (15.5 °C and 14.7 
psia). 
(2) The uncertainty of equation 3-9 is 


3-12 


• 


1M 5-858-2 


SLANT RANGE, 
m 


500 
1000 
1500 
2000 
2500 


104 
;::J:::::~==::::;::~I==:;===:;=~~::::;-===;==:::'::;==::::;::=~I~~ 102 


10-2 .l_---J._---"_----l..._---I...._--.L-_--L._....!-~~____3to..~~ 


500 
1000 
1500 
2000 
2500 
R SLANT RANGE, yd 


Figure 3-11. 
Fission Product Gamma Ray Dose Fyf as a Function ofSlant Range from d 1-kt (Fission Yield) Surface Burst, Short Ranges 
(DNA,1972) 
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Figure 3-12. 
Fission Product Gamma Ray Dose Frf as a Function ofSlant Range from a l-kt !Fission Yield) Surface Burst, Intermediate 
Ranges (DNA, 1972) 
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Figure 3-13. 
Fission Product Gamma Ray Dose Fylas a Function ofSlant Range from a 1-kt (Fission Yield) Surface Burst, Long Ranges 
(DNA, 1972) 
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Figure 3-14. 
Range-Dependent Burst Height Adjustment Factors HR for Fission Product Gamma Rays (DNA, 1972) 
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Figure 3-15, 
Yield·Dependent Burst Height Adjustment Factors H w for Fission Product Gamma Rays (DNA, 1972) 
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Figure 3-16. 
Fission-Product Gamma Ray Hydrodynamic Enhancement Factors E as a Function ofSlant Range for Relative Air DenSIty 
of 1,1 (DNA 1972) 
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Figure 3-17. 
Fission-Product Gamma Ray Hydrodynamic Enhancement Factors E as a Function ofSlant Range for Relative Air Density 
of 1.0 (DNA 1972) 
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Figure 3-18. 
Fission-Product Gamma Ray Hydrodynamic Enhancement Factors E as a Function ofSlant Range for Relative Air Density 
of 0.9 (DNA, 1972) 
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Figure 3-19. 
Fission-Product Gamma Ray Hydrodynamic Enhancement Factors E as a Function ofSlant Range for Relative A ir Density 


. of 0.8 (DNA, 1972) 
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Figure 3-20, 
Flssian,Praduct Gamma Ray Hydrodynamic Enhancement Factors E as a Function ofSlant Range for Relative Air Density 
of 0.7 (DNA,'1972J 
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Figure 3-21. 
Neutron-Induced Gamma Dose Rate Drl as a Function ofStant Range at a Reference Time of1 Hour After Bur$t (DNA, 1972) 
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Figure 3-22. 
Decay Factors for Neutron-Induced Gamma Activity, Drl (DNA, 1972) 
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Table 3-1. 
Soil Types by Composition (DNA, 1972) 


Percentage of Soil Type (by weight) 


Type IV 
Type I 
Type II 
Type III 
(beach 
sand, 
(Liberia, 
(Nevada 
(lava, clay, 
Pensacola, 
Element 
Africa) 
desert) 
Hawaii) 
Hodda) 


Sodium 
- 
1. 30 
0.16 
0.001 
Manganese 
.. 
0.008 
0.04 
2.94 
- 


Aluminum 
7.89 
6.90 
18.79 
0.006 


Iron 
3.75 
2.20 
10.64 
0.005 


Silicon 
33.10 
32.00 
10.23 
46.65 


Titanium 
0.39 
0.27 
1. 26 
0.004 


Caicium 
0.08 
2.40 
0.45 
- 


Potassium 
- 
2.70 
0.88 
- 


Hydrogen 
'0.39 
0.70 
0.94 
0.001 


Boron 
- 
- 
- 
0.001 


Nitrogen 
0.065 
- 
0.26 
- 


Sulfur 
0.07 
0.03 
0.26 
- 


Magnesium 
0.05 
0.60 
0.34 
- 


Chromium 
- 
- 
0.04 
- 


Phosphorous 
0.008 
0.04 
0.13 
- 


Carbon 
3.87 
- 
9.36 
- 


Oxygen 
50.33 
50.82 
43.32 
53.332 


• 


Table 3-2. 
Multiplying Factors (DNA, 1972) 


Multiplying 
Soil Type 
Factor 


I 
1.0 
II 
9.1 
III 
109.0 


IV 
0.024 


(2) In a particular situation, it may be necessary 
to construct a representative decay factor curve based 
on the chemical composition of the soil as compared 
to the data presented in table 3-1. For example, 
assume that a particular site has properties in which 
the percentage of aluminum most nearly matches 
that of Soil I but that manganese and sodium are most 
accurately represented by Soils II and III, respec- 
tively. Construct a decay factor curve from figure 
3-22 which approximates Soil I for times up to one- 


half hour, Soil II from ~.me-halfhour to five hours, and 
sodium after four hours. 
(3) Obtain the total gamma ray dose experienced 
by personnel from figures 3-23 through 3-26, which 
are applied to the one-hour dose rate presented in 
figure 3-21. The factor obtained from the appropriate 
soil type in figures 3-23 through 3-26 depends on the 
time of entry into a contaminated area and the 
amount of time spent in that area. These factors are 
multiplied by the appropriate dose rate in figure 3-21. 
(4) As an example, for the particular soil type 
discussed in (2) above, assume that the area is to be 
entered after two hours and to be occupied for seven 
hours. Enter figure 3-24 (for SoilIIJ for an entry time 
oftwo hours and an occupancy time of three hours 
(period during which manganese contribution is most 
important). Read the multiplying factor. Then enter 
figure 3-25 (for Soil Ill) for an entry time of five hours 
and occupancy time of two hours (period during which 
sodIum is most important) and read the multiplying 
factor. Add the factors and multiply this number by 
the"mean value of data in figure 3-21. 
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FIgure 3-23, 
Total Radiation Dose Receive'd (radJkt) in an area contaminated by Neutron-Induced Gamma Activity, Soil Type I (DNA, 1972) 
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Figure 3-24. 
Total Radiation Dose Received (radJktJ in an area contaminated by Neutron-Induced Gamma Activity, Soil Type II (DNA 1972) 
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Figure 3-25. 
Total Radiation Dose Received (radJkt) in an area contaminated by Neutron-Induced Gamma Activity, Soil Type III (DNA, 1972) 
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Figure 3-26. 
Total Rad!ation Dose Received Imd/ktl in an area contaminated by Neutron-Induced Gamma Actwity, Sol! Type IV (, 


loU 


101 
INFINITE TIME 


~«a:. 


loU 
II) 
0 
0 
a:. 
::l: 


I 


0I- 


~ 
• 


..J 
a.. 
100 
a..« 


~ 
a:. 
0 
~ 
u« 
LL. 


<.-' 
Z- 
~ 
..J 
a.. 
I- 


10- 1 
..J 
:;)x: 


PERIOD IN/ 
CONTAMINATED AREA 
- 
I 


• 


(5) As an estimate of uncertainty of the neutron- 
induced gamma activity 
in soils 
QD l' use the 
t · 
y 
.equa IOn 


TM 5-858-2 


Q2 
= 
Q2 + 0.00003 R2 Q2 
Dy! 
r 
R 


-3-30 


(3-13) 


wh~re Qr is the uncertainty of total dose and QR is 
the uncertainty in slant range (yd). Assuming that 
the dose rate is log-uniformly distributed in figure 
3-21, calculate Qr from equation 2-19 where L1 and 
L2 are the lower and ,upper bounds at any slant range 
in figure 3-21. 
• 


• 
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c. Hardened 
ground-based 
or 
shallow-buried 
systems are subject to EMP caused by high-altitude 
bursts or by contact' (or penetrating) bursts. Each 
system must be investigated to determine which burst 
condition predominates. Hardened ground-based or 
very shallow-buried systems will be susceptible to 
EMP caused by contact (or penetrating) bursts only 
if within the deposition region. The deposition region 
is that region in space where EMP fields 
are 
generated. 


4-2. Surface burst. 


a. In a surface burst, the conducting ground provides 
an effective return path for electrons traveling out- 
ward in the air from the burst point. These current 
loops produce very large azimuthal magnetic fields 
in the deposition region, especially near the ground. 
Large electric fields are also produced. The EMP for 
a surface burst will be one to two orders of magnitude 
greater than for air bursts. 
b. The magnitude and shape ·of the electric and 
magnetic fields in the deposition region cannot be 
described by simple functional relationships because 
of the complex nature of the underlying variables. For 
maximum accuracy, the currents induced in a system 
should be calculated with the recognition that the 
system, itself, influences the characteristics of the 


4-L ,'Gross phenomenology: 
. 


a. Time-varying electric and magnetic fields are pro- 
duced by nuclear detonations. Their frequency spec- 
trum varies from direct current to about 10 GHz. The 
gamma radiation produced by the nuclear detonation 
impinges on air molecules and induces a Compton 
current via'ionization ofthe air. As the Compton cur- 
rent expands away from the source, separation of the 
ions occurs and produces a radial electric field in the 
entire deposition region. Because of the atmospheric 
density gradient, the earth's magnetic field, and the 
weapon configuration itself, a magnetic field is also 
created, resulting in the radiation of an elec- 
tromagnetic wave for a medium-altitude (where the 
deposition region does not touch the earth) or a high- 
altitude burst. In a surface burst, the relatively large 
conductivity 
of the 
earth 
induces 
electric 
and 
magnetic fields (as shown in figure 4-1) whose inten- 
sities are much larger than for an air burst. 
b. The basic variables contributing to the EMP for 
exposed, flush, and shallow-buried facilities are 
- Weapon design 
- Weapon yield 
- Burst elevation 
- Atmospheric conductivity 
- Ground conductivity 


• 


• 


, 
' 
-,:. 


CHAPTER 4 


ELECTROMAGNETIC PULSE (EMP) 
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Figure 4-1. 
Surface-Burst EM? 
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EMP environment. Exact solutions for individual 
cases can be solved by using the computer programs 
described in the DNA EM? Handbook (GE·TEMPO, 
1971) and in the NWE Computer Code Directory (GE· 
TEMPO, 1974). For approximations, the relative 
changes in peak values of both the radial and polar 
electric fields and the azimuthal magnetic fields 
created in the deposition region are presented in 


figures 4-2 through 4-4 for two weapon yields and 
two ground conductivities. The time histories for 
these same fields are presented in figures 4-5 through 
4-7. Their corresponding frequency 
spectra are 
presented in figures 4-8 through 4-10. These data 
should be used for guidance only. Situations involv- 
ing other weapon yields or conductivities will require 
the execution of the appropriate computer program. • 
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Figure 4-2. 
Peak Magnetic,Field Bt Vs. Overpressure at the Air/Ground interface for Varying Ground Conductivities and Yields (DNA. 
1972) 
• 


. -4-2 


• 


• 


• 


LLJ 
....J 
<t 
u 
VI 
'"o 
....J 
.. 
E 


......... 
:> 


Cl 
....J 
LLJ 


LL. 


U 


0::: 
I- 
u 
LLJ 
....J 
LLJ 


....J 
<t 


Cl 
<t 
0::: 


::.::: 
<t 
LLJ 
a. 


0::: 


LLJ 


YIELD, 
Mt 
ag' mho/m 


W1 
10-2 
--- 
W1 


10-4 
--- 
W2 


10- 2 


P 
OVERPRESSURE, 
ps i 
LOG SCALE 
so 


Figure 4-3. 
Peak Radial Electric Field E R Vs. Ouerpressure at the 
Air/Ground Interface for Varying Ground Conduc- 
th'ities and Yields (DNA, 1972) 
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Figure 4-4. 
Peak Transverse Electric Field E e Vs. Overpressure 
at the Air/Ground Interface for Varying Ground Con- 
ductivities and Yields (DNA; 1972) 
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Figure 4-5. 
Eo· Time Waveform at the AIr/Ground Interface for Several Overpressure Levels r'SAJ. 1971) 
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Figure 4-6. 
ER-Tlme Waveform at the Air/Ground Interface for Several Overpressure Levels (SAl, 1971) 
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Figure 4-7. 
Ee-Time Waveform at the Air/Ground Interface for Several Ov~rpressure Levels (SAl, 1971) 
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Figure 4-8. 
Fourier Amplztude of B~ Wal'eform at the Air/Ground Interface for Sel'eral Ol'erpressure Lel'els (DNA, 1972) 
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Figure 4-9. 
Fourier Amplitude of E R Waveform at the Air/Ground Interface for Several Overpressure Levels (DNA, 1972) 
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Figure 4-10. 
Fourier Amplitude of E e Waueform al the Air/Ground Interface for Seueral Overpressure Leuels (DNA, 1972) 
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• 
c. Facilities buried in the earth will experience a 
decrease in the amplitude of the EMP fields that oc- 
cur at the air/earth interface. Facilities elevated 
several meters above the surface will experience a 
general (but not universal) attenuation of the EMP 
fields. 
Figures 
4-2 
through 
4-10 
illustrate 
characteristic curves. The uncertainty is derived as 
described in d below. 
d 
The uncertainties in predictions for surface-burst 
EMP depend on many variables; but for the threats 
of interest in this volume, the prime factors are over- 
pressure and weapon yield, frequency, and calcula- 
tional procedures. These uncertainties are sometimes 
referred to as tactical and calculational errors (AFWL 
EMP Handbook, 1972), which correspond to the terms 
"random" and "ignorance" uncertainties as used in 
this manual. Typical uncertainty bands for certain 
conditions are presented in figures 4-8 through 4-10. 
Typical examples are presented in figures 4-11 
through 4-13 for the uncertainty bands on the 
Fourier amplitudes of the B~, ER, and Eefields at the 
Ps02 range at a yield of WI and conductivity of 10-2 
mho/m. (Similar curves are available in the AFWL 
EMP Handbook.' for other overpressure regimes and 
yields.) These data may be used directly to assess the 
uncertainties of EMP environments. 


4-10 


e. The method used for determining uncertainty 
throughout this volume is embodied in equation 2-16. 
Because of the lack of useable information, it is 
customary to neglect the second term (relating cor- 
relation between uncertainties). Thus, 


N 
L 
n=l 


(2-16) 


where x = g (5'1,5'2' ... Yn ) has the same meaning as 
in equation 2-15. Wherever possible in this manual, 
the function g is defined in closed form so that the 
derivatives ~ ~ 
can be directly evaluated for each 
independent vaPiable (such as Pso' WI' Hz, Og)' How- 
ever, for EMP the form of the equations representing 
the data presented in figures 4-2 through 4-10 is too 
complicated for functional description. Accordingly, 
an alternate method is required, as described in f 
below. 
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Figure 4-12. 
Uncertainty Bands (0 = 1) of the E R Fourier Ampli- 
tude at the Air/Ground Interface (AFWL, 1972) 
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Figure 4-13.. Uncertainty Bands (0 = 1) of the Ee Fourier Ampli- 
tude at the Air/Ground Interface (A FWL, 1972) 


f 
Select a mean overpressure, mean weapon yield, 
mean ground conductivitY,and mean frequency 
(DNA, 1972; AFWL, 1972). For example, in'figure 4-8 
select Y1 = Pso = P~02' Y2 = W = WI' Y3 = Og = 10-2 


mho/m, and Y4 = f = fo' This defines point Ao on 
figure 4-8 in which x = I:B~ml = IB~lo' Vary Y1within 
the limits, say of PsOl and Ps03' which produces a cor- 
responding change in x(represented by the points B 
and C in figure 4-8). Thus, in equation 2-16, for the 
first independent variable y1 and for i1 
=:= 1, 


x2 = 
(IB~lo) 2 


(8g) 2::::: (dg) 2 = (IBhl-IB~31) 2 


(4-1) 


The user must select Qy 1 as the uncertainty he 
wishes 
to assign 
to the overpressure variable. 
Therefore, for the first component of uncertainty, 
equation 2-16 becomes 


Equation 4-2 is completed by performing the same 
operation for the variables W, ag, and f (if desired) 
and the summation in equation 4-2 is replaced as 
three individual terms corresponding to the uncer- 
tainties in these three variables. 
g. The uncertainty Qg in equation 4-2 is deter- 
mined from the data presented in figures 4-11 
through 4-14 (AFWL, 1972). Here, Q 
refers to the 
data scatter obtained from Monte Carro simulations 
performed by the AFWL. 
h. Note that the Monte Carlo mean values (denoted 
by the symbol e) differ from the nominal mean values 
(denoted by the solid lines). Using the nominal curves 
as a guide, construct new curves through the Monte 
Carlo means and bounds as exemplified in figure 4-11 
for 
IB~(f)1 by the dashed lines. These curves may then 
be used to estimate the uncertainties of the Fourier 
spectra at any frequency. 
i. From figure 4-11, Jor example, which cor- 
responds to figure 4-8, the uncertainty Qg is deter- 
mined from equation 2-18 where 


R = l(A'/A 
+ A IA')2 
4" 
0 
L 
u 
0 


• 


4 
-2(~ 
Q2 
+ L ynag 
2 Q2 
Pso 
n=2 _2 
8y 
Yn 
x 
n 
0 


2 


(4-2) 


n = 1 


Here, R is simply the squared averaged ratios of the 
bounds above and below the Monte Carlo mean. Note 
also that the mean value of IB~I (denoted by Ao in 
figure 4-11) was obtained from figure 4-8 and subse- 
quently used in equation 4-2, whereas Q 
was ob- 
tained froin the mean and scatter bands in tte Monte 
Carlo simulation. 
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Fourier'Amplitude ofE r Waveform ofHigh-Altitude 
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j. As noted in b above, the currents induced in a 
system should be cp.lculated by including the system 
within the EMP environment, i.e., by considering the 
'f"- 
' 
. 
,system and free"Jield environment to be coupled. 
When making EMP calculations using computer pro- 
grams, theuri~ertaintiesi~ the field calculations de- 
pend on the uncertainties in the mathematical model 
and the value of its parameters. With these.calcula- 
tions it is reccimme~ded that sufficient data be com- 
puted to allqw evaluation of equation 2':"15. This 
means that a fe~ auxiliary computations shou~d be 
made to allow for detfilrmination of the derivatives'in 
that equation as outlined in the previous paragraphs. 
The user will be required to ascertain for himself the 
nonrandom and random errors involved in such 
calculations. 


4-3. 
High-altitude burst. 


a. Any burst occurring at altitudes above 100,000 ft 
(30,480 m) is a high-altitude burst and may produce 
EMP fields that cover geographical areas as large as 
the continental United States.. ' 
. 


b. The electric field of the EMP is defined by 


, 
" 
~ 


(4-3) 


where Tl and T2 control the'rise and decay times of 
the pulse and Eo is the peak amplitude of the pulse 
(DNA EM? Handbook, GE-TEMP, 1971). It is noted 
that the amplitude of the EMP depends on (1) the 
prompt gamma ray yield, which is a function of 
nuclear weapon yield, and (2) the location of the 
system relat\ve to the detonation point. A typical 
pulse is shown in figure 4-15. The corresponding 
Fourier spectrum is presented in figure 4-14. The 
radiated magnetic field Br is proportional to Er (eq. 
4-3) according to the relation 
' 


B 
= E r . Wb/m2 ' 
(4-,-4) 


r 
Co 


where Co is the speed of light in a vacuum (3.0 x 108 
m/s). 
c. The data referenced in GE·TEMPO (1971) should 
be 
considered 
mean 
values 
of the 
free-field 
environmen't. 
. 
d. Accounting for the uncertainties of high-altitude- 
burst EMP predictions parallels the procedures 
discussed in paragraphs 4-2 d through i for determin· 
ing surface-burst EMPenvironments, and should be 
calculated accordingly. 
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Figure 4-15. 
E r-Time Waveform ofHigh-Altitude Radiated Signal 
(DNA, 1972) 
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CHAPTER 5 
FIREBALL 


log To = A + B log (R/WI /3) + C log2(R/W1I3) 


b. The mean radius is calculated from equation 5-1 
using CI = 250 for a contact burst or 220 for a sur- 
face intersecting fireball, 
nl = 0.35, and n2 = 
8.8 x 10- 6 . 
c. The uncertainty associated with equation 5-1 is 


(5-3) 


(5-2) 


D 


32.36 


C 


-235.4 


B 


566.4 


A 


-479.2 


where 
To 
Peak t~mperature in K 
R 
Range from GZ in ft 
W 
Yield in kt 
A,B,C,D are constants 
c. The mean peak fireball temperature in equation 
5-3 is given by the following coefficients for mean 
peak fireball temperature: 


where Q[ is the uncertainty in the form of equation 
5-1, Qc is the uncertainty ofthe coefficient CI , etc. 
'Assume\hat Q[::::: 0.2 and QCl = 0.15 (from equation 
2.16 in which n = 2 and R = 1.86). 


5-3. Temperature and density. 


a. The temperature and density within the fireball 
are shown in figures 5-1 and, 5-2 as a function of 
range 
for 
a 
family 
of early 
times. 
Fireball 
temperature is shown in figure 5-3 as a function of 
range for a series of later times, and in figure 5-4 as 
a function of time for a family of ranges. 
b. The peak fireball temperature vs. range can be 
represented by an equation of the form 


5-1. Gross phenomenology. 


a. In the explosion of nuclear weapons, a large 
amount 
of energy 
(kilotons 
to 
megatons 
TNT 
equivalent) is released in a time much less than a 
millionth of a second and in a small volume extend- 
ing slightly beyond the weapon. The resulting high- 
energy density produces a high-temperature (millions 
of degrees Kelvin, 
OK), high-pressure, thermally 
radiating plasma core. This core (the fireball) expands 
hydrodynamically because of the high internal 
pressure 
and 
radiates 
because 
of the 
high 
temperature. The radiation, primarily in the x-ray 
region of the electro-magnetic spectrum, is absorbed 
and reradiated by the surrounding air. Temperatures 
in the neighborhood of 1,000,000 OK persist for a few 
hundred microseconds. As the fireball expands and 
its temperature drops, radiation diffusion becomes 
less important. Hydrodynamic expansion, which pro- 
duces compression heating ofthe air behind the radia- 
tion front, overtakes the radiation front, developing 
a 
strong 
shock that 
continues to 
heat 
air to 
temperatures in excess of 10,000 OK for tens of 
milliseconds. 
b. The variables affecting the potential destruc- 
tiveness of the fireball at a ground surface location 
are the weapon yield and slant range. The larger the 
weapon yield, the greater the temperature and the 
engulfment duration will be at a given range. If the 
burst elevation is greater than the fireball radius, 
then engulfment will not occur and the effects ofther- 
mal radiation will control the design. How to deter- 
mine the effects from thermal radiation is presented 
in chapter 6. 


c. Immersion of the facility or its components in the 
fireball, or ingestion of fireball gases by the facility 
can lead to debilitating thermal loads. Knowledge of 
the fireball size and its time/temperature/density 
characteristics are required to determine the design 
strategy-avoidance or resistance-and the degree of 
protection required. 


,- 


• 


5-2. Engulfment radius. 


a. The maximum fireball radius is represented by an 
equation of the forrr-. 


where W is weapon yield in kilotons (kt); z is HOB 
in ft; and Cl , nl' and n2 are constants. Within this 
radius, facilities will be engulfed by the fireball. The 
temperatures within the fireball are presented below 
in paragraph 5-3. 
• 


R 
= C Wlll ll2Zft 
max 
1 
e 
(5-1) 


d. The uncertainty associated with equation 5-3 is 


Qi = QF+ 5.3 [A2 + B2 log2 (R/W1I3)+ C2log4 (R/W1I3) 


+ D2 log6 (R/W1I3~ Q~o + [B + 2C log (R/W1I3) 


+ 3D log2 (R/W1I3~ 2(Q2 + Q2 19) 
(5-4) 
J 
R 
W 


where Q[ is the uncertainty of the functional form of 
equation 5-3, Qc 
is the uncertainty of the coeffi- 
cients, QR is the u~certaintyof the range. Although 
5-1 
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Figure 5-1. 
Fireball Temperature Vs. Radius at Early Times Ul the Fireball ,History (adapted (ronz Brode, 1964) 


• 


.... 
• 


" 
.. 
• 


• 


;~ 
'. 
• 
• 


o 
101 I 
I 
10 
I 
I 
I 
I 
I 
I 
I 


SCALED RANGE, m/kt 1/ 3 


20 
30 
40 
I 
I 
I 
I 
I 
---l 
50 
I 
I 


7.4 


SCALED TIME, mS/kt1/ 3 


"o. 14 
0.0075 
100 
o 


~ 
a:: 


>- 
l- 


V) 
Z 
UJ 
o 


0.0 
........ 
-1 
0. 
10 
I 
I 
J 
I 
I 
I 
.<1 
I 
-I 


10-2\' 
I 
I 
I 
I 
I 
-I 
I 
I T 
I 
I, 1 
I 
I 
I 
I 
I J 
o ' 
'30 
60 
90 
120 
150 
180 


R/W 1/ 3 SCALED RANGE; ft/kt1/ 3 


QI 
, I 
,W 


Figure 5-2, 
Fireball Density Vs, Radius at Early Times in the Fireball History (adapted from Brode, 1964) 
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Latc Fircball Tcmperature Vs. Radius (adapted (rom Brode, 1964) 
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Figure 5-4. 
Temperature Vs. Time at High Peak Overpressures (adapted from Brode, 1964) 


insufficient data exist to accurately assess the uncer- 
tainty of equation 5-3 at this time, it is recommended 
that Qr ::::: 0.2 and Qeo ::::: 0.3. 
e. Duration of temperature is of interest. Since the 
fireball is like a bubble in the atmosphere, it begins 
to rise, and pulls away from the earth's surface in only 
a few tenths of a secondlktl/3 . Using a very approx- 
imate model for the effect of buoyant fireball rise on 
the temperature history at a scaled ground range of 
265 ftlkt1l3 (80.8 m1kt1l3), it appears that the late-time 
temperature environment will be reduced in approx- 
imately the manner indicated by the decreasing tail 


in figure 5-4. Thus a sharp temperature decrease 
begins after 0.4 or 0.5 slkt1l3 of exposure, and in 1.5 
to 2.0 slktl/3 the air temperature has returned to nor- 
mal ambient. Temperatures at other ranges would be 
similarly reduced at late times by the same effect. At 
the scaled ranges of 205 and 150 ftlkt 1l3 (62.5 and 
45.7 mlkt 113), temperatures increase even more but 
subsequently show a more rapid drop (at times less 
than 0.1 scaled second) because of thermal radiation 
loss, which becomes significant even before the 
fireball has begun to rise. 


• 


5-5 


• 


• 


• 


• 


TM 5-858-2 


CHAPTER 6 
THERMAL RADIATION 


a contact surface burst. The radius of the fireball at 
which the airblast shock breaks away and races 
ahead. The radiant exposure at closer points is 
negligible compared to the subsequent heat load.in 
the enveloping fireball gases, which is treated in 
chapter 5. 
(1) The mean radiant exposure is obtained from 
equation 6-1 where A = 86,000. 
(2) The uncertainty associated with the radiant 
exposure is 


Assume, Qg ::::: 0.24 for an air burst and::::: 0.3 for a 
contact burst. 
c. The transmittance T is the ratio of the radiant 
energy reaching a surface facing the fireball to the 
radiant exposure that the surface would receive if the 


f = C + Bh!Wl/3 
0 < h < 2.5 kftIMtl/3 
(6-3) 


where Qg is the uncertainty parameter associated 
with the functional form of Q, Qw is the uncertainty 
parameter associated with the weapon yield W, and 
so forth. Assume that Q2 + Q2 ::::: 0.1. The uncertain- 
ties of the transmittaJce, tt{ermal efficiency, and 
reflectance are developed in the paragraphs below. 
b. The thermal efficiency f of a weapon is the ratio 
of the thermal energy radiated by the fireball to the 
total yield of the weapon. The thermal efficiency of 
an ordinary nuclear weapon that is detonated in a 
low-altitude burst is 


Q2 = Q2 + Q2 + Q2 + Q2 + Q2 + Q2 + 4Q2 
(6-2) 
Q 
g 
A 
W 
f 
T 
F 
R 


where 
h 
= 
Height of burst, kft 
W 
= 
Weapon yield, kt 
Band C are constants 
(1) The mean thermal efficiency of an ordinary 
nuclear weapon that is detonated as a low-altitude 
burst is obtained from equation 6-3 where B = 0.078 
and C = 0.21. Assume that the mean thermal effi- 
ciency of an air burst is 0.4, and of a contact burst 
is 0.2. The thermal efficiency of a shallow-buried burst 
is small compared to that for a contact burst. A so- 
called enhanced weapon (enhanced nuclear radiation) 
exhibits lower thermal efficiencies than those cited 
above. 
(2) The uncertainty associated with the thermal 
efficiency of an ordinary nuclear weapon that is 
detonated as a low-altitude burst is 


Q2 = Q2 + h2B2 
(Q2 + Q2 /9\ 
(6-4) 
f 
g 
f2W2/3 \ 
h 
W J 


(6-1) 


T 


F 


Q = AWfTFIR2 cal/cm2 


where W 
f 
Weapon yield, Mt 
Thermal efficiency of the weapon (b 
below) 
Transmittance of the intervening at- 
mosphere (c below) 
Reflectance of clouds and the ground 
surface (e below) 
R 
Slant range, kft 
A 
Constant 
Equation 6-1 is applicable for slant ranges greater 
than about 1l00kT0.4 for an airburst or 145kT0.4 for 


6-1. 
Gross phenomenology. 


a. A thermal-radiation pulse is emitted by the ex- 
panding fireball. The thermal radiation received at 
a ground surface of design interest to this volume 
comprises energy in the ultraviolet, the visible, and 
the infrared spectral bands. The radiated thermal 
power reaches a peak in about 1 sIMtO.5. The pulse 
lasts, for all practical design purposes, about 10 
sIMtO.5. 
b. The basic variables affecting affecting the poten- 
tial destructiveness of the free-field thermal-radiation 
pulse are the: 
- Weapon design 
- Weapon yield 
- Burst elevation 
- Slant range 
- Transmittance by the intervening atmosphere 
- Reflectance of clouds and of ground surface 
The larger the weapon yield, the greater the height 
of burst, the smaller the slant range, the greater the 
transmittance, or the greater the reflectance, the 
greater will be the potential destructiveness from 
thermal radiation. 
c. The thermal-radiation pulse can cause fires and 
the degradation of the mechanical properties (possibly 
melting) of an exposed facility or facility component. 
The facility designers must either account for these 
effects in their design or preclude them by burying 
(or otherwise protecting) the facility or facility com· 
ponent. Knowledge of two characteristics of the free- 
field thermal-radiation pulse is required to implement 
design measures: Radiant exposure and radiation 
history. 


6-2. 
Radiant exposure. 


a. The total amount of radiant energy impinging on 
a unit area of a surface that is facing the fireball is 
called the radiant exposure ofthe surface. The radiant 
exposure is 


• 


• 
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intervening atmosphere were perfectly transparent 
and there were no reflecting surfaces. The transmit- 
tance for ground-surface targets is 


Q2 = Q2 + R2 exp C2bRN) [b(l + aRN) + aJ 2 
T 
f 
y2T2 


CQ~ + Q~) , 


h < 1.5 kft 
(6-6a) 


Q2 = Q2 + R2r 2 [(a + (3h + yh2)2 (Q2 + Q2) 
T 
f 
y2 
R 
Y 


+ h2«(3 +2yh)2 Q~ J ' 


1.5 < h ~ 10 kft 
(6-6b) 


• 


• 


(6-7) 


(6-8) 


y = 20 yS/4 
kft 
night 


where 


!J.y 


is 


Calculate the mean visibility V as shown in figure 
6-2. Use a base period of one year for calculating 
visibility. 
(2) The uncertainty associated with the visibility 


Coefficient of variation of the mean of 
y relative to its true value 
Coefficient of variation of y, as il- 
lustrated in figure 6-2. 
Assume that !J.y 
:::: 0.17. 
(3) The effect of prior bursts can be important. 
Bursts below about 2900 ftIMt°.4 create a dust cloud 
that will shroud the facility for many hours. Consider 
the effect of the dust cloud when calculating the 
visibility for each burst of a multiple-burst attack, ex- 
cept the first burst (chap. 13). 


e. The reflectance F is the ratio of the radiant 
energy reaching a surface facing the fireball when 
there are clouds and/or a reflecting ground surface, 
to the radiant exposure that the surface would receive 
in the absence of clouds and a reflecting ground sur- 
face. The reflectance is 


(l) Establish the facility-site, relative-frequency 
visibility function using visibility statistics obtained 
from a local weather station or from personnel at the 
site vicinity. For the seven classes of visibility 
delineated in table 6-1 the facility-site, relative- 
frequency function ~ is constructed as shown in figure 
6-2 such that the sum of ~O.8 + ~1.5 + ... b2 = 1.0. 
The subscripts refer to the seven classes of visibility. 
A typical relative-frequency factor is depicted in 
figure 6-2. Convert the weather station's night-time 
visibilities to equivalent day-time visibilities using 


(6-5b) 


(6-5a) 
T = (1 + aRN) exp (bRN) , h < 1.5 kft 


T = exp[Ca + (3h + yh2) rRN] , 


1.5 kft < h < 10 kft 


where 
R = Slant range, kft 
V = Daytime visibility, kft 
a, b, r, a, (3. and yare constants 
Equation 6-5 is based on the assumption that the at- 
mosphere is cloudless, but not necessarily clear of 
haze, fog, smoke, or dust. 
(1) The mean transmittance is obtained from 
equation 6-5 where a = 1.9, b = -2.9, r = -84.5, 
a = 3.0, (3 = -0.286, y = 0.01. Equation 6-5 is plotted 
in figure 6-1. The discontinuity in the transmittance 
shown in the figure is not physical, but is a conse- 
quence of the uncertainty in the two transmittance 
models for the two ranges of burst height in equation 
6-5. 
(2) The uncertainty associated with the transmit- 
tance for the two ranges of burst heights are 


• 


(6-9) 


(6-10) 


where 
Fca 
Reflectance of a layer of clouds above 
the burst 
Reflectance of a layer of clouds below 
the burst 
Fg 
Reflectance of the ground surface 
(1) The mean reflectance is obtained from equa- 
tion 6-9 where the mean values of Fca' Feb' and Fg 
are obtained from f, g, h below 
(2) The uncertainty associated with the reflec- 
tance is 


where Qf is the uncertainty parameter for the func- 
tional form of equation 6-5, QR is the uncertainty 
parameter for the slant range, etc. Assume Qf:::: 0.1. 
The 
uncertainty 
of the 
visibility 
is 
calculated 
subsequently. 
d. The U.S. Weather Bureau defines day-time 
visibility as the distance at which a prominent dark 
object, located near the ground and silhouetted 
against the sky, is visible and recognizable. Night- 
time visibility is defined as the distance an unfocused, 
moderately intense light source is visible. The day- 
time visibility for generic cloudless atmospheric con- 
ditions is given in table 6-1. Assume that the day- 
time and night-time visibilities are identical for iden- 
tical atmospheric conditions. 
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Figure 6-1. Mean Transmittance for a Ground-Surface Target (SA/, 1971) (l of 5) 
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Figure 6-1, Mean Transmittance for a Ground-Surface Target (SAl. 1971) (2 of 5) 
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Figure 6-1, Mean Transmittance for a Ground-Surface Target (SAl, 1971) (3 of 5) 
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Figure 6-1. Mean Transmittance for a Ground-Surface Target (SAl, 1971) (4 of 5) 
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Table 6-1. Day-time Visibility for Generic Cloudlessness (DNA,. 1972) 
• 


Type 


Very clear 


Clear 


Light haze 


Medium haze' 


Description 


This condition rare except at 
high-altitude locations. 


Sky deep blue. 
Shadows distinct, 
dark. 


Sky white; dazzling near sun. 
.Shadows visible, gray. 


Sky bright grayish-white. 
View 
sun without serious discomfort. 
Shadows visible but faint. 


Vi$ibilit~ 


mi 
(kml 


32.0 
5"J .5" 


16.0 
25.7 


8.0 
l.L9 


4.0 
6.4 
• 


Heavy haze 


Thin fog 


Light fog 


6-8· 


Sky dull gray-white. 
Sun's disc 
just visible. 
Shadows barely 
discernible. 


Sky light gray with maximum 
luminance around sun. 
Sun's disc 
not visible, no shadows. 


Sky dull gray with maximum 
luminance at zenith. 


3.0 
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0.8 


4.~ 


2.4 
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Visibility: Calculation of Mean and Coefficient of Variation (Adjust ordinate so that I~n = 1.0) 
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f 
When the burst occurs below a cloud layer, the 
reflectance is 


(6-11) 


where A and B are constants and ~ca is the fraction 
of a year that clouds appear in the site vicinity at a 
height above burst elevations of interest. 
(1) The mean clouds-above reflectance is obtained 
from equation 6-11 where A = 1.0 and B = 0.5. 
(2) The uncertainty associated with the clouds- 
above reflectance is 


(6-12) 


where 
A~ca is the coefficient of variation of. Fca 
relative to its true value. Assume that AFca == 0.1. 
g. When a layer of haze or clouds is below the burst, 
but high enough not to affect visibility, the reflectance 
Feb is as shown in table 6-2. For this condition, use 
V = 16 mi in the calculation of transmittance. 
(1) Calculate the mean value of the clouds-below 
reflectance as illustrated in figure 6-3. Use a base 
period of one year for calculating the mean value of 
Feb' The method of calculating the mean reflectance 
parallels the procedure described in d above. 
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", 


Table 6-2. Reflectance Feb for Layer of Clouds or Haze Below the Burst (DNA, 1972) 


Type 


Very clear 


Clear 


Light haze 


Medium haze 


Heavy haze 


Light cloud 


Medium cloud 


Heavy cloud 


Description 


This condition rare except at high- 
altitude locations. 


Sky deep blue. 
Shadows distinct, 
dark. 


Sky white; dazzling near sun. 
Shadows visible, gray. 


Sky bright grayish-white. 
View 
sun without serious discomfort. 
Shadows visible but faint. 


Sky dull gray-white. 
Sun's disc 
just visible. 
Shadows barely 
discernible. 


Sky light gray with maximum 
luminance around sun. 
Sun's disc 
not visible, no shadows. 


Sky dull gray with maximum 
luminance at zenith. 


Sky dark gray; brightness pattern 
gives no indication of sun's 
position. 


Re flectance 
- 
1.0 


0.95 


0.75 
• 


0.5 


0;4 


0.3 


0.2 


0.13 
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Dense cloud 
The low luminance level suggests 
the approach of nightfall. 
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'\ 
(2) The uncertainty associated with the clouds- 
below reflectance is 


...., 


• 


• 


(6-17) 


(6-15) 


Assume Qg = 0.1. 
b. The normalizing time to is 


to = AWE[l + Bh + Ch2Jo sec 


where Qw is the uncertainty parameter for the yield 
and Qh is the uncertainty parameter for the HOB. 
Assume Qg ;:: 0.17. 


where 
Q 
Total radiant exposure, cal/cm2 


t 
Time, seconds 
to 
Normalizing time, seconds (b below) 


a, (3 are constants 
(1) The 
mean fractional 
radiant-exposure 
is 
obtained from' equation 6-16 where a = -0.5 and 


(3 = 0.33. 
(2) The uncertainty associated with the fractional 
radiant-exposure is 


Q~t = Q~ + Q~ + (1 - ~t [at/to ~ 2(3(t/t)2] 2 Q~o 


(6-16) 


where 
W 
Weapon yield in Mt 
h 
Height ·of burst, kft 
A, B, C, E, a are constants 
(1) The mean normalizing time is obtained from 
equation 6-18 where A = 0.84, E = OA31,B = -0.0162 
C = 7.2 x 105 and 0= 0.764. The mean normalizing 
time can also be obtained from figure 6-5. 
(2) The uncertainty associated with the normaliz- 
ing time is 


Q2 = Q2 + (2 Q2 + o2h2[ 
B + 2Ch 
]2 Q2 
to 
g 
W 
1 + Bh + Ch2 
h 
(6-19) 


fitted to the fractional radiant-exposure history of 
figure 6-4 is 


Qt = Q [1 - exp [at/to + (3 (t/to)2JJ 
callcm2, 1 < t/to < 10 


(6-1) 


(6-13) 


(6-14) 


where 


!J.Fcb 
Coefficient of variation of the mean 
value of Fcb relative to its true value 
Coefficient of variation of Fcb' as 
calculated in figure 6-3 
Assume that!J.F 
;:: OAO - 0.35 Fcb' 
h. The groundCfurface reflectance in the facility-site 
vicinity is given by 


where!J.F is the coefficient of variation of Fg relative 
to its tru~ value. Assume that!J.F 
;:: 0.1. 
(3) Dust created by prior burstsgwill settle out and 
decrease the reflectance of a bright surface. Consider 
this effect.when calculating the reflectance for each 
burst of a multiple-burst attack, except the first burst 
(chap. 13). 


6-3. Radiation history. 


a. The history of the radiation energy is relatively 
complex, exhibiting several peaks and valleys. The 
details of the history are strongly influenced by the 
weapon yield, and the burst elevation. For purposes 
of design, the simple fractional radiant-exposure 
history shown in figure 6.-4 can be used in place of 
the radiation history. The functional form'of a curve 


where ~g is the fraction of the year that the ground 
surface appears bright., Snow and light-colored soil 
qualify as bright surfaces. 
(1) The 
mean 
ground-surface 
reflectance 
is 
obtained from equation 6-14 where K 1 = 1 and 
K 2 = 0.5. 
(2) The 
uncertainty 
associated with 
ground- 
surface reflectance is 


• 


. 6-12 
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Figure 6-4. Fractional Radiant Exposure 
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CHAPTER 7 
AIRBLAST 


• 


7-1. Gross phenomenology. 


a. Hydrodynamic expansion of bomb debris and 
fireball produces an expanding shock wave in the sur- 
rounding air. Interaction with the ground surface af- 
fects the detailed structure and geometry of the ex- 
panding shock wave. 
b. A contact burst at a scaled elevation no greater 
than 5 ftlkt1/3 (1.5 rnlkt1l3) produces a hemispherical 
shock front expanding supersonically. At scaled 
heights of burst between 5 and 160 ftlkt 1l3 (48.8 
rnlkt1l3), the incident spherical shock reflects from the 
ground (fig. 7-1); the reflected signal, traveling at a 
higher speed in air heated by the incident shock, even- 
tually overtakes the incident shock, forming a Mach 
stem that connects the ground with the triple point 
(which is the intersection o( incident, reflected, and 
Mach stem shocks). The dynamic movement of the air 
behind the Mach stem causes a second shock to form, 
producing double shock phenomena (Brode-Lewis, 
1975). Thermal-radiant heating of the ground ahead 


of the Mach stem causes the formation of a precursor 
shock, which trails back and merges with the Mach 
stem. The aspects of airblast important to the 
designer are the overpressure and dynamic pressure 
pulses. These are often described by the principal 
parameters: peak pressure, positive phase duration, 
and positive phase impulse (the time integral of the 
pressure history), Overpressure and dynamic pressure 
histories are also required; they are summarized in 
paragraphs 7-2 and 7-3. 
c. The basic variables governing airblast are: 
Weapon design 
- Weapon yield 
- Range from detonation 
- Burst elevation 
- Surface conditions 
- Surface contour 
- Structural configuration 
The larger the weapon' yield, the greater the height 
of burst up to the limit of the Mach, stem, or the 


.' 


• 


INCIDENT 
SHOCK ____ 


MACH REFLECTION 


Figure 7-1. Composition of Shocks Due to Nuclear Airburst (Brode. 1964) 


7-1 


Table 7-1. 
Thermally Near-ideal and Thermally Nonideal Surfaces (DNA, 1972) 


where R is the range from GZ in ft and W is yield in 
kt. The coefficients C 
and exponents np depend on 
scaled range, scaled HOB and the terrain characteris- 
tics. 


(1) Ground surfaces are characterized as therm- 
ally "near-ideal" (in which a precursor shock is not 
likely to form) or thermally "nonideal" (in which a 
precursor shock is likely to form). Table 7-1 presents 
a compilation of ideal and nonideal surfaces. These 
categories are used to define the parameters of equa- 
tion 7-1. 


(2) The mean values to be used in equation 7-1 
are presented in table 7-2 for ideal and nonideal sur- 
faces. The corresponding data from which the table 
was prepared are presented in figures 7-2 and 7-3. 
These data are also presented in figures 7-4 through 
7-10, which should be used directly for selecting 
mean values. 
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shorter the range, the greater will be the potential 
destructiveness of the airblast. The influence of the 
other variables is more complex; changes in weapon 
design and terrain conditions may amplify or at- 
tenuate the local airblast environment. 


d. At the ground surface or in exposed configura- 
tions, structures can be dealt a crippling or fatal blow 
by the airblast. The facility designer must either ac- 
count for the effects in the design or prevent the struc- 
tures' direct exposure to airblast. To determine the 
degree of resistance necessary and to support selec- 
tion of the design approach, knowledge of the airblast 
waveform, peak pressure, positive phase duration, 
and positive phase impulse are required. 


7-2. Free-field peak overpressure. 


a. Surface to moderately high airburst. The peak over- 
pressure PSG depends on scaled HOB (height of burst), 
scaled range from GZ (ground zero), and terrain 
characteristics. In general, the peak overpressure PSG 
may be defined as a piecewise function of scaled range 
for a given scaled HOB, as indicated in the following 
equation: 


Thermally Near-Ideal 
(precursor unlikely) 


Water 


Ground covered by white smoke 
layer 
Heat-reflecting concrete 


Frozen tundra 


Ice 


Packed snow 
Soil with sparse vegetation 


Commercial and industrial areas 


p = 1,2, ... , N 


Thermally Nonideal 
(precursor may occur) 


Desert sand 


Coral 


Asphalt 


Surfaces with low, thick vegetatio 
j 


Surfaces covered by a dark smoke 
layer 


Dark colored rock 


Most agricultural areas; 
residential areas in cities 


Dry soil with sparse vegetation 


,(7-1) • 


• 


• 


• 


• 


., 


• 


Table 7-2. Parameters for Peak Overpressure 


Scaled HOB 
Seal ed Range 


. ft/kt l /3 (m/kt 1/3) 
ft/kt l /3 
(m/ kt I /3) 
c 
n 
p 
p 
Near-Ideal Surface 


a 
(0) 
170 to 
330 
(52 to 
1011 
1.80 x 109 
-2.93 


330 to 1000 (101 to 
305) 
1.84 x 107 
-2.09 


1000 to 3800 (305 to 1158) 
1.55 x 106 
-1.73 


200 
(61) 
270 to 
330 
(82 to 
101) 
4.76 x 109 
-3.04 


330 to 1220 (101 
to 
372) 
6.85 x 106 
-1.92 


1220 to 4.400 
(372 to 1341) 
2.60 x 104 
-I. 67 


400 
(122) 
290 to 
500 
(88 to 
152) 
1.63 x 109 
-2.92 


500 to 1350 (152 to 
411 ) 
3.19 x 107 
-1.78 


1350 to 4900 (411 to 1494) 
1.04 x 106 
. -1 .63 


600 
(183) 
200 to 
520 
(61 to 
158) 
7,86 x 102 
-0,52 


520 to 1520 (158 to· 
463) 
· 
· 


1520 to 5600 (463 to 1707) 
1.06 x 106 
-I .61 


800 
(244 ) 
450 to 1500 (137 to 
457} 
6.92 x 103 
-0 58 


1500 to 6000 (457 to 1829) 
1.50 x 106 
-1.63 


1000 
(305 ) 
300 to 
800 
(91 to 
244) 
I. 65 x 10 2 
-0.42 


800 to 2000 
(244 to 
GIO) 
· 
· 
2000 to 6500 
(610 to 1981) 
5.79 x lOS 
-I. 5I 


2000 
(610) 
600 to 2500 (183 to 
762) 
8.09 x 101 
- 0.47 


2500 to 5300 (762 to 1615) 
3.17 x. 103 
-0.94 


Nonideal Surface 


0 
(0) 
ISO to 
340 
(46 to 
104 ) 
I. 18 x 109 
-2.79 


340 to 1000 (104 to 
305) 
I 95 x 107 
-2.09 


1000 to 3700 (305 to 1128) 
I. 91 x 106 
-1.76 


200 
(61 ) 
240 to 
750 
(73 to' 
229) 
4.89 x 108 
-2.68 


750 to 1200 (229 to 
366) 
2.40 x 102 
-0.48 


1200 to 4400 (366 to 1341) 
6.75 x lOS 
-I .60 


400 
(122 ) 
210 to 
480 
(64 to 
146) 
2.46 x lOS 
-I .46 


480 to 1650 
(146 to 
503) 
· 
· 
1650 to 5000 (503 to 1524) 
7.83 x 105 
-1.59 


600 
(183) 
140 to 
360 
(43 to 
110) 
7.21 x 102 
- 0.54 


360 to 
880 (110 to 
268) 
2.79 x 103 
-0.77 


880 to 5600 (268 to 1707) 
9.93 x lOS 
-I. 60 


800 
(244 ) 
300 to 
600 
(91 to 
183) 
2.19 x 102 
-0. ,12 


600 to 1500 
(1·83 to 
457) 
1,14 x 103 
-0.68 


1500 to 6000 (457 to 1829) 
4.02 x 105 
-1.48 


') 


1000 
(305) 
300 to 
800 
(91 to 
244) 
1.55 x 10· 
-0.41 


800 to 1920 (244 to 
585) 
4.16 x 102 
-0.56 


1920 to 6400 (585 to 1951) 
4.68 x 105' 
-1.49 


2000 
(610) 
550 to 2500 (168 to 
762) 
7.29 x 101 
-0.46 


2500 to 5300 (762 to 1615) 
I. 96 x 103 
-0.88 


'Note: 
Compute values from Figures 7-2 or 7-3. 
U.S. 
Army Corps of Engineers 


TM 5,:",858-2 


Figure 7-2. Peak Overpressure at Thermally Near-Ideal Surfaces 
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Figure 7-3. Pea'k Overpressure at Thermally Nonideal Surfaces 
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Figure 7-4. Peak Overpressures at a Near-Ideal Surface, Very High Overpressure Region (DNA, 1972; Brode-Lewis, 1975) • 
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Figure 7-5. Peak Overpressures at a Near-Ideal Siirfac.e, High Overpressure Region (DNA, 1972) 
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Figure 
7-6~ Peak Overpressures at a Near-Ideal Surface, Low Overpressure Region (DNA, 1972) 
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Figure 7-7, Peak Overpre"surl's at a Near-Ie/ea{ or Thermally No"ie/ea{ Sur!cu',', Very Low Overpressure R"t{io" (DNA, 1972) 
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Figure 7-8. Peak Overpressures at a Thermally Nonideal Surface. High Overpressure Region (DNA, 1972) 
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Figure 7-9~ Peak Overpressures at a Thermally Nonideal Surface, Low Overpressure Region (DNA. 1972) 


SCALED RANGE, m/kt 1/ 3 


0 
500 
1000 
1500 
2000 
l__ ~__ ~ 
5000,--, - ....-::2-'::'t.' 
d 
I 
I 
I 
I 
-,- 
I 
I 
I 
I 
If-I 500 


. 
...... 
.......'-'" 
~ 4000~ 
1-~I=~_P~l"'J -1--- r--i--l 
S~ 


I- 
" 
1000 
E~ 


VI 
3000 
l- 


e::: 
----- 
---- 
---- 
VI 
=> 
, 
e::: 
m 
, 
=> 
m 
I.L. 
, 


0 
, 
I.L. 
0 
l- 
I- 
:::c 
" 
<.:) 
2000 
~~- 
:::c 
--._- 
----~ 
<.:) 
- 
..... 
--- 
I.LI 
..... 4 
REGULAR 
---.., 
- 


:::c 
I.LI 
..... ...... REFLECTI ON 
500 
:::c 
0 
I.LI 
" 
- - 
0 
...J 
" 
/ 


I.LI 
c{ 
" 
-- 


...J 
... 
"" 
/ 
c{ 
u 


~~ 
1000 
7 
MACH 
./ 
u 


~ 
~. 
-- 
VI 
M 
REFLECTION 
./ 
........ 
I 
. 
- 
./ 
3: 
-- 


........ 
- 
.!: 


O~ 
I 
....... -...rL 
L 
I 
/ 
I 
............ 
I 
I 
...-! 
I 
I 
---.L- 
I 
I 
, 
I 
'LO 
0 
1000 
.2000 
3000 
4000 
5000 
6000 
7000 
8000 


R/Wl/~SCALED RANGE, ft/kT 1/ 3 
--- NONNUCLEAR DATA 
.~ PRECURSOR TRANSFORMING TO 
IDEAL SHOCK. 


"I-- 


• 


• 


I 
I 


10- 1 


I 
j 
I 
! 


"""1 
-< 


I 
..J 


I 
..J 
! 
I 


~ 
i 
I 
I 
i 
-J 


I 


I I I 
I 
I I I! 
1____1 
I 
I 
I 
I 


104 
105 


m/k t 1/3 


103 


I 
I II 


102 


R/W l13,SCALED RANGE, 
ft/kt l/3 


______'__..l_-'----.l-.J~.L.L'____ 
_ 
___l _ 
__'______'___'__.L...L...L.::Clo:..__ 
__._.l----1----'----l.--..L.LL.L.-'---..l 


10 1 


UJ>a: 
::::> 
u 


zo 


Cl 
UJ 
:::.:::a: 
c:x: 
~ 


UJa: 
::::> 
V1 
V1 
UJa: 
Q..a: 
w>o 


:::.::: 
<t. 
UJ 
Q.. 


o 
VI 


Q.. 


TM 5-858-2 


Figure 7-10. Peak Overpressure at the Surface from a Contact Surface Burst (DNA, 1972) 


(3) The uncertainty associated with the peak 
overpressure is: 


n 2 
Q2 
= Q2 + Q2 
+ n2Q2 + ---.E Q2 
Pso 
[ 
Cp 
p 
R 
9 
w 
(7-2) 


points other than at surface zero if the target remains 
in the region ofregular reflection, that is, if the burst 
point is no more than 30 deg (0.52 rad) offthe vertical. 
(1) The equation of the data presented in figure 
7-11 (construed to be the mean value) can be 
represented by: 


where W is the actual yield in kt, h is the HOB in 
kft, and a, 
~. and yare constants. Equation 7-4 is 


where Pso is the peak overpressure in psi, Rs is slant 
range in ft, We is the effective yield in kt, and A, B, 
C, and Dare constants. The effective yield is given 
by the Blast Efficiency Factor; 


where Q[ is the uncertainty associated with the func- 
tional form of Pso' Qc 
is the uncertainty of the coef- 
ficient, and so forth. 
P For practical use, the major 
uncertainty in equation 7-1 as a representation of ex- 
perimental and theoretical data is in the coefficient 
Cpo It is recommended that Q[:::: 0.05 be used and that 
Qcp = 0.2 (which assumes that the data are log- 
normally distributed and that the bounds of ±2 stan- 
dard deviations are located between 3/2 and 2/3 times 
the mean value). 
b. High-altitude bursts. Peak overpressures from 


high-altitud~'bursts (large HOB) are obtained from 
figure 7-11, which shows peak pressure as a function 
of scaled slant range. These data may be used at 


7-12 


W 
_e = a + ~h + yh2 
W 
(7-4) • 
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plotted in figure 7-12. It is preferable to use data 
directly from figure 7-11 and 7-~2 rather than from 
equations 7-3 and 7-4 to calculate the mean over- 
pressure and the mean effective yield. 
(2) The uncertainty of the overpressure at the 
ground surface due to a high-altitude burst as 
represented by equation 7-3 is: 


Q2 
=: Q2 + 5 3[A2 + B21og2(R /W1I3) 
Pso 
[ 
. 
5 
e 


+ C2Iog4(R /W1I3) + D21og6(R /W1I3) ] Q2 
5 
e 
5 
e 
Co 


+ [B + 2C log(R2/W~/3) 


Q2 
+ 3D log2(Rs/W~/3)]2 (Q~ + 
We) 
(7-5) 
. 
5 
9 


where Q[ is the uncertainty in the form of equation 
7-3, QR is the uncertainty of the slant range, etc. 
Assume\hat Q[;::: 0.1 to account for inaccuracies in 
fitting the function in equation 7..:.3 to the curve in 
figure 7-11. Also let 


Q2 
=: Q2 + 
1 
{[a2 + (32h2 + y2h4]Q2 
W 
W [P 
COF 
e 
a + (3h + Yh2J 


+ h2[(3 + 2yh]2 Q~} 
(7-6) 


The values ofthe constants in equations 7-3 through 
7-6 are A 
=: 12.0, B 
=: -6.08, C 
=: 0.976, and 
D =: -0.0608; for h ~ 40 kft, a =: 1.0, [3 = 0.0015, and 
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Figure 7-12. Blast Efficiency Factor for High-Altitude Bursts (DNA, 1972) 
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y = 0; for h > 40 kft, a = 0.956, [3 = 0.00101, and 
y = -3.55 
X 10-5. The uncertainty of QCOF is ob- 
tained by applying equation 2-19 to the data in figure 
7-12 where L2 and LI are the upper and lower 
bounds ofthe Blast Efficiency Factor at any particular 
HOB. Consider that QC 
= 0.13 to account for the 
uncertainty in the 
dat~ comprising Qgure 7-11 
(derived from equation 2-18 in which V R = n = 2). 


7-3. Free-field peak dynamic pressure. 


a. At GZ the ground-surface dynamic pressure is zero 
for all HOB. Outward from GZ, the peak dynamic 
pressure rises and then falls with increasing ra.'nge. 
Unlike 
overpressure, 
the 
dynamic 
pressure 
is 
relatively insensitive to the surface properties and 
depends more on the dust loading. Peak dyna~ic. 


TM5..,85.8":,~; 


pressure, a complicated function of HOB and range, 
is not easily represented by simple linear relation- 
ships, regardless ofthe manner in which the data are 
presented. 
Figures 
7-13 
through 
7-15 
present 
dynamic pressure curves in their simplest forms. In 
general, these curves may be represented by 


log Q = A(~)2 + B'(~) + C. 
I wI/3 
I w1I3 
I' 


Rj _ I ~ R < Rj 
i = 1, 2, ... , N 
(7-7) 


. 
' 
," 


where Q is the dynamic pressure in psi, R is the range 
from GZ in ft, W is yield in kt, and A, B, and Care 
coefficients to be defined in c below for various ROBs 
within certain scaled ranges. 
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Figure 7-13. Dynamic Pressure at a Perfect Surface (DNA, 1972) 
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Figure 7-14. Dynamic Pressure at a Light Dust Surface (DNA, 1972) 
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Figure 7-15. Dynamic Pressure at a Heavy Dust Surface (DNA, 1972) 
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b. It is recommended that the mean values for 
dynamic pressure be taken directly from figures 7-13 
through 7-15 or from the source data presented in 
figures 7-16 through 7-18. 


c. The uncertainty associated with equation 7-7 is 


Q~ = Q;+ 5.3HA2(RiW1I3)4 + B2(RiW1I3)2 + C2] Q~o 


Q2 
+ (RiW1I3)2 [2A(RiW1I3) + B]2 [Q~ + ;n 


(7-8) 


where Qr is the uncertainty in the form of equation 
7-7, Qc is the uncertainty of the data in figures 7-13 
o 


through-7-15, Qw is the uncertainty of the yield. Due 
to the complexity of the curves in figures 7-13 
through 7-15, the coefficients A, B, and C have not 
been calculated in advance. They can be determined 
for any range of the data by fitting equation 7-7 to 
the data within the range of interest. For example, 
if the curve corresponding to the scaled HOB of 400 
ft/kt1l3 in figure 7-14 is fitted by equation 7-7 in the 
scaled range 570 ~ RiW1I3 ~ 1200 ft/kt1l3, the approx- 
imate regression equation will be specified when 
A = 1.021 X 10-6, B = -3.485 X 10-3, and C = 4.013. 
Credible estimates of the values of Qc are not cur- 
rently available from the library of d:ta. However, 
for design purposes assume Qr = 0.1 and Qc 
:::: 0.3. 
o 
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Figure 7-16. Peak Dynamic Pressure at an Ideal Surface (DNA, 1972) 
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Figure 7-17. Peak Dynamic Pressure at a Surface with Light Dust Conditions (DNA, 1972) 
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Figure 7-18. Peak Dynamic Pressure at a Surface with Heavy Dust Conditions (DNA, 1972) 
(Dotted sections are taken from figure 7-17, for Light-Dust Conditions.) 
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7-4. Overpressure positive phase duration. 


a. The positive phase duration, that is, the time that 
the overpressure remains positive, is relatively in- 
sensitive to HOB or range. Reasonably accurate 
mathematical representations of the positive phase 
durations presented in figures 7-19 and 7-20 may be 
achieved by interchanging the variables on the 
ordinate with the variables on the family of curves 
and replotting the data as shown in figures 7-21 and 
7-22. Equations of the latter forms of the positive 
phase durations are 
• 


i = 1,2, , .. , N 


D+ 
_P_=C.+C 
~ 
R 
<R<R 
o 
D' 
i~l 
i 
W1l3 
I 
I W1I3 


(7-9) 
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where Qr is the uncertainty of the functional form of 
equations 7-9 and 7,...10, Qc is the uncertainty of the 
coefficients, and QR is the uncertainty of the scaled 
range. The value of the constants in equations 7-9 
through 7-11 are presented in table 7-3. Assume 
Qr ::::: 0.2 and Qc = 0.1 for Dt/Wl/3 < 0.25 slkt1l3 or 
Qc = 
0.17 for Dt/W1I3 > 0.25 slkt1l3 (which is 
predicated on the use of equation 2-17 where n = 2 
and R = 1.4 or R = 2, respectively). 


7-5. Dynamic 
pressure 
positive 
phase 
duration. 


a. Dynamic pressure positive phase dUrations are not 
as readily defined as the over-pressure positive phase 
durations. The usual practice is to modify the positive 
phase of the overpressure by certain factors. Thus, the 
dynamic pressure positive phase is approximated by 


for other selected scaled HOBs. In equations 7-9 and 
7-10, D; is the positive overpressure duration in sec, 
W is yield in kt, R is range in ft, and Co and CD are 
constants. 
b. It is recommended that the mean values of the 
positive phase durations be extracted from figures 
7-19 and 7-20 by interpolating the plotted data 
presented in table 7-3. 
c. The uncertainty of the positive phase duration 
for equation 7-9 is to be computed from: 


from figure 7-23, which is recast from figures 7-24 
and 7-25 where Co and CD are constants. The 
dynamic positive phase duration DJ"has units of sec, 
W is the yield in kt, and R is the range in ft. 
b. The mean values of the positive phase should be 
extracted from figures 7-24 and 7-25 by interpolating 
the plotted data. 
c. The uncertainty ofthe dynamic pressure positive 
phase duration is computed from 
• 


i = 1,2, ... , N 
(7-10) 
1 = 1,2, ... , N 
(7-13) 


where Qr is the uncertainty of the functional forms 
of equation 7-13, Qc is the uncertainty of the coeffi- 
cients of the equation, QR is the uncertainty of the 
range, etc. The values of the coefficients are presented 
in table 7-4. Assume Qr ::::: 0.2 and Qc = 0.07 for D~/1I3 


~ 0.35 slkt1l3, and Qc = 0.17 for D~/W1I3 >0.35 slkt1l3 
(from equation 2-17 where n = 2 and R = 1.3 and 
R = 1.5, respectively). 


? 


• 


+ §[CD (0.434 - log R/W1I3) _Co]2 Q~ 


(7-11) 


and for Equation 7-10: 


(7-12) 


Q2 
= Q2 + -1-1[W2/3 C2 + R2 C2] Q2 
D+ 
r 
D+2 
0 
D 
c 
u 
u 


W2/3 
f 
+__ C2Q2 
+R2C2Q2 
9 
0 
W 
D 
R 
(7-14) 
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Figure 7-19. Duration of Positive Overpressure at a Near-Ideal Surface (DNA, 1972) 
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Figure 7-20. Duration of Pusitiue Ouerpresslireai a Thermally Nonideal Surface (DNA, 1972) 
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Figure 7-21. Duration of Positive Overpressure at the Surface 
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Figure 7-22. Duration of Positive Overpressure at the Surface 
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Table 7-3. Parameters for Positive Phase Duration 


Scaled HOB 
Scaled Range 


ft/kt 1/ 3 
(m/kt 1/3) 
ft/kt 1/ 3 
(m/kt1/ 3) 
C 
CD 
0 


Near-Ideal Surface 


0 
(0) 
300 to 
750 
(91 to 229) 
1. 25 x 10-2 
2.50 x 10-4 


750 
2150 (229 
655) 
8.89 x 10-2 
1.48 
-4- 
to 
to 
x 10 


200 
(61) 
440 to 
860 (134 to 262) 
-8.08 x 10-1 
3.44 x 1C-1 


860 to 2450 (262 to 747) 
-9.12 x 10-1 
3.79 x "10-1 


400 
(122) 
450 to 
850 (137 to 259) 
-6.58 x 10-3 
2.37 x 10-4 


850 to 2075 (259 to 632) 
7.65 x 10-2 
1. 37 x 10-4 


Nonidea1 Surface 


0 
(0) 
350 to 
780 (107 to 238) 
7.89 x 10-3 
2.63 x 10-4 


780 to 1650 (238 to 503) 
1.03 x 10- 1 
1.49 x 10-4 


200 
(61) 
380 to 
500 (116 to 152) 
-9.84 x 10- 1 
4.20 x 10-1 


500 to 1300 (152 to 396) 
-5.00 x 10- 1 
2.41 x 10- 1 


1300 to 2300 (396 
701) . -1. 01 
0 
4.03 x 10-1 
to 
x 10 


400 
(122) 
150 to 
600 
(46 to 183) 
5.50 x 10- 2 
3.00 x 10-4 


600 to 2150 (183 to 655) 
1.90 x 10-1 
7.42 x 10-5 
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Figure 7~23. Duration of Positive Dynamic Pressure at the Surface 
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Figure 7-24. Duration of Positive Dynamic Pressure at a Dust-Free or Light-Dust Surface (DNA, 1972) 
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Figure 7-25. Duration of Positive Dynamic Pressure at a Heavy-Dust Surface (DNA, 1972) 
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Table 7-4. Parameters for Dynamic Pressure Positive Phase Duration 


Scaled HOB 
Scaled Range 


ft/kt l / 3 
(m/kt l /3) 
ft/k t 1/3 
(m/kt1/ 3) 
C 
Co 
0 
Near-Ideal Surface 


0 
(0) 
350 to 
550 
(107 to 168) 
O. 11 
0.00040 


550 to 1100 
(168 to 335) 
0.3] 
0.00003C, 


1100 to 2150 
(335 to 655) 
0.20 
0.00013 


200 
(61) 
450 to 
600 
(137 to 183 ) 
0.0] 
. 0.00053 


600 to 1100 
(183 to 335) 
0;31 
0.000040 


1100 
2450 
(335 
747) 
·0.24 


.- 
O.OOOlD 
to 
to 


400 
(122 ) 
450 to 
650 
(137 to 198) 
0.07 
0.00040 


650 to 1250 ,(198 to 381) 
U.31 
0.000033 


1250 to 2000 
(381 to 610) 
0.22 
0.00011 


2000 to 3100 
(610 to 945) 
0.32 
0.000055 


Nonideal Surface 


0 
(0) 
350 to 
500 
(107 to 152) 
0.063 
0.00053 


500 to 
975 
(152 to .297) 
.' 0.31 
0.000042 


975 to 1650 (297 to 503) 
o..2 3 
0.00012 


200 
(61) 
400 to 
500 
(122 to 152) 
0.29 
0.00080 


500 to 1200 (152 to 366) 
0.32 
0.000029 


1200 to 2400 
(366. to 732) 
0.27 
0.000067 


400 
(122 ) 
150 to 
300 
(46 to' 91) 
0.17 
0.00053 


300 to 2150 
(91 to 655) 
0.3] 
O. 00005'~ 


U.S:- Army Corps of 
Eng~neer5 
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7-6. Pressure waveforms. 


a. General. Under ideal conditions, i.e., airblast pro- 
pagating across snow, ice, or water, the overpressure 
and dynamic waveforms can be characterized by the 
curves presented in figures 7-26 and 7-27, where the 
peak pressure and posjtive phase durations are ob- 
tained from the data presented in paragraphs 7-2 
through 7-5. Under nonideal conditions, such idea- 
lized waveforms will not a,pply, since a number offac- 
tors act to distort and modify the blast shape. Among 
these are: 
- Location of the precursor formation region 
- Height of a structure relative to the location of 
the regular and Mach reflection regions 


7-30 


- Effect of rain, fog, or snow 
- Terrain features 
- Weapon design 
This document addresses only the precursor forma- 
tions since this probably will predominate for surface 
structures. Most structures will be located below the 
triple point and will be designed for worst-case con- 
ditions (clear day) over the actual terrain. The terrain 
(topography) must be considered in determining the 
actual waveforms. Exposed structures locally convert 
free-field dynamic pressure to overpressure; the 
degree and areal extent of conversion depend on size 
and orientation of the structures (DNA, 1972; GE- 
TEMPO, 1972 reports.) 
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Figure 7-26. Positive Overpressure Waveforms for an Ideal Shock Wave in a Standard Sea Level Atmosphere (DNA, 1972) 
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Figure 7-27. Positive Dynamic Pressure Waveforms for an Ideal Shock Wave in a Standard Sea Level Atmosphere (DNA, 1972) 
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Q2= Q2+ l.-j[C2W 2/3 + C2R2+ C2W2/3 exp(2C RlWlIS) 
I 
[ 
12 
1 
2 
3 
S 


+ C2W2/S exp(2C RIW1/3)] Q2 
4 
6 
c 


(7-16) 


+ R2 [C2 + CsCs exp(CsRlWlIS) 


+ C4C6 eXP(C6RlWlIS)] 2Q~ ~ 


near-ideal and nonideal surfaces, respectively, or from 
figures 7-33 and 7-34, or from equation 7-15 in con- 
junction with the data in table 7-5. It is recommended 
that the mean values be obtained from figures 7-31 
and 7-32 by interpolation. 
c. The uncertainty of the overpressure impulse cor- 
responding to equation 7-15 is 


b. Precursor effect-overpressure. Waveforms are 
classified in terms of the state of precursor formation 
as the five general types shown in figures 7-28. The 
type of waveform for a given situation is to be 
estimated from the data presented in figure 7-29. The 
peak values of the pulse, the positive phase duration 
and the impulse should be obtained from the data 
presented in paragraphs 7-2, 7-4, and 7-7 (GE- 
TEMPO, 1972). 
c. Precursor effect-dynamic. pressure. 
Dynamic 
pressure waveforms are also categorized as the five 
general types as shown in figure 7-30. The ap- 
propriate type for a given problem is estimated from 
figure 7-29. Actual peak dynamic pressures should 
be obtained from paragraph 7-3, and the positive 
phase should be obtained from paragraph 7-5 (GE- 
TEMPO, 1972). 
. 


7-7. Overpressure impulse. 


a. Overpressure impulse at the surface is primarily 
a function ofslant range. Reasonably accurate equa- 
tions of the impulse data presented in figures 7-31 
and 7-32 (which show parametric representations of 
impulse as functions of HOB and ground range) can 
be represented by 


_1_ = C 
+ C ~ + C 
exp (c ~) 
W1I3 
1 
2 W1I3 
3 
S W1I3 


+ C4 exp (C6 W~/3) 
(7-15) 


• 


• 


where I is the impulse in psi-s, W is yield in kt, and 
R is range in ft. The symbols C1 through C6 denote 
constants. 
b. The mean values for overpressure impulse can 
be obtained directly from figures 7-31 and 7-32 for 


where Q[ is the uncertainty of equation 7-15, Qc is 
the uncertainty of the coefficients, QR is the uncer- 
tainty of the grounl;! range, and so forth. Assume 
Q[ = 0 and Qc = 0.1 (using equation 2-17 where 
n = 2 and R = 1.5). 


• 


7-33 


TM5-858-2 


TYPE I 


PREPRECURSOR 


TYPE II 


STAGES OF DEVELOPMENT OF PRECURSOR 


TYPE III 


• 


• 


TYPE 
IV 


TYPE V 


STAGES OF EVOLUTION OF FULLY DEVELOPED PRECURSOR 


STAGES OF DECLINE OF PRECURSOR 


POSTPRECURSOR 


[ ----- 


Figure 7-28. Typical Sequence of Overpressure Waveforms over Thermally Nonideal Surfaces (DNA, 1972) 
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Figure 7-29. Type of Overpressure Waveform Developed over a Thermally Nonideal Surface (DNA, 1972) 
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Figure 7-30. 
Typical Sequence of Dynamic Pressure Waveforms over Thermally Nonideal Surfaces (DNA, 1972) 
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Figure 7-31. Positive Overpressure Impulse at a Near-Ideal Surface (DNA, 1972) 
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Figure 7-32. Positive Overpressure Impulse at a Thermally' Nonideal Surface (DNA, 1972) 
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Figure 7-33. Positive Overpressure Impulse at a Near-Ide'al Surface different burst heights 
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Figure 7-34. Positive Overpressure Impulse at a Nonideal Surface different burst heights 
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Table 7-5. 
Values of Constants for Overpressure Impulse Functions 
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Scaled HOB 
Constants 
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C4 
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CHAPTER 8 
AIRBLAST-INDUCED GROUND SHOCK 


or 


or 


(8-1a) 


(8-2a) 


~ 


Co 
avyU2 _C~ 


aH ~ 
av 
Co > U!l{2 
(8-1bJ 


where V = overhead airblast shock speed, ft/s (see fig. 
29 in Brode, 1964) 


b, Peak soil velocity. Peak vertical and horizontal 
velocities are represented by 


k zn 
av = _0_10A + B log R1W1/3 + C lolf R1W 1/3 


CoY 


8-3. Superseismic peak soil response. 


a. Peak soil acceleration. Peak horizontal and vertical 
accelerations are directly related to the peak over- 
pressure at the surface, modified by the soil proper- 
ties and the depth. The vertical component is given 
by Sauer (1964) and Crawford et al. (1974) as 


8-2. Ground Response. 


For convenience, airblast-induced ground response is 
usually segregated into two regimes-superseismic 
and outrunning. In the superseismic region, the speed 
of the overpressure shock front is greater than the 
shock speed in the ground. The ground-shock signal 
appears to be attached to and dragged along by the 
overpressme shock front traveling across the surface. 
In the outrunning region, the ground-shock signals 
arrive first, "outrunning" the overpressure shock 
front. 


where 
av 
Vertical soil acceleration, g 
Pso 
Overpressure, psi 
Co 
Soil seismic sound speed, ftls 
y 
In situ soil density, Ib/ft3 
Z 
Depth, ft 
R 
Range, ft 
W 
Yield, kt 
A, B, C, n, ko = Constants 
The horizontal component of soil acceleration may be 
taken as 


8-1. Gross phenomenology. 


a. At near-surface locations, airblast-induced ground 
shock can cause outright structmal failme offacilities 
or the failure or malfunction of internal systems. To 
achieve survival, whether by burial or other harden- 
ing techniques, the facility designer must have 
knowledge of the intensity of the ground-shock en- 
vironment and of the waveforms. 
b. The basic variables affecting the intensity ofthe 
ground-shock environment are: 
- Weapon design 
- Weapon yield 
- Burst elevation 
- Slant range 
- Earth material properties 
The larger the yield, the lower the HOB, the smaller 
the ground impedance, and the closer the facility to 
the surface, the more severe will be the airblast- 
induced ground-shock environment. 
c. The airblast signal produced by a near-surface 
nuclear explosion causes ground-shock signals to be 
formed as the shock front expands. Near the surface 
in the high overpressure region, the ground-shock 
signal is a sequence of events represented as the com- 
bination of (1) a downward, outward motion from in- 
itial loading by the overhead airblast overpressure, 
and (2) a later upward, outward motion from release 
of the overhead overpressure and from motion in- 
duced by overpressure closer in but delayed by trans- 
mission through the ground. The first resembles the 
airblast overpressme in pulse shape and dmation, the 
pulse shape modifying with depth. The second con- 
tribution, termed a "run-up" contribution because the 
effects are delayed by travel time through the 
medium, comes from the intense overpressure occur- 
ring in the region near GZ and the subsequent com- 
plex of motions known as directly induced ground 
shock. 
d. At more distance points, outrunning motion 
occurs in which the signal in the ground outraces the 
air shock on the surface. The motion is characterized 
by an upward and outward surge and often exhibits 
an oscillatory nature. In addition to these contribu- 
tions, there may also be crater-induced motion, which 
is most pronounced near the crater and diminishes 
Iwith increasing range. 
e. In the region between the ground surface and a 
depression angle of about 30 deg below the surface, 
the response is arbitrarily termed airblast-induced 
ground shock regardless of the origin of the motion. 
See figure 9-1 to visualize delineations made in this 
volume. 


• 


• 


8-1 


and 


". 
• 
Z > 20 ft 


1 ~ Z ~ 20 ft 


10.64 
-4.71 
0.487 


A' 
B' 
C' 


{ 


-0.4 
., 


n = 
-1.48, 


8-3 and 8-4), and peak stress (eqs. 8-5), use the 
following constants to calculate mean v~llles: 


A 
=·11.87 
B 
== 
-5.26 
C 
= 
0.543 
D 
= 
2.7~ x 105 
m 
= -2.09 
In calculating uncertainties at f, g, h, and i below, find 
Qpso from equation 7-2 and assume the following 
value: 
Qf 
:::: 0.13 


Qg 
:::: 0.09 
Qh 
::: 0.24 
Qk 
:::: 0.13 
Qm 
:::: 0.10 
f 
Uncertainty ofpeak soil accelerations. Using the 
values from e above, calculate uncertainties of ver- 
tical and horizontal accelerations: 


{ 


1.72 x 105 
k = 
o 
4.36 x 106 
(8-2b) 


Vertical soil particle velocity, ft/s 
Horizontal soil particle velocity, 
ftls 


1115 Z + 11 CLWl/3 


Compressional wave speed of 
soil under load, ftls 
A', B', C', D, m are constants 
Equations 8-2 are adaptations of the work presented 
by Sauer (1964) and Murrell (1974). Note that the 
horizontal particle velocity does not attenuate with 
depth. 
c. Peak soil displacement. The peak vertical and 
horizontal displacements are (Sauer, 1964) 
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Uv = 0.11 f(Z,W,CL) lOA' + B' 10gR~1I3 + C' 10g2R/W1/3 


where 


Uv 
uH 


log Pso = A + B log RIWl/3 + C log2 RIWl/3 (8-4) 


d. Peak soil stress. The peak stress in the soil is 
computed from equations 8-2 via the relationships 


where the variables have been previously defined. 
Equations 8-3 may be converted to range via the 
relationship 


Q2 + Q2 
+ Q2 
+ Q2 + n2Q2 
f 
Pso 
Co 
y 
Z 


Q2 
= 
Q2 + Q2 
+ Q2 + n2Q2 
• 


av 
f 
Co 
y 
Z 


+ [B + 2C log RIWl/3r[~ + Q~] 
(8-6a) 


Q2 + Q2 
+ Q2 + n2Q2 + 
1 
f 
Pao 
y 
Z 
(U2 _ C2)2 


0 


[U4Q2 + C4Q2 ] 
C ~U/V2 
U 
0 co' 
0 


Q2 
or 
= 
aH 
Q2+ Q2 + n2Q2+ 
1 
[U4Q2 + C4Q2 ] 
f 
y 
Z 
(U2 _ C2)2 
U 
0 
Co 
0 


+ (B + 2C log RIWl/3f (~+ Q~), 


Co~UfV'2 


Q2 
, Co >U;y2 
(8-6b) • 


av 


(8-5a) 


(8-5b) 


(8-3b) 


where 


eCLUV 
. 
0v = --, pSI 
144 


v 
. 
0H = - 
'Iv, pSI 
I-v 


dH = dv/3, in. 


1 


0 4Wl/3 
Pso - 40 
l/3~' 
dv = f(Z,W,CL) 13 Psc. -- + 280 
W, m. 
. 
C 
C2 
o 
0 
. . 
(8-3a) 


y lb sec2 


'1= g 
ft4 
g = 32.2 ft/sec2 


v = Poisson's ratio 
e. Constants in mean and uncertainty calculations 
ofpeak soil response. For peak accelerations (eqs. 8-1, 
peak velocities (eqs. 8-2), peak displacements (eqs. 


8-2 


(8-7b) 


(8-8b) 


(8-7a) 


(8-10b) 


(8-11b) 


(8-10a) 


(8-11a) 


and 


[ 
W4/3 
ZW] 
dv = Bv---B'v- 
R3 
R3 


[ 
W4/3 
ZW] 
dH = 
BH--- B'H- 
R3 
R3 


b. Peak rock velocity and displacement. Utilizing an 
adaptation of work performed by Cooper (1973), the 
peak horizontal and vertical particle velocities. and 
displacements can be represented by equations of the 
form: 


where Av, AH, Av,Ait, Bv, BH , Bv,and Bit are con- 
stants and the other variables have definitions as 
previously defined. 


c. Peak rock stress. The peak horizontal and vertical 
stresses are identical to the definitions presented in 
equation 8-5. 
d. Mean peak 
rock 
response. 
The mean peak 
horizontal and vertical accelerations are evaluated by 
using equations 8-1 (variables redefined for rock 
where 
necessary) with the values presented in 
paragraph 8-3e. The corresponding velocities and 
displacements are obtained from equations 8-10 and 
8-11 where the values ofthe constants from the data 
presented in table 8-1 should be used. The mean peak 
horizontal and vertical stress may be obtained from 
equation 8-5 where the proper mean values for the 
properties of particular rocks are considered. 
e. Uncertainty afrack response. Equations 8-6 and 
8-9 may be utilized to compute the uncertainty of 
rock acceleration and stress. The values of Qf and Qrn 
presented in paragraph 8-3e are still applicable. The 
computation of velocity and displacement uncertain- 
ties produces 
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Q2 
= Q2 
dH 
dV 


These equations may be expressed in terms of range 
and yield by utilizing equation 7-2. 
i. Uncertainty of peak soil stress. Using values in 
e above, calculate uncertainties of vertical and 
horizontal peak stresses: 


[ 
2230 Z + 11 C W1I3)2 Q2 
+ 
L 
~ 


1115 Z + 11 CLW1I3 
9 


Q2 + [ 
1115 Z 
]2[~ + Q2 + Q2] 
g 
1115 Z + 11 CLW1I3 
9 
CL 
Z 


+ 0.803Q~ 


so 


[ 
5.2 J1lo4 Co + 280 Pso 
]2 
2 
+.4 
Qp 
. (8-8a) 
13 ~o Co + 280(Pso - 40) 
so· 


Q2 
= 
Q2 + [ 
1115 Z 
]2 [Q2 
+ Q2] 
uv 
g 
1115 Z + 11 CLW1I3 
cL 
Z 


+ [B' + 2C' log R!W1/3rQ~ 


r 


g. 
Uncertainty ofpeak soil velocities; Using values 
from e above, calculate uncertainties of vertical and 
horizontal peak velocities: 


+ [ 
1115 Z 
_ B' 
1115 Z + 11 CLW1I3 


r- 2C'log RiW1I3rQ~ 


Q2 
= Q2 + m2[Q2 + ~] 


UH 
h 
R 
9 


h. Uncertainty of peak soil displacements. Using 
values from e above, calculate uncertainties of ver- 
tical and horizontal peak displacements: 


Q2 
Q2 + [ 
1115 Z 
]2 [Q2 
+ Q2] 
dV 
k 
1115 Z + 11 CLW1I3 
cL . 
Z 


.. 


.. 


.. 


• 


• 


8-4. Superseismic peak rock response. 


a. Peak rock acceleration. The peak horiontal and ver- 
tical accelerations in rock are described by equations 
8-1. 
• 


(8-9a) 


(8-9b) 
(8-12a) 


(8-12b) 


8-3 


TM 5-858-2 


Table 8-1- Mean Parameters for Calculation of Velocity and Displacement Due to Airblast on Rock 
- 


~ 


. ' 144 
A' 
0 
o < Z/Wl / 3 <' 31 ft/ktl / 3 • 
= 
, 
= 


~424· 
H 
Z/Wl (3 -:- 61 ft/kt l / 3 
2.16 
, .' 
31 
< 


I 


Av 
1193 
A " 
·1.57 
0 
< Z/Wl / 3 
< 31 ft/kt l / 3 
= 
= 


) 0.62 
~134 
V· 
Z/Wl / 3 -:- 
ft/kt l / 3 
31 
< 
61 


BH 


f5.03 
X 106 
B ' 
0 
0 
< Z/Wl /3 
< 31 ft/kt l / 3 
= 
= 
\ 
106 
H 
104 
Z/Wl / 3 < 61 ft/kt l / 3 


~1 
17.82 
X 
., 8.99 
X 
'31 
< 


1'06 
/ 
) 


BV 
15.03 
X 
B 
I 
18.99 X 104 
o < Z/Wl / 3 
< 31 .ft/kt1/3 
= 
= 
, 
'3.35 
X 106 
V 
X 104 
'Z/W1/ 3 < 
ft/kt 1/ 3 
13.59 
31 
< 
61 


where Qp :::: 0.36 and Qq :::: 0.36. 


8-5. Superseismic coil and rock waveforms. 


a. Soil waveforms. According to the data presented by 
Murrell (1974), vertical velocities at a constant range 
maintain a relatively constant waveform but decrease 
in amplitude with depth. Conversely, at a constant 
depth, the vertical velocities undergo distinct changes 
in both waveform and amplitude with variations in 
range. Horizontal velocities qualitatively display the 


u.s. Army Corps of Engineers 


[ 
B' ZW12 
Q2 
= Q2 + 9Q2 + 4 + _Y_ 
Q2 
Dv 
q 
R 
R3dV 
w 


+ [B'YZW]2 Q2 
R3dv 
Z 
(8-13a) 


(8-13b) 


same type of behavior although they are markedly 
different in character than their vertical counterparts. 


(l) Typical examples taken from Murrell (1974) 
have been normalized and are presented at various 
ranges for the vertical velocity in figure 8-1 and for 
the horizontal velocity in figure 8-2. The data 
generally apply at all depths in the superseismic 
region. 
(2) Parenthetically, it is noted in figure 8-2 that 
outrunning motion appears at the 340 ftlkt1l3 range 
and is characterized by the rolling motion after a time 
of about 1.0 slkt1l3. The outrunning is well developed 
at the 535 and 680 ftlktl/3 ranges. The amplitudes 
presented in paragraph 8-3 apply to the motion prior 
to the time of the arrival of the outrunning motion. 
After the outrunning motion occurs, the amplitudes 
of paragraph 8-6 will apply. The reader is encouraged 
to 
study 
Murrell 
(1974) prior 
to 
constructing 
waveforms having the generalized waveform charac- 
teristics of figures 8-1 and 8-2 but adjusted to match 
the amplitude statistics presented in paragraphs 8-3 
and 8-6. 
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Figure 8-1. Vertical Velocity Waveforms for Soils at All Depths(A dapted from Murren 1974) 
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Figure 8-2. Horizontal Velocity Waveforms for Soils at All Depths (Adapted from Murrell, 1974) 
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8-6 


• 


b. Rock waveforms. Utilizing normalized data from 
Murrell-Carleton (1973), figures 8-3 and 8-4 present 
selections of vertical velocity histories at various 
scaled ranges but at the constant specified scaled 
depths. It is noted that at a given range, the depth 
influences the character ofthe waveforms (in contrast 
to the data presented in a above in which depth plays 
a minor role in modifying velocity signatures in soil). 
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Conversely, figure 8-5 presents normalized horizon- 
tal velocity histories for a selection of scaled ranges 
and for all depths of interest. The reader should study 
Murrell-Carleton (1973) prior to adapting the data in 
figures 8-3 through 8-5 to suit his own needs. 
Derived signatures should be compatible with the 
peak amplitude data presented in paragraph 8-4 
above for airblast-induced response in rock. 
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Figure 8-3. Vertical Veloczty Waueforms for Rock, Scaled Depth of 4.3 [tIkt1l3 0.31 m1ktl /3) (Adapted from Murrell-Carleton, 1973) 
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Figure 8-4. 
Vertical Velocity Waveforms for Rock, Sc'aled Depths Below 38.8 ftlkt?13 (11.83 mlkt1l3) (Adapted fro~ Murrell-Carl~ton,1973) 
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Figure 8-5. Horizontal Velocity Waveforms for Rock for All Depths (Adapted from Murrell-Carleton, 1973) (l of2) 
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Figure 8-5. 
Horizontal Velocity Waveforms for Rock for All Depths (Adapted from Murrell-Carleton, 1973) (2 of 2) 
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8-6. Outrunning peak ground response. 


a. Peak ground acceleration. 
Peak vertical and 
horizontal outrunning accelerations (adapted from 
Sauer, 1964) are given by: 


TM 5-858-2 


1 


-3.5, 150 ~ R!W1I3 <::; 800 , 1 <::; Z ~ 16 
-3.5, 150 <::; R!W1I3 <::; 800 , 
Z > 16 
m = 
-2.0, 800 < R!Wl/3 <::; 3000, 1 <::; Z <::; 16 


. 
-2.0, 800 < R!W1I3 <::; 3000, 
Z> 16 


where 
av 
Vertical acceleration, g's 
aH 
Horizontal acceleration, g's 
R 
Range, ft 
W 
Yield, kt 
Z 
Depth, ft 
Co 
Seismic velocity, ftls 
ko' n, m are constants 
b. Peak ground velocity. Equations for the calcula- 
tion ofpeak vertical and horizontal oscillating outrun- 
ning particle velocities designated as urn and ap- 
plicable at all near-surface depths of interest in design 
can be represented by the equation 
(8-17b) 


D = 6400 
s = -1.98 
The mean peak values for the paramet~rs in equa- 
tion 8-15 are presented in table 8-2 (Higgins, 1975). 
For convenience, the data in table 8-2 are cross 
referenced to the appropriate figure (Higgins, 1975). 
The ratio of the initial peak, ul' to the oscillatory 
peak, urn' is presented in figure 8-6 as a function of 
scaled range. 
e. Uncertainty ofpeak acceleration. The uncertain- 
ties of equation 8-14 for peak outrunning accelera- 
tions are: 


(8-14b) 


(8-14a) 


• 


(8-15) 


where 
u 
= Peak oscillatory outrunning velocity, 
m 
ft/s 
c. Peak ground displacements. Peak vertical and 
horizontal out-running displacements for all near- 
surface depths (adapted from Sauer, 1964) are given 
by 


where Qf :::: 0.8. 
f 
Uncertainty of peak velocity. The uncertainty 
associated with the amplitude of the oscillatory com- 
ponent ofthe vertical outrunning ground shock or for 
the maximum horizontal motion for contact or near- 
surface bursts is 


(8-18) 


where 
dv = Vertical displacement, in. 
dH = Horizontal displacement, in. 
D, s are constants 
d. Mean peak ground response. Equations 8-14 and 
8-16 may be evaluated for their mean values from 
the Quantities: 


where Q 
is the uncertainty of equation 8-15. The 
values ol'Qg are presented in table 8-2 in which the 
data scatter in (Higgins, 1975) is assumed to be log- 
normally distributed. The uncertainty, QUI' of the in- 
itial peak vertical velocity ul (depicted in figure 8-7) 
is given by 


(8-19) 


(8-20b) 


(8-20a) 


QOH = Qov 


where Qk :::: 0.6. 


where the values of Qh are presented in table 8-3 
from data presented in (Higgins, 1975). 
g. Uncertainty ofpeak displacement. The uncertain- 
ties of the displacements in equations 8-16 are 


Q~v = Q~ + s2Q~ + (1~S)2 Q~ 


(8-16a) 


(8-16b) 


dvW 1l3 = D(R!Wl/3)s 


3 X 1014 
-0.44 


3.75 x 1015 
-1.35 
k = 
n= 
° 
1.5 X 1010 
-0.44 
• 


1.88 x 1011 
-1.35 
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CD' 
I- 
~. 


Table 8-2. Values of Parameters for Calculating Outrunning Mean Ground Shock Velocity and Associated Uncertainties (Adapted from 
Higgins, 1975) 


Direction 
Material 
Coefficient 
Exponent 
Higgins 
FigU~e No. ' 
Uncertainty 
C 
r 
in 1975 Reference 
Slg, 


Wet clay/shale 
3 
-1. 2 
32 
1.33 x 10 
0.36 


Dry clay/shale 
2 
-1.05 
34 
0.60 
3.91 
x 10 


Vertical U 
Sand/sandstone 
2 
-1. 05 
35 
0.48 
m 
3.91 x 10 


, 


Tona1ite 
3 
-1.2 
33 
0.60 
1.15 x 10 


Clay 
3 
-1. 35 
36 
3.08 x 10 
0.40 


Wet c1 ay/ shale 
4 
5.41 x 10 
·"'1.8 
50 
0.60 


Dry clay/shale 
4 
2.57 x 10 
-1.65 
52 
0.30 


Horizontal u 
Sandi sandstone 
5 
2.79 x 10 
-1. 95 
53 
0.36 
m 
4 
Tonalite 
2.69 x 10 
-1.7 
51 
0.50 


Clay 
3 
5.66 x 10 
-1.4 
54 
0.36 
. 


u.s. Army Corps of Engineers 
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Figure 8-6. Ratios of Initial to Oscillatory Peak Particle Velocities 
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Figure 8-7. Characteristic Vertical Outrunning Waueforms (Adapted from Higgins, 1975) 


• 


8-14 


• 


TM 5-858-2 


Table 8-3. Uncertainties of the Initial Peak Velocity 


Uncertainty 
Geology 
~h 


Wet clay over shale 
0.18 


Dry clay over shale 
0.37 


Sand over sandstone 
0.32 


Tonalite 
0.33 


Clay 
0.16 


U.S. Army Corps of Engineers 


(8-21) 


(8-23) 


(8-22) 


Acceleration 
for 
a 
particular 
HOB and range 
Overpressure at that same HOB 
and range, from chapter 7 
Overpressure at that same range 
for· a contact burst, from chapter 7 
Acceleration at that same range 
for a contact burst, from equation 
8-1. 


a,:HOBl 


Pso(HOBl 


where 


equations 8-1, 8-2, 8-5, and 8-10 must be adjusted 
for HOB effects. At any scaled range in these equa- 
tions, the ratio of the overpressure at the appropriate 
scaled HOB to the overpressure for a contact burst 
from chapter 7 establishes a correction factor used to 
adjust the peak velocity, acceleration, and stress. 
Thus, 


Pso!HOBl 
a(HOBl = 
P 
a(o) 
solo) 


Velocity and stress for a par- 
ticular HOB 
Velocity and stress at that 
same range for a contact 
burst using equations 8-2, 
8-5, and 8-10 
c. Displacement. The peak displacement data of 
equations 8-3 and 8-11 also must be adjusted for 
HOB effects. At any scaled range in these equations, 


and 


Pso(HOB) 


u(HOBl = 
P 
u(o) 
so(o) 


Pso(HOBl 
0(HOB) = 
P 
0(0) 
solo) 


where 


u(HOB), °(HOBJ 


8-8. Height of burst effects. 


a. Correction factors. Detonation at positive heights 
of burst will alter overpressure and initial coupling 
of energy to the ground, hence the airblast-induced 
ground shock in the superseismic region. To a first 
approximation, 
HOB-corrected 
accelerations 
and 
velocities can be obtained by applying HOB peak 
overpressure correction factors; displacements can be 
modified by applying the HOB-correction factors for 
the overpressure positive-phase impulse. 
b. Velocity, acceleration, and stress. The velocity, 
acceleration, and stress amplitudes characterized in 


8-7. Outrunning waveforms. 


a. Velocity. Outrunning is a complicated phenomenon 
that is not easily described by simple waveforms. 
Nevertheless, over the years, certain waveform 
characteristics have been identified that describe the 
principal elements of the velocity signature. The 
reader will be required to adapt the signatures recom- 
mended below 
to 
the 
amplitudes 
presented 
in 
paragraph 8-6. 
b. Vertical. 
Figure 8-7(a) presents the recom- 
mended waveform for vertical outrunning motion in 
layered soft rock and soil sites. Figure 8-7(b) shows 
the recommended waveform for homogeneous sites. 
c. Horizontal. Horizontal velocity waveforms vary 
continuously with range, and, near the crater, with 
depth as well. Composite waveforms in figure 8-7 are 
more difficult to develop. Figure 8-8 presents typical 
waveforms for wet and dry clays over shale, and sands 
over sandstones. These data are presented for three 
crater-volume scale ranges based on apparent crater 
volumes. Refer to (Higgins, 1975) for definition oftime 
scales. The reader may also obtain some guidance 
from paragraph 8-5a on how outrunning develops. 
Consult chapter 12 for the computation of crater 
volumes. 


• 


• 
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the ratio ofthe impulse at the appropriate scaled HOB 
to the impulse for a contact burst establishes a factor 
that is to be used to adjust the peak" displacement. 
Thus, 


ImOB) 


~HOB) = -- ~o) 


1(0) 
(8-24) 


~o) 


Impulse at that same HOB and 
range, from chapter 7 
Impulse at that same'range for a 
contact, burst, from chapter 7 
Displacement at that same range 
for a contact burst, from equations 
8-3 and 8-11 
• 


where 


~HOB) 
Displacement at a particular HOB 
and range 
The ratio of I(HOB}!ICo) as a function of range and HOB 
is presented in table 8-4. 
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Figure 8-8. Normalized Horizontal Velocities for Clay over Shale an'd Sand over Sandstone Sites (Adapted from Higgins, 1975) 
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, ,Table 8-4, Impulse IIHOBj!(o) for Displacement Computation' 


Scaled Rim ge 
HOB 
ft/kt 1/ 3 
(m/kt 1/ 3) 


ft/kt 1/3 
(m/kt 1/3) 
200(61) 
400(122) 
600 (188) 


Near-Ideal Surface 


300 (91) 
1.2 
1.1 
1.4 
.' 


400 (122) 
1.3 
1.0 
1.1 


600 ( 183) 
1.2 
1.0 
1.0 


800 (244) 
1.2 
1.1 
1.0 


1000 (305) 
1. 1 ' 
1.0 
1.0 


" 


1250 (381) 
1.1 
1.0 
1.0 


1500 (457) 
l.0 
1.0 
1.0 


2000 (610) 
1.1 
1.0 
1.0 


Nonideal Surface 


300 (91) 
0.8 
0.8 
- 


400 (122) 
1.0 
1.1 
1.1 


600 (183) 
1.0 
1.2 
1.2 


800 (244) 
1.0 
1.2 
1.2 


1000 (305) 
1.0 
1.2 
1.2 


1250 (381) 
1.0 
1.1 
1.1 


1500 (457) 
1.0 
1.1 
1.1 


2000 (610) 
0.9 
1.0 
1.0 


2500 ( 762) 
0.9 
1.0 
1.0 


U.S. Army Corps of Engineers 


• 


d 
Uncertainty. 
For all practical purposes, the 
uncertainties associated with acceleration, velocity, 
and displacement are equivalent to those established 


for contact bursts, as defined in paragraphs 8-3f, g, 
h for soils, and paragraph 8-4e for rock. 
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CHAPTER 9 
CRATER-INDUCED GROUND SHOCK 


\ 
U.S. 
Army Corps of Engineers 


Figure 9-1. 'Regions of Ground Shock 
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.. 
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9-1. Gross phenomenology. 


a. Direct coupling of weapon energy from near- 
surface, contact, and shallow underground bursts 
drives shock waves through the ground and leads to 
crater formation. The relative importance of airblast- 
induced signals, unloading signals from the ground 
surface, and signals resulting from crater formation 
depends on the location of the facility relative to 
ground zero (GZ). As shown in figure 9-1, the region 
between 45 deg(0.79 rad) and 30 deg (0.52 rad) is in- 
fluenced by composite signals of both airblast-induced 
and crater-induced shock waves. In the central con- 
ical region below depression angles of ""45 deg. crater- 
induced signals are dominant and are the subject of 
this chapter. 
b. The basic variables affecting the intensity and 
duration of the crater-induced ground shock are: 
- Weapon design 
- Weapon yield 
- Burst elevation or burial 


CRATER- 
INDUCED 
REGION 


- Properties of the site media (rock, soil) 
- Slant range to the facility (or any other point of 
interest) 
c. Crater-induced ground shock can be predicted 
from (1) mathematical simulation, (2) scaling and ex- 
trapolation of near-surface burst data, and (3) scaling 
and extrapolation of fully contained burst data. 
Method 3 currently offers the most credible basis for 
estimating the crater-induced ground shock environ- 
ment and its uncertainty. In this method, under- 
ground nuclear ground-shock environments from 
near-surface (contact or shallow-buried) 
nuclear 
bursts are estimated from the ground-shock data of 
fully contained bursts. The regions for which these 
data apply include the crater-induced region in figure 
9-1 and, with some reservations, the composite region 
as well. 
d. Of importance to design of facilities are ground- 
shock intensity and time histories and their pulse 
shape. These features are delineated in paragraphs 
9-2 and 9-3. 


pso 


--+---tl.. U 


DEG 


AIRBLAST-INDUCED REGION 
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9-2. Ground-shock intensity. 


a. Decay curves for peak acceleration, velocity, and 
displacement are presented in the general form 
(Perret-Bass, 1975). 


values of the 
coefficients in equation 9-2 
are 
presented in tables 9-2 and 9-3. 
d. The uncertainties associated with the crater- 
induced ground shock for deeply buried systems are • 


b. In using equations9-1, which pertain to fully 
contained 
bursts, 
a 
weapon-yield 
correction 
is 
necessary to cc:mvert to.the ground-shock response for 
near-surface bursts. This correction utilizes the 
equation 
• 


(9-3c) 


(9-3b) 


(9-3a) 


where Qr are the uncertainties of the functional forms 
of the decay laws, Qc are the uncertainties of the, coef- 
ficients, and so forth. Assume Qr :::: 0.05 and compute 
QCa from equation 2-18 where L2 and L1 are the up- 
per and lower bounds on the data at any given slant 
range in figures 9-2through 9-5. Similarly, compute 
QCuandQCd from. figures 9-6 through 9-9, and 9-10 
through 9-13, respectively. Assume n = 1.65. 
e. The uncertainties associated with the yield cor- 
rec,tion presented.in equation 9-2 are 


Q2 
= Q2 + r~ l~2[1 + 2az/3W1I3 + ~z2/3W2/3 
we 
r 
W2, 


e 


(9-1c)" 


(9-1b) 


(9-1a) 


Peak 
acceleration 
(g), 
particle 
velocity (ft/s), displacement (in) 
Site geology parameters (see table 
9-1) 
Slant range from GZ (ft) 
Effective weapon yield (kt) 


( 
R )na 
Wl/3 
e 


where 
a, u, " 


C,n 


"W1I3 
e 


aWl/3 
= C 
e 
a 


u 


WWe = [1 + 
z 
z2 
z3 
aW1I3 + ~W2/3 + YW 


where We and Ware the effective and true yields in 
kt; z is depth (z > 0) or height of burst in ft; a, 
~, y, 
and ~ are coefficients of HOBIDOB influence and rrn 
accounts for scaling of deep-buried bursts to contact, 
bursts as a function of site material properties and 
the ground motion parameter of interest. 
c. The mean peak of the crater-induced ground 
shock environment is calculated from equations 9-1 
and 9-2, using values of the site-dependent coeffi- 
cients and exponents presented in table 9-1. The 


Z4 
+ ~-- + 
W4/3 


~. '.. ] rrn , _z_ < 20 ftlkt1l3 
W1I3 
(9-2) 
+ ~2z6W-2/3rQ~ 


+ z2[aW2/3 + 2~zW1I3 + 3yz2 


+ 4~z31W1/3]2 
Q2! + Q2 
(9-4) 
z 
rrn 


where Qr is the uncertainty of equation 9-2, Qw is 
the uncertainty in the yield, Qc are the uncertain- 
ties' of the 
coefficients, 
and so forth. 
Assume 
Qr = Qc = Qr 
:::: 0.2. - 
rn 
• 
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Table 9-1, Mean Values of Parameters for Calculations of Direct-Induced Ground Motion (equation 9-1) 
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Coefficient 
Exponent 
Range 


Parameter 
Geology 
C 
n 
ft/kt 1/ 3 
(m/kt1/ 3) 
a 
a 


Alluvium 
2.15 x 1014 
-5.78 
66-262 
(20-80) 


Alluvium 
6.02 x 105 
-2.13 
197-1148 
(60-350) 


aWel/3(g'kl/3 
Dry tuff 
1.. 42' x 1013 
-4.77 
131-492 
(40-150) 
Dry'tuff 
7.55 'x 105 
-1.92 
328-1641 ._(100- SO0) 


- 
108 , 


'. 


Eq. -(9-1a) 
Wet tuff,. 
9.58 :X 
-2.61 . 
98-1969 
(30-600) 


: 
Hard rock 
1. 09 
1013 
-4.35 
131-656 
(40- 200) 
x 


Hard rock 
1.46 x 108 
-2.32 
295-7218 
(90-2200) 


C 
n 
u 
·u 


Alluvium 
2.91 x 109 
-3.27 
98-492 
(30-150) 


Alluvium 
6.03 x 103 
-1. 16 
328-1148 
(100-350) 
u(ft/s) 
Dry tuff 
7.66 x 106 
-1. 98 
131-1641 
(40-500) 


Wet tuff 
1. 66 
6 
-1.56 
98-1969 
(30-600) 
x 10 . 


Eq. 
(9-1b) 
x 106 
(40- 2200) 
Hard rO.ck 
5.50 
-1. 72 
131-7218 


C6 
n6 


Alluvium 
5.02 x 107 
-3.04 
131-492 
(40-150) 


O/w1/ 3(in/kt 1/ 3) Alluvium 
3.27 x 102 
-1.11 
328-1148 
(100-350) 
e 
Dry. tuff 
2.04 x 106 
-2.20 
- 328-1641 
(l00-500) 


Eq. 
(9-1c) 
Wet tuff 
4.39 x 107 
-2.63 
164-1969 
(50-600) 


Hard rock 
3.20 x 105 
-1. 88 
230-7218 
(70-2200) 


U.S. Army Corps of EngIneers 
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Table 9-2. Factors in Calculating HOB/DOB Effective Yield (equa- 
tion 9-2) 
, 
Coefficient 
Value 


'"'( 
3.39 x 10- 2 


S 
-3.33 x 10-4 


y 
2.36 x 10- 7 


E, 
-1.19 
x 10-8 


U.S. 
Army Corps of Englneers 


Table 9-3. Factors r m to Account for Geological Properties (equation 9-2) 


Acceleration 
Geology 
and Velocity 
Displacement 


Alluvium 
0.04 
0.0025 


Tuff 
0.04 
0.0025 


Hard rock 
0.04 
0.01 


U.S. Army Corps of Englneers 
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(Perret-Bass, 1975) 
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Figure 9-4. Attenuation ofScaled A cceleration-Wet Tuff(Perret- 
Bass, 1975) 


• 
SLANT RANGE, m/kt1/ 3 
SCALED 


102 
103 


104. 
I 
I 
I IlL 
I I II 
-, 
, 
\ 
MEAN 
\, 
--- BOUNDS 
\ 
\ 
\ 
\, 


103 
\, 
\ 
\ 
\ 
, 


M 
\ 
, 
...... 
\ 
- 
\ 
\ 
.... 
.:>t. 
\ 
\ 
. 
\ 
Cl 
\ 
\ 
\ 


z 
\ 
\ 
0- 
\ 
\ 
I- 


102 
\ 
\ 
• 


<t:a: 
UJ 
, 
-J 
UJ 
\ 
ww 
\ 
<t: 
\ 
\ 
0 
\ 
\ 
UJ 
-J 
\ 
\ 
<t:w 
\ 
\ 
1/'1 
\ 
\ 


M 
\ 
\ 
....... 
\ 
\ 


-3Q) 10' 


ra 


• 


9-7 


• 
RANGE, m/kt 1/ 3 


103 


I 
I 
I 
I 
I lid 


SCALED 
SLANT 


102 


I 
I II d 


1M 5-858-2 


---MEAN 
--- BOUNDS 


C""I 
......-+J 
.::t. 


102 
. 
0'1 


Z 
0 
I- 
<:( 
a:: 
UJ 
....J 
UJ 
U 
U 
<:( 


10] 
0 
UJ 
\ 
....J 
<:( 
\ 
u 
V') 
\ 
\ 


l""'I 
\ 
........ 
\ 
Q) 
3 
\ 


rc 
\ 
10° 
\ 
\ 
\ 
\ 
\ 
\ 
\, 


] 0- 1 I-_L.-~~""""".J..L.:.._---I._.L-...I--'-'\:L....L.J""'" 


]02 
103 
104 


R/W ]13 
SCALED SLANT 
RANGE, ft/k t 1/3 
e 
' 


Figure 9-5. Attenuation ofScaled Acceleration-Hard Rock (Perret- 
Bass, 1975) 
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Figure 9-6. Attenuation of Particle Velocity-Dry Alluvium (Perret-Bass, 1975) 
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Figure 9-7. Attenuation of Particle Velocity-Dry Tuff (Perret-Bass, 1975) 
• 


9-10 


1M 5-858-2 


\ 
\ 
• 
SCALED SLANT RANGE, m/kt1/ 3 


102 
103 


104 ;":=::==;:=~I ::1:::;::::;::;:::;:::::;:::;:;==~~=:=,:I 
~I:::;:::::;:~:;:::;:~ 


Figure 9-8. Attenuation of Particle Velocity-Wet Tuff (Perret-Bass, 1975) 
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Figure 9-9. Attenuation of Particle Velocity-Hard Rock (Perret-Bass, 1975) 
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Figure 9-10. Attenuation ofScaled Displacement-Dry Alluvium 
(Perret-Bass, 1975) 
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Figure 9-11. Attenuation of Scaled Displacement-Dry Tuff (Perret-Bass, 1975) 
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Figure 9-12. AttenuatIOn of Scaled Displacement-Wet Tuff (Perret-Bass, 1975) 
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Figure 9-13. Attenuation of Scaled Displacement-Hard Rock (Perret-Bass, 1975) 
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9-3. Pulse Shapes 


a. Ground-shock 
histories 
are 
required 
for 
the 
analysis of some aspects of system survivability. For 
underground locations that are below a declination 
angle of 30 deg, it is assumed that particle velocity 
can be represented by simple triangles, as defined in 
figure 9-14a. The corresponding acceleration and 
displacement plots are shown in figure 9-14b and 
9-14c. The symbols in the figure are defined as 
follows: 


TM'5-858-2 


b. The mean pulse shape is derived by utilizing the 
above relationships and the data in table 9-1. Also 
utilize 
a = 0.694 
b = -7.05 X 10-3 
c = 3.15 X 10-5 


c. The uncertainty of a, u, and d are presented in 
equation 9-3 and in conjunction with equation 9-4. 
In addition, assume that 


Peak acceleration from equation 9-1a, g 
Peak displacement from equation 9-1c, ft 
2t+, sec 
u 
= rise time, sec 
32.2 
a 
Positive phase duration = 2d/u, sec 
Negative phase time = 4t+, sec 
Peak velocity from equation 9-1b, ft/sec 
Peak negative velocity = u/m, ft/sec 


• 


• 


• 


m= 


3 


3 


~ ~ 261 ftlkt1l3 
1/3 '" 
W 
. 
e 


~ > 261 ftlkt1l3 
W1l3 
e 


(9-5) 


Q2 
= Q2 
tf 
t+ 


Q2 
= Q2 + Q2 
, tr 
u 
a 


Q2 
= Q2 
t_ 
t+ 


Q 
:::: 0 2 + Q2 + 1 Q2 
m 
. 
R 
9" we 
(9-6) 
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Figure 9-14. Pulse Shape Parameters 
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CHAPTER 10 
CRATER 
• 


• 


• 


10-1. Gross phenomenology. 


a. A crater is formed when a nuclear weapon is 
detonated below ru 35 ft/ktl/3 (ru 10.7 m/ktl/3). The 
significant features ofthe crater include the apparent 
crater (i.e., the visible contour), the true crater as 
defined by the limits of dissociation of the natural 
material, ejecta that falls back into the true crater, 
a rupture zone in which cracking and crushing occurs, 
a plastic zone (or displacement zone in rock) in which 
significant permanent deformation occurs, and the 
ejecta that, along with ground deformation, form the 
crater lip. These zones are depicted in figure 10-1. 
b. The basic variables affecting the cratering 
phenomenon are: 
- Weapon design 
- Weapon yield 
- HOBIDOB 
- Site condition (stratification, water-table eleva- 
tion, etc.) 
- Material properties (principally shear strength) 
The greater the yield, the greater will be the DOB 
(or the lesser the HOB), or the weaker the site 
medium, the larger the crater volume. The crater will 
be flatter the higher the burst or the weaker the site 
medium. Crater size increases with DOB up to the 
optimum cratering DOB (::::131 ft/ktl/3, or 40 m/ktl/3). 
Further increases of DOB reduce crater size. 


c. Currently, crater dimensions are estimated by 
scaling from the cratering data obtained during the 
atmospheric nuclear tests (perhaps a score of crater- 
ing events) and from the large number of atmospheric 
HE tests conducted with charges weighing from a few 
grams to 500 tons. The credibility of mathematically 
simulated craters has not been generally established; 
the prediction methods described in this volume do 
not include these data. Moreover, fairly recent data 
from the AEWES, DNA-PR-0012, Phase 2 and DNA- 
PR-0028, Phase 3 (1976, 1977), which included crater- 
ing experiments in soils with a high water content, 
have not been included in this manual since the 
results were reported after the manual was prepared. 
d. The volume of material excavated by a nuclear 
blast is of little direct concern to designers except as 
it relates to the formation of ejecta material which 
produces an impact phenomenon and buildup of ejec- 
ta on the ground surface beyond the crater lip. It also 
contributes to the fo~mationof the dust cloud. These 
subjects are discussed in Chapters 11, 12, and 13. A 


detailed derivation of crater volume is presented in 
Chapter 12. A system cannot survive if its original 
position places it within the crater proper or even 


. somewhat beyond the crater. Moreover, only very sim- 
ple, compact, inherently rugged facilities can be ex- 
pected to survive in the contiguous rupture zone. The 
crater (crater proper and crater rupture zone) pro- 
duced by a single, relatively small-yield weapon poses 
a threat to a near-surface facility. The crater produc- 
ed by a single super-yield weapon, a laydown pattern 
of weapons, or a sequential detonation of weapons 
delivered to a single point poses a threat to a deep- 
underground facility. The facility designer will thus 
need to know the dimensions and uncertainties ofthe 
true crater, the rupture zone, and the plastic zone 
under the three patterns of attack. 


10-2. Single burst. 


a. 
True crater. The true crater is estimated from the 
apparent crater dimensions shown in figures 10-2 
through 10-5 for a variety of soils and rocks. The 
generalized equations for the apparent crater dimen- 
SIOns are 


Ra = a'Wl/3 + {3'z + Y'Z3W_2/3~ 
z 
ft, -- < 16.4 
Da = aWl/3 + {3z + yz3W-2/3 
WI/3 


and 


Ra = A'WI/3+ B'z + C' z2W-l/3 ~ 


ft, 16.4 <_z_· < 100 
D = AWl/3 + Hz + Cz2W-l/3 
Wl/3 
a 
/ 
(10-1) 
and the true crater dimensions are 


(10-2) 


where R is radius, D is depth, z is HOBIDOB (all in 
feet, positive for DOB and negative for HOB), and W 
is weapon yield (in kilotons). The quantities a, {3, y, 


~, A, B, and C are defined below. 
(1) The mean true crater radius and depth are 
calculated from equations 10-1 and 10-2 using the 
values presented in table H)-I. 
. 


10-1 


- 
o 
~ 


ORIGINAL 
GROUND 


LEGEND: 


~..,- ~jj EJECTA 
~ FALLBACK 


ITIIIIIIIIII UPTHRUST 
I 
I LIMIT OF APPARENT RUPTURE 


~ LIMIT OF SEVERE RUPTURE 


[:/jJ{tfJ LIMIT OF PLASTIC 
DEFORMATION OR 
DISPLACEMENT ZONE 


DP 


PLASTI C-ELASTI C 
ENVELOPE 


ELASTI C ZONE 


ZONE OF PLASTI C 
DEFORMATION OR 
DISPLACEMENT 


CHARGE 


Da 
...L Dt 
~ 


.... 
3: 


CI'I 
I00 
CI'I 
GO 
It.) 


Figure 10-1. 
Typical Crater from a Near-Surface Explosion (Rooke et ai., 1974) 
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Figure 10-2. 
Apparent Crater Dimension for Dry Sailor Soft Rock (";;'10% Water Content), (DNA, 1972) 
• 


10-3 


TM 5-:-858-2 


DEPTH OF BURIAL, m/kt1/ 3 


-10 
0 
10 
20 
30 


240 ::'=1;=:~=:;=~I=~~=:::;::I :::::;:::::;::::::;;::::1::::;:::=;:::::::::=1~:=;:~ 
• 
70 


10 


60 


401-----I--H'------+------+---'-----i 


200 


M 
........ 
M 
... 
........ 
..... 
50 
........ 
.::t/. 


R /W1/ 3 
........ 
160 ---- 
.::t/. 
':-,' 
..... 
= 
........ 
4- 
a 
E 


~ 
V'l 
V'l 
Z 
. Z 
0 
0 


V'l 
40 
V'l 
z 
z 
UJ 
UJ 
.X: 
120 
x: 


0 
0 


0::: 
0::: 
UJ 
UJ 
I- 
D /W 1/ 3 
I- 
et: 
30 


et: 
0::: 
= 
0::: 
• 


W 
a 
.W 
l- 
I- 
Z 
z 
UJ 
80 


UJ 
0:::. 
0::: 
et: 
et: 
0- 
0- 
0- 
0.. 
et: 
20 
et: 


Z/W 1/ 3, DEPTH OF BURlAL, ft/kt 1/ 3 


Figure 10-3. 
Apparent Crater DimensIOns for Wet Sailor Soft Rock (>10% Moisture Content), (DNA, 1972) 
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Figure 10-4. 
Apparent Crater Dimensions for Dry Hard Rock (';;3% Moisture Content), (DNA, 1972) 
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Figure 10-5. 
Apparent Crater Dimensions for Wet Hard Rock (>3% Moisture Content), (DNA, 1972) 
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(2) The uncertainties associated with the true 
crater dimensions are 
• 


Q2 
= Q2 + Q2 + Q2 
Rt 
g 
~ 
Ra 


Q2 = Q2 + Q2 + Q2 
Dt 
g 
~ 
Da 


where 


+ z2 ((3' + 3y'z2W-2/3)2 Q2z 


W2/3 
t 
+ -g- (a' - 2Y'z3W-l)2 Q~\ 
Iw:131 < 164 


(10-3) 


(10-4a) 
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(2) The uncertainties associated with the rupture 
zone are 
Q2 
Q2 + 
Q2 
+ 
Q2 + 
Q2 
RR 
g 
9 
~ 
Ra 


Q2 
Q2 + Q2 
+ Q2 + 
Q2 
(10-6) 
DR 
g 
9' 
~ 
Da 


where Qe, Qe', and Q~ are the uncertainties of the 
coefficients in equation 10-5, and QRa and QD are 
. 
a 
defined III a (2) above. Assume that Q2 + Q2 + Q2 = 
Q2 + Q2 + Q2 :::: 0 05 
g 
9 
~ 
g 
9' 
~ 
.. 
c. Plastic or displacement zone. The plastic zone (for 
soils) and the displacement zone (for rocks) are 
estimated from the apparent crater, and their dimen- 
sions are 


(10-7) 


(10-5) 


+ z2 (B' + 2C 'zW-l/3)2 Q2 
z 


and 


Q2 
= Q2 + l..-1[A'2W2/3 + B'2z2 + C'2z4W-2/3]Q2 
Ra 
[R2 
CO 
a 


(10-8) 


Q2 + Q2 
+ 
Q2 
-+- 
Q2 
g 
r 
~ 
Ra 


Q2g + 
Q2 
+ 
Q2 + Q2 
r' 
~ 
Da 


Q2Rp 


Q2Dp 


Values of rand r' are given below. 
(1) The mean plastic or displacement zone radius 
and depth are calculated from the values presented 
in table 10-1 and equations 10-1, 10-2, and 10-7. 
(2) The uncertainties associated with the plastic 
or displacement zone are 


R' 
ERa 
D' 
E'Da 
a 
a 
R' 
ERt . 
D' 
E'Dt 
t, 
t 
R' 
ERR 
DR 
E'Dr 
R 
R' 
ERp 
D' 
E'D 
(10-9) 
p 
p 
p 
10-7 


~here Qr. Qr', and Q~ are the uncertainties in equa- 
tIOn 10-7, and QRa and QDa are defined in a (2) above. 
Assume that Q2 + Q2 + Q2 = Q2 + Q2 + Q2 :::: 005 
g 
r 
~ 
g 
r' 
~ 
.. 
10-3. 
Multiple, pattern burst. 


Close spacing of charges (less than 1.4 single-charge 
crater radii) enhance both radius and depth of crater- 
ing. Figure 10-6 shows the enhancement of single- 
charge crater dimensions obtained by reducing the 
spacing between charges in a row-charge at optimum 
DOB. The enhancement is proportional to lIvs where 
s is the charge spacing in ft. As the charge spacing 
decreases, the optimum DOB for each charge must 
be increased by the enhancement factor in figure 10-6 
in order to achieve these results. Since the enemy may 
not be able to achieve penetration to the optimum 
DOB, figure 10-6 is an upper bound. Nevertheless, 
it is recommended that the enhanced crater dimen- 
sions be taken as 


(l0-4b) 


W2/3 
t 
+ -g- (A' - C 'z2W-2/3)2 
Q~\, 


z 
16.4 < Wl/3 < 100 


Values of e, e', and ~ are defined below. 
(1) The mean rupture zone radius and depth are 
calculated from the values presented in table 10-1 
and equations 10-1, 10-2, and 10-5. 


In these equations Q[ is the uncertainty associated 
with the functional. form ofRa, and Qc is the uncer- 
tainty associated with the constants in ~he equations. 
Assume that Q[ :::: 0.36, Qi + Q~ :::: 0.01; assume that 
Qco = 0.07 for reasonably homogeneous sites and 
0.12 for layered sites or for sites where a water table 
intersects the crater. These uncertainties were ob- 
tained from equation 2-17 where n = ~ and R = 1.35 
for the former case (homogeneous) and R = 1.67 for 
the latter (layered). The calculation of Q2 parallels 
equation 10-4 where the notation Ra is r~placed by 
Da · 
b. Rupture zone. The rupture zone is also estimated 
from the apparent crater, and its dimensions are 


•• 


• 


• 


o 
I 
CD ' 


Table 10-1. 
Mean Value of Parameters (equations 10-1 & 10-2) 


Parameter 


Site 
a' 
B' 
y' 
a 
B 
y 
A' 
B' 
C' 
A 
B 
C 
, 
8 
8 ' 
r 
r' 


Dry Sailor 
-5 
-7 
-3 
-1 
1.D9xl02 
_3.43xlO- 3 
-1 
3.74xlO l 
Dry Soft Rock 
48.6 
4.7 
-7.43xlO 
22.3 
1.9 
-2.42xlO 
-3.43xlO 
7.75xlO 
9.50xlO 
1.4 
2.0 
1.5 
6.0 
3.0 


Dry Hard 
-6 
-6 
-3 
7.44xlO l 
-4.06xlO- 3 
9.50xlO- l 
3.16xlO l 
Rock 
39.0 
6.3 
-3.80xI0 
17.4 
2.5 
-1.2IxI0 
-5.94xI0 
1. 35 
1.4 
2.0 
3.0 
4.0 
4.0 


Wet Sailor 
-7 
_2.23xlO- 7 
-6.25xlO- 3 
9.80xlO l 
-4.69xlO- 3 
4.14xlO I 
Wet Soft Rock 
65.0 
4.2 
-5.57xlO 
24.6 
1.9 
1.68 
1.13 
1.4 
2.0 
1.5 
6.0 
3.0 


Wet Hard 
-7 
-7 
-3 
9.41xlO l 
3.14xlOl 
Rock 
45.9 
4.4 
-6.13xlO 
22.3 
1.7 
-1. 30xlO 
-5.31xlO 
' 
1. 25 
-4.69xlO- 3 
1.18 
1.4 
2.0 
3.0 
4.0 
4.0 


u.s. Army Corps of Engineers 
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Figure 10-6. 
Enhancement ofSingle-Charge Apparent Crater Dimensions in a Row Crater as a Function ofCharge Spacing at Optimum 
DOB in Soil (Rooke et aL, 1974) 
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Similar equations apply to the depth dimensions, i.e., 


TM 5-858-2 
where the enhancement factors are of the form 


'.. 
• 


• 


C10-13) 


ClO-12) 


(10-14) 


and the enhanced dimensions are 


where E" is the increase in true crater, dimensions 
and N is the number of successive shots. The 
parameters A and Bare defined in b" below. 
b. The mean increase in true crater dimensions can 
be obtained from figure 10-7 or from equations 10-12 
and 10-13, with the constants defined in table 10-2. 
c. The uncertainties of the true crater dimensions 
from the nail-driving attack are 


Q2,,::;: Q2 + CA2 + B2N2JRt)2 Q2 
+ N2B2 Q2 + Q2 
Rt 
r 
\Rt' 
CR 
R,,2 
N 
Rt 


. 
,'.' 
,,' t 


formed by preceding shots to enhance cratering. 
Bounds of the true crater data for tests in rock are 
shown in figure 10-7. These data are represented by 
equations of the form 


E" 
:= A + BN 


(l0-10a) 


C10-10b) 


:= ACslR )1/2 
a 


:= ACsiD )1/2 
a 


E 


Q~, 
Q2 
+ Q2 
+ 1 Q2 + 
1 Q2 
a 
r 
A 
:I" 
s 
:I" 
Ra 


Q~t 
Q~~ + 
Q2 


t 


Q~R 
Q~t + 
Q2e 


Q~p 
Q~a + 
Q2 
+ 
Q2 
(lO-lla) 
r 
t 


E' 


where Qr is the uncertainty of equations 10-9, and 
QA is the uncertainty of the constant in equation 
10-10. Assume Qr :::: 0.04 and calculate QA from 
figure 10-6 and equation 2-17 where n := 2 and L2 
and L1 are the upper and lower bounds on the data 
at any selected burst spacing. Similarly, the uncer- 
tainties of the crater depths are calculated from 


a. Use the mean value of the data range presented 
in figure 10-6 to determine the mean enhancement 
factor. This occurs when A := 1.2 in equation 10-10. 
b. The uncertainties associated with equations 
10-9 and 10-10 are 


Q~, 
Q2 
+ 
Q2 
+ 
1 Q2 
+ 
Q2 
a 
r 
A 
:r 
s 
Da 


Q~t 
Q~, + 
Q2 
a 
t 


Q~R 
Q~t + 
Q~, 


Q~p 
Q~, + 
Q~, 
+ 
Q2 
ClO-11b) 
a 
t 


10-4. 
Repeated attack (nail driving). 


a. Repeated attack, sometimes called "nail driving," 
is the successive firing of shots in the center of craters 


where Qr is the uncertain~y of the form of equation 
10-13, QCRand QeD arethe uncertainties of the nail- 
driving enhancement factors applied to the radius and 
diameter, respectively. Assume Qr :::: 0.2 and calculate 
QCR and QCD from equation 2-17 where n := 2 and 
L2 and L1 are the upper and lower bounds on the 
enhancement factors in figure 10,...7 for any selected 
number of shots. 


10-5. 
Crater volume. 
Crater volumes in homogeneous and in layered 
geologies are calculated in chapter 12, where the 
depth of burial by crater ejecta is defined. 
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Figure 10-7. 
Enhancement in TrucLCrater Dimensions {or Nail-Driving (Rooke et aI., 1974) 
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10-12 


Table 10-2. 
Parameters for Increase of True Crater Dimensions 
Due to Nail Driving 


Geology 
Dimension 
A 
B 


Rock 
Radius 
0.88 
0.12 


Depth 
0.45 
0.55 


Soil 
Radius 
u.8l 
0.19 


Depth 
0.99 
0.01 


U.S. 
Army Corps of Engineers 
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CHAPTER 11 
EJECTA/DEBRIS IMPACT 


(11-3) 


(11-6) 


(l1-4a) 


(11-4b) 


Range. ft 
RNl/3 
, 
a 
Constants 


Therefore, the maximum particle diameter is 


where 
Va 
= Apparent crater volume, ft3 
it follows that 


where the particle configuration is assumed to be 
spherical, that is, 


In turn, the maximum particle weight at any range 
(from unpublished work by Post) is determined from 
the total weight to be 


A ={ 1.08 (hard rock) 
n = j-0.0477 , W ~ 1 kt 


0.54 (other geologies) 
I 0 
, W < 1 kt 


b. The mean value for equations 11-1 through 11-5 
can be obtained from the values of the constants 
presented below: 


where Qf = 0.36 and Q 
= 0.05. (Q2\', 
is defined in 
Yp 
a 
ch. 12.) 


where 
R 
Rs 
B,m,n,a 


B 
0.1094 
m 
0.8 


a 
-1.644 
c. The uncertainty associated with equation 11-5 is 


Q2 
= Q2 + (m-1)2 Q2 
(mn)2 Q2 
am 
f 
9 
l'p + -9- 
w 


2 
~ 
)') 
+ ~ Q2 + 
m _ ~ - Q2 
9 
R 
3 
9 
V 


(11-1) 


(11-2) 


Apparent ejecta weight, lb 
Total ejecta weight, lb 
Yield, kt 
Constants 
depending 
on 
geologic 
properties 


11-1. Gross phenomenology. 


a. Material is thrown out of the crater region as a 
part of the unloading process following the initial 
strong shock driven into the ground by the explosion. 
Ejected particles produced by a near-surface burst ex- 
tend over a range from the crater lip to remote points. 
Discrete particles are thus part of the crater forma- 
tion, 
extending through 
the 
continuous 
ejecta 
coverage, on out into the discontinuous region. The 
resulting projectiles of debris create the ejecta/debris 
impact hazard and may pose a substantial threat to 
exposed structures. (Projectiles may also result from 
fragmented structures and trees between the point 
of interest and ground zero, GZ. These are of second- 
ary importance and are not treated here.) Ejecta im- 
pact encompasses the problems of estimating the 
initial mass of ejecta and the' size of its consequent 
particles, the number of impacting particles, their in- 
coming velocity and angle, and their distribution. 
Designers must also consider the damage potential 
of ejecta impact in its timing relative to other nuclear 
weapon effects. 
b. The basic variables affecting the number, size, 
velocity, and angle of approach of the impacting 
ejecta/debris projectiles at target points are: 
- Range to GZ 
- Weapon yield 
Weapon design 
Height of burst 
Geology 
- Target elevation 


11-2. Ejected mass and particle size. 


a. The effects of ejecta impact are determined from 
calculations of total ejected mass, number of impact 
particles, mean maximum particle size and number, 
and mean impact velocities and angles. Using the 
data of Post (1974) and Crawford et al. (1974) state 
that the total ejected mass can be calculated from the 
apparent crater mass with the equation 


where 
Ma 
Mt 
W 
A,n 


Since 


• 


• 


11-1 


l 


0.00147 for unconsolidated alluvial soil 


s = 
0.5 for hard rock, cohesive soils, or layered 
geologies 
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11-3. Number of impact particles. 


a. It is useful to express the number of particles equal 
to or greater than a specified size that will occur at 
a specified range from GZ. Such an equation is of the 
form 
• 


• 


where 


v·L 
CD 
R 


Impact velocity, ft/sec 
Drag coefficient of particles 
- 
Range from GZ, ft 
. 


Preattack density of ejecta, Ib/ft3 
90a . 
~sin 2e, ft 
CD 
ao 
Arbitrary particle diameter, ft 
e 
Ejection elevation angle, deg 


~ fJJ 
Constants 
b. At an arbitrary range, the incoming particle size 
may be selected in the range a 
.:s:;; a 
where a 
is 
b 
. 
df 
0 
m 
m 
o tame 
rom equation 11-5. The ejection elevation 
angle is assumed to be distributed between 15 deg and 
75 deg, but strongly concentrated near 45 deg, as 
shown in figure 11-1. The drag coefficient is assumed 
to lie between 0.3 and 1.2, with a mean value of 0.6. 
The mean impact velocity is computed from equation 
11-9 where A = 5.673 and $= 23.08. 


c. The uncertainty in the impact velocity is given 
by 


Q2 = 
v·I 


(11-7) 
a8- 3 _ a8- 3 
m 
0 


8 
am 


K 


where 
I 
Number of particles per unit area (ft2) 
greater in size than ao(ft) (an arbitrary 
particle size) but less than am (ft) 
s 
= Constant depending' on geology 
6y 
nyp 
. 


y 
Postattack unit weight of ejecta, Ib/ft3 
Equation 11-7 was derived by Post (again, in un- 
published work) based on a review of empirical data; 
it implicitly contains the range parameter from the 
variables De and am where De is the continuous ejec- 
ta depth, as defined in chapter 12. 
b. The mean number ofimpacting particles per unit 
area is obtained from equation 11-7. For ao ~ 4 in., 


KsDe 
I=-- 
(s-3) 


c. An estimate of the uncertainty associated with 
equation 11-7 is given by 


Ql ~ Q; + Qt, + }~::;55-_~:512 Q;m 
,(11-8) 


where Qr;:::; 0.2, Qam is obtained from equation 11-6, 
.and QD is presented in equation 12-2. 
e 


11-.4. Particle impact velocities and angles. 


a. Computational and experimental evidence support 
the adoption of a model of ejecta trajectories based on 
the assumption that particles are ejected at a random 
elevation angle eand at a velocity dependent on size, 
and that the trajectory is determined by gravity and 
aerodynamic drag forces. Analysis of the trajectories 
obtained from calculations have led to the simplified 
description of impact velocity reported by Crawford 
et a1. (1974): 


~ 


R 
~1/2 
v· = -- 
exp (-RIR) 
I 
sin 2e 
0 


~fB(~:"~ 


11-2 


, 
R>RJ2 


(11-10) 


where Qr = 0.36 and QeD = Qe = 0.3. 


11-5. Impact damage 1P0tentiai. 


a. Impact damage potential can be measured in 
terms of particle kinetic energy (momentum), ex- 
pected target penetration, or expected crater volume. 
For hardened facility applications, use the damage 
measures-penetration or crater volume. Penetration 
has been examined as a continuing part of armor/ 
armor-piercing projectile development and, more 
recently, as part of the safety program for nuclear 
power reactor development. Available data for pro- 
jectiles striking concrete targets are utilized to ex- 
trapolate the effect ofrack fragments striking similar 
targets. 
b. The penetration depth, d, is adapted from ASCE 
(1976): 


',' 
t 


o.04 FD(KNYp)1I2 apVio.9 
, 
dlao.:s:;; 2 
d = 
2 X 10-6 FDKNYpapVi 1.8 
, dlao > 2 (11-11) • 
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Figure 11-1. Ejection Angle Distribution for High-Yield Surface Bursts (Seebaugh, 1975) 


1 
2 
1 
2 
O'SQ2 
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Yp 
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c. The uncertainty associated with equation 11-11 
is 


(d-a)2[ 
__0 
Q2 
+ Q2 + Q2 + Q2 
d 
FD 
KN 
~ 


+ 3.24 Q2VJdlao >2 
(11-13) 


1 
where Qf::::: 0.5, QFD ::::: 0.25, QK 
=:: ~;., QN ::::: 0.15, 
Q 
=:: 0.05,andQv·isobtainedfromequati~mll-10.The 
yP 
I 
quantity Qfc is the uncertainty of the unconfined com- 
pressive strength of concrete. 
(11-12) 


f'c 


Vcrater 
=:: 6.5 d3 


where 
FD 
Particle deformability factor 


{ 
0.7 for hard rock 


::::: 0.25 for soft rock or alluvium 


K 
=:: Concrete penetrability index =:: ~ 


Ultimate compressive strength of con- 
crete, psi 
N 
=:: Particle nose shape factor::::: 1.0 
ao • =:: Arbitrary particle diameter ~ am' in. 
Vi 
=:: from equation 11-9 based on particle 
size ao' ftlsec 
Yp 
=:: Preattack particle density, lb/in.3 
The estimate of penetration crater volume is based 
on an assumed conical crater depth d and radius 2.5d 
• 


11-3 


• 


• 


• 
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CHAPTER 12 
EJECTA/DEBRIS DEPTH 


12-2. Debris depth. 


a. The depth of crater ejecta debris from a near- 
surface burst is related to the apparent crater volume. 
(See chapter 10 for calculation of apparent depth and 
radius.) In turn, the crater volume is implicitly depen- 
dent on weapon yeild, height of burst, weapon type, 
and geology (Post, 1974). The scaled debris depth as 
a function of scaled range can be represented by an 
equation of the form 


12-1. Gross phenomenology. 


a. To prevent burial of a facility and interruption of 
its function or to plan active removal of debris, the 
designer must know the expected depth of debris, in 
addition to the expected particle size distribution 
delineated in chapter 11. Deposition of crater ejecta 
as surface debris becomes significant when the height 
of burst is less than 35 ftlkt1. 3 (10.7 mlkt1l3). The 
depth of debris in the region of continuous coverage 
of the ground depends on range from ground zero (GZ) 
and varies (usually unpredictably) with azimuth, 
probably the result of fracture patterns and in- 
homogeneities in the medium. Outside the continuous 
region (greater than five crater radii), the debris con- 
sists of separated discrete particles unlikely to bury 
the facility. The amount of ejecta threatening burial 
of a facility depends on the crater volume and is ad- 
dressed in this chapter for both homogeneous and 
layered geologies. 
b. The basic variables effecting the ejecta thickness 
at a point on the ground are: 
- Range to GZ 
- Weapon yield 
- Weapon design 
- Height of burst or depth of burst 
- Geology 
Ejecta depth is enhanced by large yields, reduced 
ranges, weapon penetration, and certain geological 
features described below. 


b. The 
mean 
average 
ejecta debris 
depth 
is 


~alculatedfrom equation 12-1 where the constants 
have the values: 


{1 


-0.423 
a 
-2.3 
B 
0.017 
A 
0.75 


Q2 
= Q2+ _1_ [D - AaRea R/Va 1l3 
De 
[ 
9D2 
e 
e 


- B {1 Re f.l RNa113rQ~a 


where Q[ is the uncertainty of the functional form of 
equation 12-1, Qv is the uncertainty of the apparent 
a 
crater volume, etc. Calculate Q[ from equation 2-20 
where L1 and L2 are the lower and upper bounds 
from figure 12-1, multiplied by V~/3. The uncertainty 
of Va is obtained from paragraphs a (2) and a (3) 
below. 


Equation 12-1 is plotted in figure 12-1 as the median 
value with uncertainty bounds. These data are ap- 
plicable for depths of burst less than 17 ftlkt1l3 


(5 mlkt1l3) and for most of the layered geologies en- 
countered within the continental United States. 
Equation 12-1 should be multiplied by the quantities 
in table 12-1 for the other geologies shown therein. 
The mean apparent crater volumes are calculated 
from the procedures presented in paragraph 12-3 
below. 
c. The uncertainty associated with equation 12-1 is 


12-3. Crater volumes. 


a. Homogeneous geologies. 
Normalized apparent 
crater volumes are presented for both "high-yield" 
and "low-yield" nuclear devices. These terms from 
Cooper (1971) refer to "hot" devices having a high 
yield-to-mass ratio and to "cold" devices having a low 
yield-to-mass ratio. The relationship of volume to dif- 
ferent ranges of HOBIDOB are also considered for 
both high-yield and low-yield devices. Typical crater- 
ing efficiencies are shown in table 12-2 for half-buried 
HE charges. 


(12-1) 


where 
De 
Average ejecta thickness, ft 
Va 
Apparent crater volume, ft3 
R 
Range from ground zero, ft 
A, B, a, and {1 are constants 


~ = A exp (a ----.lL) + B exp ({1 ----.1L) 
V 1/3 
VI/3 
V 113 
a 
a 
a 


," 


• 
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Figure 12-1. 
Ejecta Depth Prediction Curves for Nuclear Events, 1 ~RlVIl3 ~ 12 (Post, 1974) (l of2) 
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Figure 12-1. 
Ejecta Prediction Curues for Nuclear Euents, U";; R/V 1/3.,;; 40 (Post, 1974) (2 of 2) 
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Table 12-1. 
Geology Correctir:n Factors 


~1edium 
Correction Factor 


Hard Rock 
2.0 


Porous Rock 
0.75 


U.S. Army Corps of Engineers 


Table 12-2. 
Normal Cratering Efficiencies E for HalfBuried HE' 
Charges (Post, 1974) 
. 
o . 


Idealized Geology 
E 
ft 3/kt 
a 


Wet Soil 
4 x 106 


Wet Soft Rock 
2 x 106 


Dry Soil 
1.4 x 106 


Dry Soft Rock 
1.0 
6 
x10 


Hard Rock 
0.6 x 106 


• 


• 


• 


(1) The mean apparent nuclear crater volumes for 
100 < W < 10 kt can be described by the following: 


1M 5-~~8".2 


for the DOB conditions. The values of the constants 
A and yare defined as: 


Values of Constants 
• 


"High-yield" weapons: 


VaN0 = 0.04 e-55·1Hs 


VaNo = 0.05 e- 166 Hs 
VaN0 = 4.07 (0.0033 + Ds)0.77 


"Low-yield" weapons: 


Hs~ 0.002 
o ~ Hs~ 0.002 
o ~ Ds~ 0.02 
(12-3) 


Assumed 
Condition 


"High-yield" weapons 
Hs 
~ 0.002 
o ~ Hs < 0.002 


O~ Ds < 0.02 


A 


0.0033 


y 


-55.1 


-166.0 


0.77 


VaNo ~ 4.39 (Ds + 0.0216)°·88 
VaN0 = 4.07 (0.0033 + Ds)O.77 


O~ Ds<0.055 


0.055 ~ Ds< 0.2 
(12-4) 


"Low-yield" weapons 
o ~ D s < 0.055 


0.055 ~ Ds < 0.2 


0.0216 


0.0033 


0.88 


0.77 


The quantity .Q[ is the uncertainty of equations 12-3 
and 12-4, QHOB is the uncertainty of the HOB, etc. 
Assume Q[ :::: 0.2 and let 


(12-9) 


where QE = 0.3 based on the use of equation 2-20 
for an 
u~certainty band of ± 50% on the data 
presented in table 12-2. 
(3) An 
approximation 
of the 
uncertainties 
associated with equations ·12-5 and 12-6 is 


HOBN 1/3 
o 
DOBN 1/3 
o 


WTEo. 
. 


Cratering efficiency in ft3lkt 


Yield in kt 


(12-5) 


"High-yield" weapons: 


[ 
.. ]0.77 
VaN 0 = 4.07 Wa 0.0033 + DsWaJ3 


in which HOB and DOB are both expressed in ft. For 
yields greater than 10 kt, the DOB cases of equations 
12-3 and 12-4 should be replaced by the. following: 
• 


"Low-yield" weapons: 
[ 
YD WaJ3 
. ]2 
. 
s 
Q2 
+ 
(A + D W1/3) 
DOB 
s 


. [1 
YDsW 
aJ3 
]2 
Q2 
+ 
- 3(A + D Wl/3) 
Yo 
s 
(12-10) 


(12-6) 


where a = -0.3 DsWO.02 
(2) The uncertainties associated with equations 
12-3 and 12-4 are specified as 


where Q[ is the uncertainty of equations 12-5 and 
12-6 and Qy is the uncertainty ofVo' The constants 
for the variot.s conditions are: 


y 
0.77 


Value of Constants 


A 
0.0033 


Assumed Condition 


"High-yield" bursts: 
"Low-yield" bursts: 
o < Ds 
~ 0.055 
0.0216 
0.88 
0.055 < Ds < 0.2 
0.0033 
0.77 


Assume Q[:::: 0.2 and use equation 12-9 to compute 
Qy 
. 
o 


(12-7) 


(12-8) 


for the HOB conditions and 
• 
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Equation 12-11 can be written for a layered site of 
n different materials via the recursive equations: 


• 


• 


(V1,2,...,n-1,n-V2,...,n_1,n)/(V1-V2,., .,n-1,n) 


= 1 - exp [-5.4d1(V1,2,...n_1,n)1I3] 
(12-12n) 


Equations 12-12 must be solved sequentially, begin- 
ning with the first equation and working upward 
through n materials. 
(1) For convenience, the mean for equations 
12-12 coupled with the equations and figures dis- 
cussed in previous subparagraphs are presented in 
figures 12-2 through 12-6 for a few interesting 
geologies. In figures 12-2 through 12-4, cratering ef- 
ficiencies can be determined for three different 
geological conditions. Based on these cratering effi- 
ciencies, 
the 
mean ejecta-debris depth 
can 
be 
calculated from figure 12-5 for airblast pressures of 
600 and 1000 psi and in figure 12-6 for a larger selec- 
tion of overpressures. (Refer to figure 7-10 to deter- 
mine overpressure/range relationship.) 
(2) The uncertainty of ejecta depth for layered 
sites should be calculated from equation 12-2. Equa- 
tion 12-2 in turn requires the uncertainty Qv . 
Depending on the HOB or DaB and the yield and type 
of weapon involved, Qv should be determined from 
one of equations 12-7 tlfrough 12-10 (omitting equa- 
tion 12-9, which is used to compute the uncertainty 
of Yo). Equations 12-7, 12-8, and 12-10 require a 
value of VO' It is recommended that V0 be calculated 
from the applicable condition in equations 12-3 
through 12-6 by iterating the value of V0 until the 
ratio VaN0 on the left-hand side of the equation is 
compatible with the right-hand side that involves 
scaled HOB's and DaB's. 


(12-11) 


(12-12b) 


(12-12a) 


Layer thickness in ft 
Crater volume (at the HOB or DaB 
of interest) when d = 00 


Crater volume (at the HOB or DaB 
of interest) when d = 0 


where 
d 
Vu= 


b. Layered geologies. Apparent crater volumes can 
be determined for layered geologies from the equation 


(vn-1,n-Vn)/(vn-1-Vn) 


(Vn-2,n-1,n -Vn-1,n)/(vn-2-Vn-1,n) 


• 
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Figure 12-2. 
High-Yield Nuclear Surface-Burst Cratering Efficiencies in Layered Geologies, Soil over Soft Rock (Shale or Sandstone) 
with a Water Table (Post, 1974) 
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Figure 12-3. 
High-Yield Nuclear Surface-Burst Cratering Efficiencies in Layered Geologies, Soil over Hard Rock (Limestone) with a 
Water Table (Post, 1974) 
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Figure 12--4. 
High-Yezld .vuclear Surface-Burst Cratering Efficienc!es In La)ered GeologIes, AllUL'ial Soil Ol'er Hard Rock Baser 
u'ith a Water Table 'Post, 1974) 
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Figure 12-5. 
Scaled Ejecta Depth Vs. Cratering Efficiency at 600- and lOOO-psi Levels (Post, 1974) 
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Figure 12-6. 
Scaled Median Ejecta Depth vs. Cratering Efficiency at Varying psi Levels Post, 1974) 
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CHAPTER 13 
DUST 


(13-1) 


(13-2) 


(13-3) 


(13-5) 


(13-4) 


Q2 = Q2 + Q2 
+ n 2Q2 
R 
r 
RI 
W 


where Qr is the uncertainty associated with the func- 
tional form ofR, QRI is the uncertainty parameter 
associated with RI , etc. Assume 
Qr 
:::: 0.2 and 
QRI = 0.2. 
(4) The cloud height is given by 


is 


where 
R 
Cloud radius, m 
RI 
Cloud radius for a 1 Mt burst, m 
W 
Weapon yield, Mt 
n 
Constant 
(2) The mean radius is obtained from equation 
13-1 where RI = 3000 m and n = 0.4. 
(3) The uncertainty associated with the radius R 


(1) The cloud radius is given by 


Assume Qr :::: 0.2, QHI = 0.2. 
b. Initial dust concentration. Assume that the top 


~ meter of dust cover is uniformly dispersed into the 
cloud. Accordingly, the initial dust concentration is 


is 


Thickness of dust layer on the earth's 
surface that is dispersed into the cloud, 
m 


Q 
Density of the in-situ earth layer, g/m3 


!Po 
Initial dust concentration, g/m3 
H 
Dust cloud height (eq. 13-3) 


(1) The mean initial 
dust concentration is' 
obtained from equation 13-5 where 
~ = 0.5 m or 
whatever depth applies to a particular site, and 
e= 1.0 x 106 g/m3. 


~Q 
!P =- 
o 
H 
where 


~ 


where 
H 
Cloud height, m 
HI 
Cloud height for a 1 Mt burst, m 
m 
Constant 


(5) The mean height is obtained from equation 
13-3 where HI = 100 meters and m = 0.2. 
(6) The uncertainty associated with the height H 


13-1. Gross phenomenology. 


a. Dust is lofted into the air by the deposition of the 
initial radiation on the ground surface, by the crater- 
ing process, by the outward sweep of the airblast 
across the ground surface, and by the inward sweep 
ofthe afterwinds produced by the rising fireball. The 
lofted dust is quickly partitioned into a low-lying or 
base surge cloud and a high-altitude or ejecta dust 
cloud. 
b. Only recently have dust clouds been accorded at- 
tention by the nuclear weapon effects community. The 
material presented in this chapter is a synthesis of 
the literature on the subject (see, for example, Young, 
1965; Tewes, 1970; Crawford, 1970; Allison, 1971; 
Swatosh, 1971; Powers, 1975). 
c. The basic variables 
affecting the 
potential 
destructiveness of the free-field dust environment 
received at a ground surface location of interest here 
are: 


- Weapon yield 
- Burst elevation 
- Dustiness of the ground surface 
- Dustiness of the crater ejecta 
- Prevailing winds 
- Height of temperature inversion (mixing depth) 
Dust effects increase with larger weapon yields, lower 
heights of burst, dustier ground surfaces, calmer 
winds, or lower inversion bases. 
d. Personnel will rapidly become incapacitated if 
they breathe the heavily dust-laden air. Appropriate 
design measures must be taken either to remove the 
dust from the air or to allow for functioning in a 
sealed facility. Knowledge of two aspects of the dust 
cloud are required 
to implement dust removal 
measures or button-up procedures: 


. Particle-size distribution vs. time 
Load vs. time 
For calculations the ground-surface dust and the 
ejecta dust are treated separately as though they 
formed two noninteracting clouds. 


13-2. Base-surge dust cloud. 


a. Cloud model. The model of the base-surge dust 
cloud is developed on the assumption that it is instan- 
taneously formed and initially centered at the weapon 
aim-point. The model is configured as a right circular 
cylinder of constant radius R and constant height H 
and 
in 
which 
the 
initial 
dust 
concentration 
throughout the cloud, !Po' undergoes no dilution or 
diffusion of the cloud with the surrounding air but 
decreases in concentration due to vertical settlement 
ofthe dust particles. The cloud will move horizontally 
at the prevailing wind velocity. 
• 


• 


• 
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(2) The uncertainty associated with the initial 
dust concentration is 


(13-6) 


where Q;:::: 0.4, Q~ :::: 0.4, and Qe = 0.2. 
c. Vertical settlement speed The vertical settlement 
speed of a dust particle is 
t = 
a 


o 


H1WO.2 


VO(1:00) 


a> 1000 !-1m 


100 ~ a ~ 1000 )..1m • 


• 


(13-10) 


(13-9) 
a < 10 jJm 


10 ~ a: ~ 100 !-1m 


co 


H1wo.2 


( 
) 


3 
V' 
a 
o 100 


where QH and Qv are calculated from equations 13-4 
and 13-8. Assume Qr :::: 0.36. 
f Load us. time. The dust load that is entrained in- 
to the facility is ofcritical importance to thedesigners. 
The load is determined by integrating the. product of 
the dust concentration and the air-intake-volume rate 
over the elapsed time of interest: 


only particles smaller than diameter a are in the 
cloud. 
(1) The mean value ta is obtained from equation 
13-9 where the mean value of H and V are obtained 
from a (5) and c (1) above, respectively. The mean 
value of equation 13-9 is plotted in figure 13-2. 
(2) The uncertainty associated with ta is 


(13-8) 


a < 10 !-1m 


(13-7) 


a> 1000 !-1m 


10 ~ a ~ 100 !-1m 


100 ~ a ~ 1000 !-1m 


co 


VI(~)3.0 


0 100 


o 


V= 


where 
a 
Vo' V~ 


Particle diameter, )..1m 
Settling speeds for a 1000 !-1m and 
100 !-1m particle, respectively, mls 
The meaning of V = co for a > 1000 !-1m is that parti- 
cle sizes larger than 1000 !-1m are not lofted. 
(1) The mean vertical settlement speed of a par- 
ticle of diameter a is obtained from equation 13-7 
where V0 = 55 mls and V~ = 5.5 mls. Mean settle- 
ment speed is plotted in figure 13-1. 


(2) The uncertainty associated with the settle- 
ment speed is 


Q~ = 0.13 


For the assumption of a constant air-intake-volume 
rate and of a particle diameter distribution of a-4, the 


L __ f.a2fat 
1fI(a,t)~(t)dtda 


al 
0 


d. Wind 
transport. 
Assume that the cloud is 
transported as a rigid body by the prevailing wind. 
Assume that the cloud is initially centered above 
ground zero and thereafter moves with the prevail- 
ing wind. Use local wind statistics (contact the 
National Weather Records Center, operated by the 
Environmental Services Administration at Ashville, 
NC, for meteorological data) to determine transport 
speed and direction of the cloud. Track all clouds that 
originate in the general vicinity of the point of in- 
terest to develop probabilities of their passing over 
the point of interest. Appropriately account for the 
"contribution" of each cloud. 
e. Particle size us. time. With the passage of time, 


progressively smaller and smaller particles settle out 
ofthe Cloud. In effect, the cloud progressively exhibits 
a finer distribution of dust. At time 


where 


L 
t 


~ 


1fI 


(13-11) 


Dust load, g 
Time, s 
Air-intake-volume rate, m3/s 
Dust concentration at an arbitrary 
time t, g/m3 
Size of smallest and largest particles 
in cloud, m 


co, 


• 
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Figure 13-1, Mean Vertical Settlement Speed ,ofDust Particles. 
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Figure 13-2. Mean Settle-Out Time of Dust Particles-Base Surge 
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load attributable to particles of diameters al 
~ a 


~ az is 


10H 
-- 
s 
Vo 
' 


(13-14) 


(13-15) 


v = v (2Ra)2 
a R' 


R' = (8R3 + 12v RZt)1I3 
a 
a a 


where 
Ra 
Apparent crater radius, m 
R' 
Cloud radius at time t, m 
t 
Time, s 
va 
Cloud initial radial surge velocity, mls 
(1) The mean cloud radius at time t is obtained 
from equation 13-15 where the mean value of Ra is 
defined in chapter 10 and Va = 200 mls. The cloud 
will stabilize when the mean surge velocity is very 
small compared to the mean prevailing wind velocity. 
(2) The uncertainty associated with the cloud 
radius at time t is 


where Ra = apparent crater radius, m; R' = cloud 
radius at time t, m; v = cloud radial surge velocity 
at R I, mls; va = cloud initial radial surge velocity, 
mls. The cloud is also assumed to be initially centered 
at the weapon aim-point and has a uniform initial 
dust concentration throughout the cloud, 
4J~, which 
undergoes no dilution or diffusion of the cloud with 
the surrounding air, except for the horizontal expan- 
sion stated in equation 13-14. The cloud maintains 
a uniform, horizontal dust concentration at all times 
but the dust concentration decreases due to vertical 
settlement. of the dust particles. The cloud moves 
horizontally at the prevailing wind velocity. 
b. Radius. The cloud radius is 
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apparent crater 
radius 
and thereafter 
expands 
radially according to the equation 


(13-12) 


V'o 


1000 H , s 


azH ,m 
Vot 


10 a (~)l/3 ,m 
o V't 
o 


Smallest particle that settles, m 


H 
- 
s 
Vo ' 


= 


where 


a* 


a** 


ao 


tl 


tz 


t3 


4J 


• 


• 


(1) The mean load is obtained from equation 
13-12 by substituting the appropriate mean values 
for variables where al = 1 !Jm az = 1000 !Jm. 


(2) The uncertainty associated with the load L is 


13-3. Ejecta dust cloud. 


a. Cloud modeL The ejecta dust cloud model assumes 
that the cloud instantaneously reaches a height equal 
to the base of the upper atmosphere temperature in- 
version (H') and remains at that height, and that it 
instantaneously reaches a radius equal to twice the 


where Qf is the uncertainty associated with the func- 
tional form of R', Qv 
is the uncertainty associated 
with va' etc. Assume tlf ::::; 0.36 and Qv '= 0.36. Refer 
. 
a 
to chapter 10 for calculatIOn of QR . 
c. Height. The cloud height is eq~al to the height 
of the base of the temperature inversion. The mean 
height (based on a one-year interval) and the uncer- 
. tainty parameter of the inversion base can be deter- 
mined from weather statistics. 
d. Initial dust concentration. Assume that a frac- 
tion ofthe earth material that comprised the apparent 
crater volume is uniformlYdispersed into the initial 
ejecta dust cloud. Accordingly, the initial dust con- 
centation is 


(13-16) 


(13-17) 


(13-13) 


0.4, and QH is 


QZ = QZ + QZ + QZ + QZ + QZ 
L 
f 
Q 
~ 
H 
~ 


where 
Q f 
::::; 0.4, 
QQ = 
0.2, 
Q~ 


calculated from equation 13-4. 


,<I 


• 


13-5' 


f1000 


HOB 
H' 


where 
fHL 


". 


• 


1000 H'N~ , s 


10 H'No' s 


H'No' s 


fHL f1000 Vay 


h. Load us. time. The form of equations 13-11 and 
13-12 remain unchanged; however, utilize the follow- 
ing values of the parameters: 
Fraction of the ejecta dust which is 
partitioned into the dust cloud 
0.5 - HOB/500 , -250 < HOB < 250 
m1Mt1l3 
Fraction of the crater ejecta compris- 
ing 1000 !Jm diameter particles and 
smaller 
Scaled height of burst, m1Mt1l3 
Height of the inversion layer above 
the ground surface, m 
Va 
Apparent crater volume, m3 


y 
In situ earth density, g/m3 
(1) The 
mean 
initial 
dust-concentration 
is 
obtained from equation 13-17 where the mean values 
of Va and Ra are defined in chapter 10. Mean values 
of f1000 and H' are site specific. 
(2) The uncertainty associated with the initial 
dust-concentration is 
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Q2 = Q2 
+ Q2 
+ Q2 
+ Q2 
+ Q2 + Q2 
+ Q2 


tp~ 
ftp' 
[HL 
[1000 
Va 
y 
Ra 
H' 


(13-18) 
12 V 113 
. 
a 
! 


and in which Q[tp' ::::: 0.36, Qy = 0.17, and Q[[::::: 0.36. 
e. 
Vertical settlement speed. Paragraph 13-2c (for 
base-surge dust) is applicable. 
f 
Wind transport. Paragraph 13-2d is applicable. 
However, remember that the ejecta dust cloud radius 
is not constant but is increasing with time as defined 
in the model a 'above. 
g. Particle size us. time. Paragraph 13-2e is ap- 
plicable when H is replaced by H . 
• 


(13-20) 


In (a2/al) 


In (ae/al) 
f1000 


Q2 = 
Q2 + Q2 
+ Q2 + Q2 + lQ2 
L 
[ 
HL 
H 
y 
9 
va 


where Q[::::: QHL = 0.36, QH' ::::: 0, Qy= 0.05, and Qv 
is obtained from equations 12-9 or 12-10. 
a 


(1) The mean value of the ejecta dust cloud is ob- 
tained from equation 13-12 when ae = largest lofted 
ejecta particle ::::: 18 m, a1 = 1 !Jm, 
~2 = 1000 !Jm. 


(2) The uncertainty associated with the ejecta 
dust cloud is 


(13-19) 


where 


Q2 
+ Q2 + ( 
HOB/500 
)2 Q2 
[HL 
[[ 
0.5 - HOB/500 
HOB 


• 


• 
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CHAPTER 14 
FIRESTORM 


(14-1) 


(14-2) 


(14-3) 


a 


ec 


Ignition threshold radiant exposure, 
cal/cm2 


Volumetric heat capacity, cal/oC/ 
cm3 


Thickness of the material, cm 


Absorptivity of material 
= ~ atmax 
L 


a 
Thermal diffusibility ofthe material, 
cm2/s 


T 


L 


He 
Critical irradiance, cal/cm2/s 


tmax 
Time 
to 
maximum 
temperature 


;::: 0.032 WO.5, s/kt1l2 


W 
Weapon yield, kt 


K1, K2, K3, and n are constants 


ignitions will occur; if the conditions specified in it 
above are present, the fire will become a firestorm. 
The value of Qt in equation 14-1 is obtained from 
equation 6-16 in chapter 6 of this manual for a time 
of tmax as defined below. The value of Qo is also deter- 
mined below. 


A near-surface burst of > 10 kt will produce 10 
cal/cm2 (4.184 x 105 3/m2) radiation intensity at suf- 
ficient distances to ignite 300 acres for virtually all 
dry-season preattack conditions. 


b. Spontaneous ignition. Spontaneous ignition of 
kindling occurs when the incident energy density ex- 
ceeds specific values. Thus when 


where 
Qo 


(1) Ignition thresholds for dry kindling materials 
are expressed with equations of the form (Kerr, 1971) 


For moist fuels, equation 14-2 is modified by the fac- 
tor 1 + h-ver so that 
, 


where h = relative humidity in percent. 


14-3. Fuel density. 


a. Conditions. Prerequisites for a firestorm are: 
- Gross area aflame > 300 acres 
- Net fueled area -;- gross area> 25 percent 
- Firebreaks < 50 ft wide 
- Fuel density (based on gross area) > 175 tons/acre 
- Sufficient kindling/tinder available so that 10 
cal/cm2 of incident nuclear radiation will cause 
ignition 


14-1. Gross phenomenology. 


a. A firestorm, which can develop from the uniting 
ofthe many small fires started by the thermal radia- 
tion pulse, represents the most violent type of mass 
fire. The thermal convection column and strong 
ground-level winds blowing inward restrict growth of 
the fire, but intensify the burning. Firewhirls or fire 
cyclones may develop during the intensified burning 
phase. Full development of a firestorm can occur in 
an hour; burnout may not occur for several hours. 
b. A firestorm results in oxygen depletion, eleva- 
tion of ground and air temperatures (several hundred 
K), reduced or enhanced atmospheric pressures, and 
massive air pollution. Appropriate design measures 
must be taken, either to alter the site fuel-density 
characteristics sufficiently to reduce the probability 
of a firestorm to a safe level, or to provide the facility 
with heat shields and semi-closed-cycle ventilation, 
life support, and combustion air systems for the dura- 
tion ofthe firestorm. The designer must know the den- 
sity of the site fuel, the types of fuel, and the condi- 
tions that will ignite the fire in the thermal radia- 
tion environment (chap. 6). 
c. After a nuclear burst, the likelihood that a 
firestorm will develop is a function of: 
-' Weapon yield 
Burst elevation 
- Slant range 
- Weather conditions 
- Fuel density and kindling type 
The first three factors primarily affect ignition 
through their influence on incident radiant energy 
(see paragraph 6-2, Radiant Exposure). The other fac- 
tors account for the availability and condition offuel. 


14-2. Weather conditions. 


Firestorms can occur in the dry season and if wind 
velocities are less than 8 ft/s (2.4 mls), a calm wind. 
Climatological data of a facility site are required to 
estimate need for protection and level of protection 
against firestorms. 


• 


• 


• 
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(3) The uncertainty associated with equation 
14-2 can be reasonably expressed by 


.. 


• 


• 
(14-4) 


Vatmax 
in which r = 
- K 2 
. 
L 


and where Qr is the uncertainty of the functional 
form of equation 14-2 or 14-3, QK1 is the uncertain- 
ty of the constant K 1. For moist Kindling, and to a 
first approximation, the uncertainty associated with 
equation 14-3 is 


Q2 
= Q2 
+ eh2Q~L (Q~o + Q2) 
(14-5) 
QOm 
Qo 
Q2 
2 
h 
Om 
where Qis taken from equation 14-4 and Qh is the 


uncertal~ty of the value of the humidity, 
(4) Using existing experimental data, it is not 
practical to evaluate the 
various uncertainties 
associated with the various parameters in equations 
14-2 and 14-3. Instead, it is recommended that the 
calculated values for Qo and QOm be increased by the 
factor 1.15 and that Qr :::: 0.2 be used based on the 
data presented in table 14-2, in which the mean and 
standard deviation (exclusive of the wheat-straw data) 
were used to estimate the bias and scatter between 
calculated and experimental data. Assume Qc :::: 0.2. 
o 


~ 


O' 


K 3 = 
0, 
2,6,_ 


~ 


1, 
K2 = 
0, 
0, 


T <0.25 
0.25 < 
T <0.52 


T ~0.52 


(2) Mean values for parameters related to wild- 
land fuels are presented in table 14-1 for use in equa- 
tions 14-2 and 14-3. In addition, the following restric- 
tions prevail: 


~ 


-750 
-2250: 
-1522, 


1, 
1, 
0.4, 


The value of He is generally not a fixed value for all 
fuels but is dependent on the opacity and geometry. 
However, a value of 0.4 callcm2/s (1.7 x 104 J/m2/s) 
can be used in lieu of specific values measured in the 
laboratory. 


n 


. 


• 


l.4-2 


• 


Table 14-1. 
Physical and Thermal Properties of Wild Land Fuels (Kerr, 1971) 
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'When two values are given in this column and the one following it, the lower value of the 
absorptlvity 
(and hence the larger value of 
pcL/a) 
are for newly fallen leaves. 
All 
other values are for weathered materials. 
Few of these absorptivities have actually' been 
measured. 
Most are estimated from the appearance of the material. 


• 


..' 


• 


Fuel Materials 


Grasses and Other Fine Fuels 


Horsehair lichen 
(A~ectoria jubata) 


Cheatgrass seedpod (Bromue 
tecto1'UJ7l) 


Wiregrass (Aristrida stricta) 


Desert stipa (Stipa spp.) 


Cheatgrass leaves (Bromus tectorum) 


Harding grass 
(Pha~aris tuberoea 


var stenoptera) 


Sedge (Carex geyeri) 


Wheat straw (Triticum spp.) 


Cheat grass stalks (Bromus tectorum) 


Hardwood and Shrub Leaves 


Beech (Fagus spp,) 


Chestnut oak (Q. montana) 
Scarlet oak (Q. coccinea) 


Lauraceae~s.E. Asia 


(Lauraceae spp,) 


Litsea (Litsea spp,) 


Rhododendron (Rhododendron 
cataLJbiense) 


Madrone (Arbutus menziesii) 


Arrowleaf balsam root 


(Ba~samorhiza sagittata) 


Oregon grape (Berberis repans) 


Service berry 
(Ame~anchier spp.) 


Messmate stringybark 


(Euca~yptus ob~iqua) 


Conifer Leaves and Needles 


Western larch (Larix 
occidenta~is) 


Redwood (Sequoia sempervirens) 


Western red cedar (Thuja 
p~icata) 


Chamise (Adenostoma 
fascicu~atum) 


Douglas-Fir (Pseudotsuga 
menziesii) 


Short leaf pine (Pinus echinata) 


Lodgepole pine (Pinus contorta) 


Coulter pine (Pinus 
cou~teri) 


Ponderosa pine (Pinus ponderosa) 


Engelmann spruce (Piaea 


enge~mannii) 


Thermal 
Uiffusivity, 


(J 


2 
cm 
s 


1.90 


2.12 


1. 88 


1-.98 


2.35 


2.10 


1. 96 


2..14 


2..12 


2.10 


2.12 


1. 92 


2.03 


2.16 


1.97 


2..04 


2.37 


2.15 


2.14 


1.92 


1. 90 


2.06 


2,06 


2.14 


1. 94 


1. 96 


1. 94 


2.03 


1. 96 


1. 93 


pcL 


0.39 


0.33 


1.41 


1.11 


0.82 


1. 76 


2.60 


3.89 


5,73 


1. 05 


2.00 


3.19 


3.04 


2.18 


3.75 


3.79 


2..02 


3.16 


3.36 


5.57 


4.1 


5.79 


6.15 


5.78 


9,75 


9.05 


10.4 


9.12. 


10,6 


12.4 


Absorptivity" 
a 


Dimensionless 


0.4 
0,3 


0.5 


0.4 


0.3 


0.5 


0.5 


0.6 


0.3 


0.6 


0.8 


0.8 
0.7 


0.7 


0.6 


0.8 


0.8 
0,6 


0.7 


0.6 


0.6 


0.8 


0.7 


0.8 


0.8 


0,6 


0.8 


0.8 


0.8 


0,8 


0.8 


0.7 


(pL)1/2 
C- 


cm 


9.75 x 10-4 


1.11 
x 10-3 


2.8 
x 10- 3 


2.8 
x 10- 3 


2.7 
x 10- 3 


3.5 
x 10-3 


5,2 
x 10- 3 


6.5 
x 10- 3 


1.9 
x 10- 2 


1. 76 x 10-3 


2.5 
x 10-3 


4.0 
x 10- 3 


4. 6 
x 10-3 


4,3 
x 10-3 


3.6 
x 10- 3 


4.7 
x 10- 3 


4.7 
x 10-3 


6.3 
x 10-3 


2..9 
X 10- 3 


5.3 
x 10-3 


5.6 
x 10- 3 


7.0 
X 10-3 


5.9 
X 10- 3 


7.2 
x 10-3 


7.7 
x 10-3 


9.6 
x 10- 3 


12..2 
x 10- 3 


11.3 
x 10- 3 


13.0 
x 10-3 


11.4 
x lC- 3 


1.-32 x 10- 2 


1.8 
x 10- 2 


14-3 


-I,:=. 


(11. 
- 
I 
~ 
~ 


I 
(11, 
~ 
~' 


I 
~ 


Table 14-2. 
Comparison of Weapon- Test Data with 'Predictions Based on Correlation Model (Kerr, 1971) 


f>.1aterials 


Beech 
Coulter 
Ponderosa 
Horsehair 
Desert 
Leaves. 
Madrona, 
Pine 
Pine 
Lichen 
Cheatgrass 
j'oleedlegrass 
Wheatstraw 
Weathered 
Weathered 
Needles 
Needles 


Ignltion Thresholds, cal/cm- 2 


Yield. 
t max ' Test 
Calc. Test 
Calc. 
Test 
Calc. 
Test 
Calc. Test 
Calc. 
Test 
Calc. 
Test 
Calc. 
Test 
Calc. 
Weapon Test 
kt 
s 
Data* Values Data· 
Values** 
Data'" Values·'" 
Data* 
Values Data* 
Values 
Data* 
Values 
Data* 
Values f 
Data* 
Values§ .. 


Buster 
Easy 
31 
0.190 
3.7 
7.5 
4.6 
9.5 
II. I 
10.4 
9.6 
(4 ) 
(6) 
(7) 


Snapper 
3 
31 
0.185 
5 
5.1 
4.3 
5.2 
3.0 
8.3 
>8.9 
9.6 
3.4 t 
5.2 
6.1 to 7.3 
8.3 
7.3 to 8.9 
9.6 
7.3 to 8.9 
10.4 
9.7 
(7.9) 
(7. I) 
(#) 
(6.3) 
(9.5) 
(9.0) 
(7.4) 
(7.4) 


4 
19 
0.155 
5.3 
5,7 
6,6 
5.4 
8.2 
9.8 
5.3 
5.4 
5.3 to 6.6 
8.6 
10.8 
9.6 
8.2 to 10.8 
10.8 
9.9 
(14) 
(10.6) 
(10. I) 
(13.6) 
(11.9) 
(II. 4) 
(12.6) 


Upshot/Knothole 
4 
II 
0.103 
4.8 
9.4 
6.8 
(5.4) 


9 
26 
0.179 
6.3 to 7.7 
9.5 
8.8 
(4.4) 


10 
15 
0.138 
9.5 
9.2 
8,0 
(5.8) 


'" Estimated moisture contents (expressed as percent of dry weight) are gIven in parentheses. 


Calc~lated values are based on finest particles.. 
-2. 
._ 
t PossIbly a spuriOUS value. 
Other samples located to receIve exposures up to 7.3 cal cm 
dIdn't IgnIte. 


(N) Unknown but assumed to be 
B% for purposes of calculation. 


t Based on average thickness of needles 


§ Based on thinnest needles 


• 


:, 
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CHAPTER 15 
RESIDUAL NUCLEAR RADIATION 


• 


15-1. Gross phenomenology. 


a. Radiation emitted later than one minute after 
detonation is defined as residual. Residual radiation 
is composed of fission-product reactions and neutron 
reactions from fallout and of the neutron-induced 
gamma rays discussed in chapter 3. Fallout is the 
most significant contribution for most burst condi- 
tions of interest; it is separated into the early phase 
(which settles within 24 hr) and the delayed phase 
(which settles after 24 hr). Early fallout has military 
significance because exposure to such fallout may 
debilitate humans-they must be protected during 
those periods when the radiation intensity poses a 
threat. The design requires both dose rate and ac- 
cumulated dose (Glasstone, 1977). Neutron-induced 
soil activity (chap. 3) should be considered for burst 
elevations less than 75 ftlkto.'35. 


b. Effects of residual nuclear radiation depend on: 
- Weapon design 
- Weapon yield 
- Burst elevation 
- Range 
- Climatic conditions 
- Soil composition 


15-2. Early fallout. 


a. Air bursts. Except for the phenomenon of rain-out 
(para. 15-4), dispersion of radioactive debris from an 
air burst before it reaches the ground is usually suf- 
ficient to preclude a hazard of military significance. 
b. Surface bursts. Deposited radioactivity during 
the first day after detonation of a surface burst can 
extend for several hundred miles. Dose-rate contours 
for the environment one hour after burst over ideal 
surfaces can be prepared, as shown in figure 15-1, 
from the data presented in figures 15-2 through 
15-11. For reasonably level terrain, the values given 
in these figures should be multiplied by 0.7; rough 
and hilly terrain requires the data to be multiplied 
by a factor of 0.5 to 0.6. The effective wind is the 
average of the net wind vectors at the top and bot- 
tom of the stabilized cloud, as defined in figures 15-12 
and 15-13. The net wind vectors should be obtained 
from a plot of the wind vectors (hodograph) as a func- 
tion 
of altitude (fig. 
15-14, for 
example). 
The 
hodograph should be plotted for equal increments of 
altitude. The net wind vectors are the vectors extend- 
ing from zero elevation to the points on the hodograph 
corresponding to the top and bottom of the stabilized 
cloud. For example, the net wind vectors and effec- 
tive wind for a typical blast are superimposed on the 
hodograph of figure 15-14. 
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Figure 15-1. 
Idealized Early Fallout Dose-Rate Contour (DNA, 1972) 
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Figure 15-2. 
Downwind Distance and Fallout Radiation at t = 1 Hr After Burst as a Function ofYield, 10-Knot Effective Wind (DNA, 1972) 
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Figure 15-3. 
Downwind Distance and Fallout Radiatwn at t = 1 Hr After Burst as a Function of Yield, 2G-Knot Effective Wind (DNA, 1972) 
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Figure 15-4. 
Downwind Distance and Fallout Radiation at t = 1 Hr After Burst as a Function ofY~Id, 40-Knot Effective Wind (DNA, 1972) 
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Figure 15-5. 
Maximum Width and Fallout Radiation at t = 1 Hr Afier Burst as a Function of Yield, 10-Knot Effective Wind (DNA, 1972) 
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Figure 15-6. 
Maximum Width and Fallout Radiation at t;= 1 Hr After Burst as a Function of Yield, 2D-Knot Effective Wind (DNA, 1972) 
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Figure 15-7. 
Maximum Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 40-KnotEffective Wind (DNA, 
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Figure 15-8. 
Distance to Maximum Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 100Knot Effective Wind 
(DNA, 1972) 
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Figure 15-9. 
Distance to Maximum Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 2D-Knot Effective Wind 
(DNA,1972) 
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Figure 15-10. 
Distance to M=imum Width and Fallout Radiation at t = 1 Hr After Burst as a Function ofYield, 4Q-Knot Effective Wind 
(DNA, 1972) 
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Figure 15-11. 
Ground Zero Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield (DNA, 1972) 
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Figure 15-12. 
Height of the Stabilized Cloud Bottom as a Function of Yield (DNA, 1972) 
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Figure 15-14. 
Typical Hodograph of a Vertical Wind Structure 


c. Transition height bursts. For burst heights be- 
tween the surface and the scaled height 100 ft/kt1!3 
(30.5 m/kt1.'3), fallout contamination will be a 
variable problem. The elevations below which fallout 
should be considered are presented in figure 15-15. 
An estimate of the dose rate should be obtained by 
multiplying the data in figures 15-2 through 15-11 
by the adjustment factors presented in figure 15-16. 
For bursts in the upper one-quarter ofthe transition 
zone, the neutron-induced soH activity described in 
paragraph 3-3frepresents a significant contribution 
to the residual radiation environment that must be 
considered in addition to the fallout. 


15-3. late fallout 


a. Decay rate. The decay rate of radioactivity nor· 
malized to the quantitative data of radiation dose 
rates presented in figures 15-2 through 15-11 is sum- 
marized in figure 15-17 for any time after detonation. 
The total dosage that will be accumulated for ex- 
posure starting one hour after burst can be calculated 
from figure 15-18 and the one-hour dose rate data of 
figures 15-2 through 15-11. 
b. Mean decay rate. For practical purposes, the 
fallout radiation can be assumed to obey the t-1.2 law 
in which the dosage rates at t = 1 hour are modified 


15-14 


by the factor t-1.2 to obtain the dosage rates at an ar· 
bitrary time. 
. 


15-4. Precipitation effects. 


In inclement weather, air bursts (which usually-pro- 
duce residual radiation effects of negligible military 
significance) may cause dangerous contamination on 
the ground. The effect is called preCipitation scaveng- 
ing, or rain-out, and its importance depends on 
whether the burst occurs before or during precipita- 
tion; on the weapon yield, size and shape of the 
nuclear cloud, and radioactive decay rate; on the rain 
cloud size and type and duration of precipitation; and 
on the effects of precipitation on the contamination 
after it reaches the ground. Specific details regarding 
precipitation effects should be obtained from the DNA 
(1972) report. 


15-5. Uncertainty. 


a. The data presented can be considered only a guide 
for estimating residual contamination. The greatest 
uncertainty originates with the weather itself. 


Th~refore, the degree of surface residual radioactive 
contamination is largely dependent on precipitation 
and wind conditions during and after detonation. 
Other effects that should be considered are weapon 
yield and HOB. 
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Figure 15~15. 
Minimum Height of Burst Below Which Fallout Occurs as a Function of Yield (DNA, 1972) 
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Figure 15-16. 
Height of Burst Adjustment Factors for Various Yields (UNA. 1972) 
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Figure 15-17. 
Fission Product Decay Factors Normalized to Unity at 1 Hour After Detonation (DNA, 1972) 
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Figure 15-18, 
Normalized Dose Accumulated in a Fallout Contaminated Area from H + 1 Hour to H + 1000 Days (DNA, 1972) 
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b. The effect of yield and HOB uncertainties on pro- 
ducing the dose rate uncertainty can be computed 
from 


Q2 
_ 
2 
1 
D' - 
QF + ( 
Z 
) 4 
log -- y 
Wn 
. 


l~n + B~Og ~n - ~~Qiv + P2Q~} 


(15-1) 
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., 


• 


where 


B 
== ~1.824, 
to.824, 
W > 9.2 R2.96 kt 


TM 5-858-2 
. 
" 


'0.273, 
Z1Wo.35 ~ 337 ftlkto.35 


B 
== lo 
Z/WO.35 < 337 ftlkto.35 


y 
== 
2.11 
n 
== 
0.35 
Qz and Qw are .the uncertaiJ;lties of HOB and yield, 
respectively; and QF is associated with the uncer- 
tainty in the wind. It is suggested that radiation con- 
tours similar to that shown in figure 15-1 be con- 
structed for wind vectors corresponding to the mean 
value of the winds aloft data and for the means plus 
their corresponding standard deviations. From these 
data two radiation dose rates can be determined-one 
for the mean wind condition (call it D), and another 
for the mean plus the standard deviation wind condi- 
tion (say Do)' Therefore, compute 


D 
- D 
Q F 
== 
0 
(15-2) 
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Airblast: 


Airblast-induced Ground Shock: 


Alpha Particle: 


Azimuthal Magnetic Field: 


Beta Particle: 
Burst: 
Burst Configuration: 


Burst Elevation: 


Compton Current: 


Continuous Coverage Region: 


Crater: 
Apparent Crater: 


True Crater: 


Crater-induced Ground Shock: 


Cratering Efficiency: 


Debris: 


Decay: 


Delivery Accuracy: 
Deposition Region: 


Depth of Burst (DOB): 


Direct-induced Ground Shock: 


Dynamic Pressure: 


Ejecta: 


Ejecta Dust Cloud: 


Electromagnetic Pulse (EMP): 


Engulfment Radius: 
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The overpressure signal following the shock wave 
driven through the air by an explosion. 
Ground response caused by passage of the airblast 
over the ground surface. 
A helium atom stripped of its two electrons; a helium 
nucleus. 
A magnetic field whose field lines are circles in planes 
parallel to the ground surface. 
. A free electron. 
Explosive release of energy, explosion. 
Description of the location of the explosive center of 
growth relative to the ground plane, and to the 
target or observation point. 
Height of burst; distance from the ground surface up 
to the explosive CG. 
The current resulting from the collision of gamma 
rays with air molecules and the outward displace- 
ment of electrons detached by the collisions. 
The area where sufficient ejecta is deposited to pro- 
vide complete coverage of the ground surface. 


The open hole left after a cratering explosion, smaller 
than the true crater due to partial filling by 
fallback. 
The excavation or hole defining the limit of dissocia- 
tion (separation) of the natural site material. 


Late time stress and particle motion signals at- 
tributed to crater formation. 
The ratio of crater volume to weapon yield for crater- 
ing explosions. 
Ejecta deposited on the ground surface; wreckage 
from damaged structures or vegetation. 
Decrease in amplitude; also, radioactive decay, the 
natural fission process which leads eventually to 
the formation of stable (nonradioactive) substances. 
The CEP of a weapon system or delivery system. 
The source region for EMP; the region where the 
primary gamma ray interactions with the at- 
mosphere occur. 
The distance from the ground surface down to the ex- 
plosive center of gravity; the depth of burial. 
Ground response resulting from the coupling of 
energy in the ground in the vicinity of the crater. 
One-halfthe product of the density and square of the 
particle velocity of air behind the shock front. 
Material thrown out ofthe region surrounding ground 
zero during crater formation. 
The portion of ejecta particles forming a dust cloud 
subject to prevailing wind forces. 
The "radio flash" 
on broad band electric and 
associated magnetic fields resulting from nuclear 
explosions. 
The maximum radius of the fireball. 
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Enhanced Weapon: 


Equivalent Ejecta Thickness: 


Equivalent Yield: 


Fallback: 


Fallout: 


Fireball: 


Fission: 


Fluence Spectrum: 


Free Field: 


Fuel Density: 


Fusion: 


Gamma Rays: 


Ground Shock: 


Ground-Surface Dust Cloud: 


Ground Zero (GZ): 


Height of Burst (HOB): 


"High Yield" Weapon: 


Hydrodynamic: 


Impedance: 


Impulse: 


Laydown Pattern: 


Load: 


A-2 


A weapon designed to enhance some asp~ct of bomb 
performance, e.g., the neutron bomb or an earth 
penetrating reentry vehicle/warhead, which respec- 
tively enhance neutron radiation at the expense of 
blast and ground shock and ground shock at the ex- 
pense of airblast and radiation. 
An ejecta depth based on uniform azimuthal distribu- 
tion which ignores the normally observed ray pat- 
tern and azimuthal variation. 
The yield of a weapon detonated under different con- 
ditions that would produce the same effect at the 
point of observation. 
Material lofted by the cratering process which does 
not escape but drops back into the crater. 
Radioactive bomb case debris and dust which settles 
out of the dust cloud. 
The high temperature, high pressure expanding 
plasma core produced by a nuclear explosion. 
The splitting ofunstable heavy element atoms accom- 
panied by the release of energy. 
A spectrum relating radiation intensity (fluence) to 
energy level. 
Locations in the air or ground unaffected by the 
presence of man-made structures. 
A measure of the availability of fuel for evaluating 
firestorm potential, i.e., weight of fuel per unit 
ground surface area. 
. . 
The combination of light nuclei into heavier nuclei 
accompanied by the release of energy; the basis of 
thermonuclear (hydrogen) bombs. 
Radiation similar to visible light or x-rays but with 
very short wavelengths; hard x-rays. 
The ground stress and ground motion resulting from 
direct energy deposition and/or airblast' loading 
from an explosion. 
The portion of the low:-lying dust cloud comprising 
dust that originates from the ground surface (out- 
side the crater). 
The point on the ground surface directly above or 
below the burst point. 
The elevation of the explosion's center of gravity 
above the ground surface. 
A general class of nuclear weapons that release a 
major portion of their energy as radiation; high 
yield to mass ratio devices. 
Fluid like; pertaining to the expansion of the hot 
plasma core during the early phases of a nuclear 
explosion. 
. 
The 
product 
of material density, e, 
and 
a 
characteristic wave speed, c. 
The integral of the overpressure, usually for the 
positive phase duration. 
An attack strategy utilizing the delivery of a number 
of weapons for subsequent simultaneous detona- 
tion. Delivery may be in an array to attack point, 
line, or areal targets. 
The integral ofthe product ofthe dust concentration 
and the air-intake-volume rate over a lapsed time. 
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Low-lying Dust Cloud: 


"Low Yield" Weapons: 


Mach Reflection Region: 


Mach Stem: 


Million Electron Volts (MeV): 
Multiple Burst Attack: 


.Laydown Pattern: 


Sequential Attack: 


"Nail Driving:" 


Neutrons: 


Neutron-induced Gamma Activity: 


Nuclear Radiation: 


Nuclear Weapons: 


Outrunning: 


Overpressure: 


Pattern Bursts: 


Plastic Zone: 


Polar Electric Field: 


Positive Phase (duration): 


Precipitation Scavenging: 


Precursor (shock): 


Rad: 


Radial Electric Field: 


Radiant Exposure: 
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The dust cloud that remains in contact with the 
ground surface. 
A general class ofweapons which release a major por- 
tion of their energy as kinetic energy of bomb case 
debris; low yield to mass ratio devices. 
The ground surface region below a low-altitude burst 
where reflected shock waves overtake and merge 
with the incident shock waves to form a Mach stem. 
The nearly vertical shock front formed by the merg- 
ing of reflected and inCident shock waves in the 
Mach reflection region. 
Temperature or energy unit for high energy particles. 


An attack using simultaneous detonation of an array 
of previously delivered weapons. 
Repeated attack of a target, especially when crater- 
ing is a primary kill mechanism, or can be used to 
dig out a buried target; also known as Nail Driv- 
ing attack. 
A sequential multiple-burst attack designed to dig out 
and kill a buried target. 
A particle of neutral electrical charge found in most 
atomic nuclei; atomic weight::::: 1. 
Residual radioactivity induced in soil by neutron bom- 
bardment during a nuclear explosion. 
The alpha, beta, and gamma radiation resulting from 
a nuclear explosion. Commonly divided into intitial 
or prompt radiation (0 < t < 1 min) and residual 
radiation (t > 1 min). 
Weapons whose explosive energy is the result offis- 
sion or fusion processes. Includes both atomic (fis- 
sion) and thermonuclear (fission/fusion) weapons. 
The arrival of ground signals propagated through 
underlying higher velocity layers at the ground sur- 
face before the arrival of the airblast shock front. 
The airblast pressure (in excess of· the ambient 
pressure) resulting from an explosion. 
Multiple burst attacks utilizing simultaneous detona- 
tion of an array of weapons. 
The zone of soil/rock surrounding a crater in which 
significant plastic deformation has occurred. 
An electric field normal to the earth's surface; a part 
of the EMP signal associated with an azimuthal 
magnetic field. 
The length of time the airblast signal shows a positive 
overpressure. 
Removal of radioactive dust from the dust cloud of a 
nuclear explosion by rain; rainout. 
A shock wave running ahead of the normal shock 
front, caused by radiative heating of thermally 
nonideal surfaces. 
A unit of absorbed radiation dose equivalent to 100 
ergs of radiation per gram of absorbing material or 
tissue. 
An electric field directed radially from the center (of 
a nuclear explosion). 
The amount of thermal radiation received per unit 
area (cal/cm2). 
A-3 
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Radiated Pulse: 


Radiation Dose: 
Rainout: 


Reflectance: 


Regular Reflection Region: 


Rupture Zone: 


Secondary Gamma Radiation: 


Sequential Detonation: 


Shake: 
Shock Wave: 


Slant Range: 


Spontaneous Ignition: 


Superseismic: 


Temperature Inversion: 


Thermal Efficiency: 


Thermal Radiation: 


Thermally Ideal Surfaces: 


Thermally Nonideal Surfaces: 
Transition Height: 


Transmittance: 


Triple Point: 


Waveform: 
Weapon Design: 


X-ray: 


Yield: 


A-4 


That portion ofthe EMP disturbance radiated by the 
induced electric currents. 
Total quantity of ionizing (nuclear) radiation. 
Removal of radioactive dust from a nuclear dust cloud 
by precipitation; precipitation scavenging; rain- 
induced fallout. 
Ratio ofthe radiant energy reaching a surface facing 
the fireball to that which would reach the same sur- 
face in the absence of clouds, the ground, and other 
reflective surfaces. 
The region on the ground surface below a low altitude 
burst where the reflected wave does not overtake 
and merge with the incident shock; the region 
within a circle ofradius equal to the height ofburst, 
centered at ground zero. 
A region surrounding the true crater in which crack- 
ing and crushing occur. It, in turn; is surrounded 
by the plastic zone. 
Gamma radiation produced by neutron interaction 
with materials outside the nuclear weapon. 
An attack strategy in which a target is repeatedly at- 
tacked, such as the "nail driving" attack to dig out 
and kill a. buried target. 
10-8 second 
The air-blast wave propagated by an explosion, 
characterized by 
discontinuous 
rise 
to 
peak 
pressure. 
Spherical range from the explosive center of gravity 
to the point of interest. 
Ignition as the result ofradiant exposure to the ther- 
mal pulse from a nuclear explosion. 
Ground signals moving faster than the airblast shock 
front. 
Layer of air in which the temperature increases with 
height. 
Ratio of the thermal energy radiated by the fireball 
to the total energy release or yield of the weapon. 
Radiation from the fireball in the ultra-violet, visi- 
ble, infrared spectral bands. 
Perfectly reflective, nonabsorbing surface, e.g., water, 
Ice. 
Absorbtive, nonreflective surface. 
The minimum burst elevation for which dust pickup 
and fallout do not occur. 
Ratio ofthe radiant energy reaching a surface facing 
the fireball to that which would reach the surface 
through a perfectly transparent atmosphere. 
The point at which the Mach stem and incident 
reflected shocks intersect. 
The time-dependent amplitude history of the signal. 
Weapon construction factors influencing the perfor- 
mance of a nuclear weapon; cf enhanced weapon. 
A gamma-like radiation not associated with nuclear 
explosions. 
Explosive. energy released by a nuclear explosion, 
commonly expressed in terms of the weight of an 
energy equivalent TNT charge. 
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